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INTRODUCTION 


The  work  under  this  contract  essentially  consisted  of  adapting  the  existing 
Grumman  Specific  Linear  Optimal  Control  Program  to  the  hydrofoil  case  and  of  illus¬ 
trating  the  utllitj^  of  that  adaptation  by  application  to  the  AG  (EH)  hydrofoil  craft.  The 
computer  program  itself,  the  adaptations  required,  and  the.  application  to  the  AG  (EH) 
are  all  discussed  in  Volume  II  of  this  report.  This  volume  is  limited  to  consideration 
of  the  foil  and  craft  hydrodynamic  characteristics  required  by  the  Specific  Linear 
Optimal  Control  Program. 

The  hydrodynamic  study  was  heavily  influenced  by  O'Neill's  demonstration  in 
Reference  5-9  that  unsteady  loading  effects  can  be  significant  to  hydrofoil  craft  of 
larger  sizes.  Much  of  this  volume  is  devoted  to  evaluation  of  the  magnitude  and 
confidence  level  for  unsteady  loading  effects.  Another  major  area  of  consideration 
was  the  limitation  imposed  by  cavitation  on  the  linear  range  for  foil  and  flap  lift. 

The  hydrod5mamic  and  control  studies  were  of  comparable  magnitude  and  pro¬ 
ceeded  concurrently.  Consequently,  the  hydrodynamic  characteristics  employed  in 
the  control  study  were  preliminary  estimates  and  are  not  consistent  with  the  hydro- 
dynamic  conclusions.  Similarly,  the  hydrodynamic  characteristics  are  carried  in 
full  detail  because  the  control  optimization  program  was  not  available  for  the  evaluation 
of  reductions  in  form  for  those  characteristics. 

Because  Volumes  I  and  11  of  this  report  concern  distinct  disciplinary  fields,  no 
attempt  has  been  made  to  employ  a  common  nomenclature  for  the  two  volumes. 


CONCLUSIONS 


NOTE:  All  conclusions  conerning  the  flap  are  reserved  for  Appendix  E. 
Steady-State  Loadings 


•  Theory,  towing  tank  measurements,  and  prototype  neasurements  for  the  AG(EH) 
forward  foil  hinge  moment  are  contradictory.  The  experimental  measurements 
are  inconclusive  with  respect  to  the  aerodynamic  center  and  residual  hinge 
moment  individually  and  reflect  doubt  upon  the  confidence  level  for  foil  hinge 
moment  generally.  Further  model  testing  with  symmetric  and  cambered  sections 
on  identical,  simpler  planforms  would  be  desirable. 

Unsteady  Loading 

•  For  foil  motions,  the  unsteady  loading  responses  of  this  report  are  basically 
classic  aerodynamic  reopv/..res  though  some  of  Theodorsen's  ’T”  coefficients 
have  been  replaced  by  practical  values,  Jones'  3-dimensional  foil  motion  function 
has  been  employed  in  place  of  Theodorsen's  function,  and  the  S-dimensional 
hydrodynamic  lift  curve  slope  nas  been  employed.  The  experimental  data 
examined  for  this  report  do  not  require  any  more  detailed  accountability  for  the 
free  surface 

•  The  confidence  level  for  lift  and  foil  hinge  moment  response  to  foil  motion  is 
generally  adequate,  subject  to  SLOCOP  studies  of  the  sensiti\ity  of  the  autopilot 
design  to  these  responses,  though  experimental  evidence  is  inadequate.  The 
principal  limitation  upon  the  prediction  of  the  foil  moment  response  is  the  predic¬ 
tion  of  the  steady-state  foil  moment. 

•  The  celeriiy  of  gravity  waves  presents  a  significant  distinction  between  the  aero¬ 
dynamic  and  hydrod5mamic  orbital  motion  (gust)  cases.  Correlation  of  theory  and 
experiment  is  poor/fair  for  lift  and  inadequate  for  foil  hinge  moment.  SLOCOP 
studies  under  way  are  intended  to  establish  the  accuracy  requirement  in  this  area. 

Cavitation  Characteristics 

•  Optimum  cavitation  characteristics  for  the  AG  (EH)  forward  foil  with  a  20%  chord 
flap  result  when  the  incidence  is  fixed  at  an  angle  producing  a  lift  coefficient  of 

.  167  (1185  psf  foil  loading  at  50  knots).  Varying  the  foil  incidence  with  speed  to 
trim  the  flap  at  Ig  is  disadvantageous  for  cavitation 

•  Conclusions  concerning  flap  effectiveness  and  cavitation  effect  boundaries  are 
reserved  for  Appendix  E. 
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SYMBOLS 

NOTES;  1.  All  units  in  ft. /# /sec. /rad.  unless  otherwise  noted, 

2.  Parenthesis  distinguish  force  or  moment  source  from  slope;  e.g. , 

(C  )  is  C  due  to  incidence. 

Li  L 

3.  Primes  indicate  exposed  foil,  or  phase  referenced  to  orbital  angle  at 
half  chorvl,  or  characteristic  measured  in  plane  perpendicular  to 
quarter-chord  line. 

4.  Refer  to  Table  5-1  for  summary  of  unsteady  load  symbols. 

5.  For  dual  signs  (+,  +);  upper  sign  refers  to  upper  surface,  lower  sign 
to  lower  surface. 

6.  In  general,  unsteady  characteristics  are  distiaguished  from  steady-state 

characteristics  by  symbol  identifying  the  unsteady  motion;  e.g. ,  C  is 

^h 

C  response  to  heave  motion;  is  normalized  response,  C_  /C_  . 

L  n  L,  L 

n 

C  is  ?  scalar,  C  is  a  vector. 

L  Lj, 

A  Aspect  ratio 

a.c.  Aerodynamic  center 

B  Buoyancy 


Chord.  May  be  MAC,  MAC’,  MGC,  etc. 
Buoyancy  coefficient,  B/qS 
Drag  coefficient,  D/qS 
Flap  chord 

Incidence  hinge  moment  coefficient,  H/qSc 


Flap  hinge  moment  coefficient,  h/qS^c^ 


dCR/dC^ 


dCj^/dl 


1 


3 


C(K) 


'■'I'  ''L'l 


dC„/d» 

n 

3-D  dC. /dff 
n 

2- D  or  3-D  dC^/da 

h 

3- D  dC^/d6 

h 

2-D  or  3-D  dC^/d6 
n 

Jones’  3-D  foil  motion  function.  2-D  form  is  virtually  identical 
with  Theodorsen's  function.  See  Eq.  (3.4) 

Lift  coefficient,  L/qS 

O’Neill’s  C, 

""i 

Section  lift  coefficient 

Additional  (type)  lift  contribution  of  flap 


Basic  (type)  lift  contribution  of  flap 

Total  lift  contributi'/n  of  flap,  c  +c- 

a  b 

Incidence  lift  curve  slope,  dC  /di, 

L 

Design  lift  coefficient  (camber) 

Effective  design  lift  coefficient,  c 

0=0 

Total  lift  less  flap  lift 
Zero  flap  C 

L 

Pattison’s  C 
Section  lift  curve  slope 
Foil  lift  curve  slope 
Flap  lift  oiirve  slope 

Max.  or  min.  c^  /C^  on  foil  for  flap  lift 

Max.  or  min.  c  /C  on  foil  for  Incidence  lift 
X  L 

Max.  or  min.  c./C  on  foil  for  pitch  lift 
X  L 
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,  for  design  lift  coeff. 


a«c. 


a.  c. 


eff. 


ac. 


B 


"M. 


■'M 


a 


■'M 


a 


a.c. 


OL 


"M 


OM 


w 

D 

d 

da/d6 


n 


G(k) 


g 

H 

h 


Moment  coefficient  about  a.  c. ,  C,, 

M 

Actual  ,  see  Eq.  (4.4.1) 
ac 

C,,  due  to  flap 

a.c. 

Moment  coefficient  about  C.G.  due  to  buoyancy,  M_/qS-^ 

B 

Pattison’s  C 

Moment  coefficient  about  C.G.  due  tc  craft  pitch,  M/qSA  (not 
slope) 

Moment  coefficient  slope  about  a.c.  in  orbital  motion 


Moment  coefficient  abe-t  r.G.  due  to  residual  lift 

Moment  coefficient  about  C.G.  due  to  residual  moment 

Wave  celerity, v/Tg^  rr 
Pod  diameter 
Depth 

Flap  effectiveness 
Froude  number,  vA/gJ 

Normalized  lift  response  to  orbital  motion,  C  /C  =  C  /C  , 

L  L  L  L 
o  o  w 
w 

or  centroid  for  flap  basic  load  distribution  referenrad  to  c/4  and 
expressed  as  fraction  of  chord 

Jones’  gust  response  function  which  is  C,  /C,  referenced  in 

■  JLi  ij 

w 

in  phase  to  orbital  angle  at  L.  E. .  Virtually  identical  with 
Sears’  fimction  except  that  Sears’  function  is  referenced  to 
orbital  angle  at  mid-chord.  See  Eq,  (3.5).  Note  that  Jones’ 
function  is  3-D  and  Sears’  function  is  2-D. 

Acceleration  of  gravity 

Distance  between  L.  E.  and  incidence  hinge,  or  normalized 

incidence  hinge  moment  response,  C„  /C„ 

n  H 
a  a 

V 

Distance  between  L.  E.  and  flap  hinge,  or  heave,  +  up,  or  occa- 
sinally  depth 
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(-h/v) 

ht 

I 

i 

J(k) 

k 


k’ 

^iV(3)’  ^B(W) 

^7(3)’  ^*B(W) 
L 

I 


"M 


s 

M 

MAC 

MAC’ 

MBC 

Mf 

MGC 


M 

M, 


H 


Angle  of  attack  due  to  heave  motion 
Wave  height,  crest  to  trough 
Normalized  Incidence  moment  response, 


K 


i 


Incidence  angle,  +  nose  up,  or  imaginary  part  of 
Bessel  function 

C.  /2  TT ,  or  section  area/tc,  or  normalized  inverse  wavelength, 

Xf 

o 


k(l  +  C^/V)  =  (dc/2  V 


Wing/body  interaction  fac;,ors  from  Chart  1,  Reference  4-1 
Lift,  +  up 

Foil  base,  distance  between  fwd  and  aft  foil  hinges 
Steele’s  L 

w 

Distance  between  pod  L.  E.  and  pod  a.  c.  (c.  p. ) 

Distance  between  pod  L.E.  and  shoulder  (max.  dia.) 

Moment,  +  nose  up 

Mean  aerodynamic  chord 

Mean  aerod3Tiamic  chord  for  exposed  foil 

Mean  buoyant  chord 

Moment  due  to  flap  deflection,  +  nose  up 
Mean  geometric  chord 
Moment  about  incidence  hinge,  +  nose  up 
Moment  about  flap  hinge,  +  nose  up 
Local  pressure 

Atmospheric  pressure,  2116  psf 
P  -  P  =  P  *  +  pgd  -  P 
Ambient  pressure,  ^  Pgd 
Vapor  pressure,  72  psf 
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warn 
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*s 

v/V 

Av/V 

Av  /V  ' 
a 


(Av/V), 


(W/S)^ 

jasov) 


Dynamic  pressure,  |  p  V 
Area,  or  distance,  or  total  foil  area 
Pressure  coefficient,  1  +  (P^  -  P)/q 
Exposed  area 

2  / 

Nose  cross-sectioh  area,  tt  D  /4 
Theodorsen's  coefficients.  See  Table  3-1 
Thickness 

Speed,  +  for  craft  fwd 
Speed,  knots 

Velocity  of  encounter,  V  +  C 

K 

Volume  of  pod  nose 

Local  velocity  increment  due  to  thickness 

Local  velocity  increment  due  to  camber 

Local  velocity  Increment  due  to  angle  of  attack  for  unit 


Local  velocity  increment  due  to  flap  basic  loading  for  unit  6, 
=  1  of  Allen’s  P.  /C 


Weight 

Orbital  velocity,  ^  y/2  ngA 
Foil  loading 

Foil  loading  for  zero  flap 

Distance  from  exposed  foil  root  chord  L.  E.  to  body  lift  c.  p.  (a.c. ) 

Angle  of  attack  in  general,  0,  i,  (-  h/V),  w  /V.  +  nose  up 

8 

Flap  deflection,  +  nose  up  . 

A  flap  basic  (type)  load  parameter,  C__  /G  c ,  ^ 

n.  I  do 
6  o 

Span  station,  or  normalized  flap  moment  response,  C,  /C, 

h  h 


%  “ 


Pitch  an^e,  +  nose  up 


Sweep  angle,  quarter-chord  unless  otherwise  specified,  or 

normalized  lift  response,  C  /C_ 

L  L 

“u  “ 
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X 


X 


h 


a 


P 

a 

’^UT 


U) 


a 

B 

C 

c/4 

= 

F 

f 

foil 

H 

h 

i 

L,E. 
^  max 


Taper  ratio  or  wavelength 

Normalized  flap  moment  response  to  heave  motion, 


\/(- v)/\ 

(c  /C  )./(c  /C  )  -  1,  a  cavitation  parameter 

^  Li  1  L  L  r 

(c  /C_  )  /(cVc  )  -  1,  a  cavitation  parameter 

Z  L  O  Z  L  r 

Density,  1.  9905  slugs/ft.  ^  (#  sec.  ^/ft. 
Cavitation  number  (P^  -  P^^)/q 
PatUson's  phase  for  C 


V  A  V  — 

A  cavitation  parameter,  -  +  +  C 

V  V  z 


A  V 


V 


eff 


Reduced  frequency,  2  tt  V  A 

6 


2V 

k’,  or  a  cavitation  parameter 

0 


SUBSCRIPTS 

H 

About  Theodor  sen's  point  a  where  a  =  2  —  -  1 

L' 

Body,  or  due  to  buoyancy 
Circulatory  component 
Referenced  to  quarter-chord 
For  =  0 

For  foil  motion,  or  for  flap  motion 
Flap 

of  or  for  the  foil 
For  foil  hinge  moment 

For  flap  hinge  moment,  or  for  unsteady  heave 

Due  to  steady  or  unsteady  incidence 

For  incidence,  due  to  heave 

Leading  edge 

Maximum 
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Non-circulatory  (inertial)  component 

Orbital  motion,  or  residual  (value  for  o  =  i  =  6  =  0) 

At  zero  lift 

For  unsteady  pitch 

Root 

Tip 

Steady-state 

For  unsteady  foil  motion  generally,  represents  P  and/or  h 

For  orbital  motion 

Due  to  flap  angle 

For  unsteady  pitch 

Fwd 

Aft 

At  infinite  depth  (aerodynamic) 


9 


f 

I 


it  I 


i  f 


i\ 


'  f 


li 


Mtta 


1.  REFERENCE  CRAFT  CONFIGURATION 

The  AG(EH)  (’’Plainvlew”)  has  been  employed  as  the  model  throughout  this 
report  to  illustrate  the  generalized  hydrodynamic  reladonsliip,  A  craft  weight  of  320 
tons  with  a  90/10  distribution  is  assumed  throughout  this  volume  of  the  report.  The 
foil  position  references  are  the  foil  incidence  hinge  lines  and  the  craft  geometry  is 
shown  cn  Figure  1-1.  Figoi^es  1-2  and  1-3  present  the  foil  geometries. 

The  net  foil  lift  and  moment  characteristics  of  this  report  are  derived  by  the 
methods  of  Reference  1-1.  Reference  1-1  refers  to  the  gross  lift  and  moment  charac¬ 
teristics  of  die  complete  foil,  neglecting  the  pod,  and  of  the  ’’exposed  foil",  which  is 
a  hypothetical  foil  consisting  of  le  two  oxt^osed  panels  of  the  actual  configuration 
joined  together  at  their  root  choi'ds.  The  gross  lift  curve  slopes  for  the  complete  and 
exposed  foils  for  the  AG  (EH)  forward  and  aft  configurations  are  derived  in  Tables  1-1 
and  1-2  by  the  methods  of  Reference  1-2.  The  complete  foil  and  exposed  foil  mean 
aerodynamic  chords  and  aerodynamic  centers,  also  derived  from  Reference  1-2,  are 
shoi^ii  on  Figures  1-2  and  1-3. 

1.1  REFERENCES 

1-1  Pitts,  Nielsen,  a^u  ivaattare,  "Lift  and  Center  of  Pressure  of  Wing- 
Body-Tail  Combinations  at  Subsonic,  Transonic,  and  Supersonic 
Speeds",  NACA  Report  1307,  1959. 

1-2  DeYoung  and  Harper,  "Theoretical  Symmetric  Span  Loading  at  Sub¬ 
sonic  Speeds  for  Wings  Having  Arbitrary  Plan  Form",  NACA  Report 
921,  1948. 

1-3  Lindsey,  Stevenson,  and  Daley,  "Aerodynamic  Characteristics  of  24 
NACA  16  -  Series  Airfoils  at  Mach  Numbers  Between  0. 3  and  0.  8," 
NACA  TN  1546,  September  1948. 

1-4  Richter  and  Palmer,  "AG (EH)  Hydrofoil  Reserach  Ship  Hydrodynamic 
Report  of  Model  Test  Program,"  Grumman  Aerospace  Corp.  Report 
M23.  77,  8/21/62,  Rev.  1/15/63. 
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TABLE  1-1  GROSS  LIFT  CURVE  SLOPES-FWD 


CHARACTERISTIC 

— 

REFERENCE 

COMPLETE  FOIL 

EXPOSED  FOIL 

Aspect  Ratio,  A 

3 

2.805 

Taper  ratio,  X 

.3 

.329 

c/4  sweep.  A,  deg. 

35.2 

35.2 

Section 

16  -  (.353)  08 

16  -  (.  353)  08 

Section  k  =  c^^  /2tt 
a 

NOTE  1. 

.9 

.9 

aA 

3.34 

3.12 

0  A 

a 

Fig.  4,  Ref.  1-2 

.0556 

.  0540  • 

Lift  curve  slope  C 

a 

2.87’ 

2.79 

Ref.  2  C  Prediction 
L 

a 

3.01 

2.92 

NOTES;  1.  From  anunpublished  compilation  of  measured  characteristics  of  *  ^ 

Reference  1-3.  Section  lift  curve  slope  was  measured  at  Reynolds 

6  f  1 

number  of  .  85  x  10  and  Mach  number  of  ,  3,  Slope  is  function  of  j  j 

section  thickness  ratio  only  and  this  table  is  valid  for  model  and 

prototype.  pi 

2.  Predictions  of  this  report  differ  from  those  of  Reference  1-4  in  use 
of  .  9k  rather  than  1.  Ok  and  in  interpolation  in  various  charts 

employed.  j  J 

n 

n 

□ 
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Figure  1-1.  Craft  Geometry 
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Figure  1-2.  Forward  Foil  Geometry 
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Figure  1-3.  Aft  Foil  Geometry 


2.  NOMENCLATURE 


The  principal  items  of  romenclature  employed  in  this  volume  are  illustrated  on 
Figure  2-1;  <>11  symbols  are  defined  in  the  list  of  symbols. 

The  sign  sense  of  this  volume  differs  from  Theodorsen  for  lift  and  heave  which 
are  both  taken  as  positive  up  in  this  volume.  Note  that  the  half-chord  is  never  employed 
in  this  volume  and  that  the  symbol  ”b"  is  employed  for  foil  span. 

In  discussions  concerning  the  total  craft  configuration,  subscripts  are  employed 
to  distinguish  forward  and  aft  characteristics;  e.g. , 

C  =  fwd  lift  coefficient 

C  =  aft  lift  coefficient 

^2 

.  C  =  craft  lift  coefficient  =  C  +  C 
L  S  b  Lig 

In  discussions  of  a  single  foil,  the  foil  location  not  identified  by  subscript. 

Exponential  representation  for  vectors  is  avoided  in  this  volume,  vectors  being 
presented  as  complex  numbers  and  phase  shifts  symbolized  by  [  +  <t>.  All  angles  are  in 
radians  unless  symbolized  0  . 

The  lift  curve  slope  for  incidence  motions  on  the  side  of  a  fixed  pod  differs 

from  that  for  a  foil  rigidly  attached  to  the  pod.  The  two  cases  are  distinguiph<^d 

verbally  here  as  "pitch  lift"  and  "incidence  lift"  and  are  distinguished  symbolically 

as  C  and  C  .  In  seas  "Incidence  lift  control"  systems  (including  flap  lift  control 
L  L, 

a  i 

systems)  respond  to  craft  and  orbital  motions  on  the  pitch  lift  curve  slope,  it  is  only 
incidence  motions  (Including  flap  motions)  which  produce  lift  at  the  incidence  lift 
curve  slope  rate. 
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Hous-i  t  vAiUfc  m  >Homi 
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c/1,  /  (f  :  f\i’o  a^rp. 


Figure  2-1.  Nomenclature 


3.  GENERAL  FOIL  LIFT,  PITCHING  MOMENT,  AND  FLAP  HINGE  MOMENT 


The  inotantaneous  foil  lift.  Incidence  hinge  moment,  and  flap  hinge  moment  are 
the  sum  of  foil  motion  and  orbital  velocity  Induc'^d  components  superimposed  on  a 
steady-state  or  trim  value. 

t 

Theodorsen  presents  the  foil  motion  components  for  the  two  dimensional  foil  In 
Reference  3-1.  Theodorsfn's  equations  are  accommodated  to  the  three  dimensional 
foil  In  this  report  by  the  following  modifications: 

•  Theodorsen's  function,  C(k),  is  employed  in  Jones'  three  dimensional 
form  of  Reference  3-2 

•  Theodorsen’s  lift  curve  and  flap  hinge  moment  slopes  (vdiich  are  2-dimen- 
slonal,  thin  airfoil  theory)  are  replaced  by  a  practical,  3-dlmensional 
estimate.  The  estimates  of  Reference  3-3  are  employed  in  this  report 

•  Theodorsen's  flap  effectiveness,  T^^/n  (which  is  2-dlmensional,  thin  air- 
-  foil  theory),  is  identified  as  dtaf/d6  for  evaluation  by  practical  methods 

•  Inertial  terms  which  are  functions  of  incidence  or  flap  motion  are  reduced 
by  the  ratio  of  exposed  to  total  foil  area 

•  Jones'  gust  function,  G(k),  is  added  to  account  for  orbital  motion 

•  A  trim,  steady- state  value  is  added. 


The  result  of  these  modifications  to  Theodorsen's  Eqs.  (XVni)  -  (XX)  of 
Reference  3-1  may  be  written: 


iJ 

-1-k' 

h  +  fl  +  1  2(1 

(  ’"'-A 

e  +  ^  iU 

L  TT  n  , 

J 

c 

L  '4 

c/  J 

\ 

a 

Oi 


r*L  liiJ 
[d6  ^2  n  J 


w 

-I 


H 


Cl  C(k)  +  CL  -^G(k)  2^k+CL 
o  a  L_  trim 


=c  +c  +c  +c 

Sc_  ^C_  '"w  ^trim 

F  F 


(3.1) 
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where  the  reference  phase  is  the  orbital  angle  at  the  incidence  hinge. 

M. 


C„  = 


H 


2t7  S’  MAC'/H 


H  q  S  MAC  c  S  MAC 


(!-l) 


k'^  h 


^  S’  MAC  . .  Ir  1  /H  l\hi,  .1/3  H\) 
^  S  MAC  ^  ([32  U  ■  2/  “^2(4"  c/| 

.  iL £  mm:  (Zi  ^  ^  +  f !i  ^  o(h  .  H\Iilk.2 

2  S  MAC  (tt  d6  L  TT  *Vc  c/  nJ 

,  2(!!.  H)Il  .llu],..}  6 


-  i'lc 

V  c 

4/^L„ 

\  e 

- 

4  / 

^  Uj_^  \  , 

w  w  »■'  trim 


=  C„  c  +  c  +  c 

"nC^  "c„  "0  ».rto 

r  Jb 


where  two  important  observations  are  in  order: 


(3.2) 


•  Eq.  (3. 2)  yields  the  hinge  moment  presented  to  the  incidence  control  system; 
in  the  craft  equations  of  motion  (S’/S)  X  (MGC’/C)  factor  becomes  unit  in  the 
non-oirciilatory  heave  and  pitch  terms 

•  The  orbital  motion  term  is  doubtful. 

The  reference  phase  is  the  orbital  angle  at  the  hinge 
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+  C  <-  i  -  k'  h  + 

[i  + 1 2(7  - 

fi)k 

li  iik'l 

6 

h  1  c 

or  \ 

L 

c/ 

J  1^6 

2  n  J 

w 

.Ch  fOBzfK-ik.c 

or  u— trim 


=  c  +c  +c  +c 

^C_  V  '^trlm 


where  the  orbital  motion  term  is  doubtful  and  the  reference  phase  is  the  orbital  angle 
at  the  incidence  hinge. 

Theodorsen’s  "T"  coefficients  are  tabulated  in  Table  3-1  for  certain  flap  chords. 

Jones'  foil  motion  function  as  given  by  Eq.  (39)  and  by  the  inverse  wavelength 
transform  of  E^q.  (30)  of  Reference  3-2  may  be  written 

i  . 165  k  i  . 335  k 
^  “  .045  +  ik  “  .3+  ik 


C<k>A  =  6  =  '  - 


i  .361k 
.381  +  i  k 


V  ,  i  .  283  k  ^ 

.  =  3  ■  ^  ■  .54  +  ik 

where  the  2-dimen8ional  equation  is  Theodorsen's  function. 

Jones'  gust  functions  as  given  by  the  inverse  wavelength  transform  of  Eq. 
of  Reference  3-2  may  be  written 


=  1  - 


i  .  236  k 


i.513k  i.lTlk 


,058  +  i  k  .364  +  i  k  2.42  +  i  k 


=  1  i  >448  k  i  .  272  k  i  .  193  k 

•  =  6  ^".29+ik“.725  +  ik”3  +  ik 


i  . 679  k  i  .227  k 
=  3  .558  +  ik  '  3.2  +  ik 


where  the  2-dimensional  equation  is  Stark's  function  and  Bisplinghoff's  Eq.  (5-376), 
Reference  3-4,  except  for  phase  reference;  Jones  references  phase  to  the  orbital 
angle  at  the  leading  edge  while  Sears  and  Bisplinghoff  reference  the  orbital  angle  at 
the  half-chord. 
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k  is  the  reduced  inverse  v.^velength,  defined  as 


■  X/c 

k’  is  the  reduced  frequency,  defined  as 


(3.6) 


(3.7) 


Heave,  pitch,  flap  deflection,  and  orbital  angle  are  vectors.  The  orbital  angle 
at  the  incidence  hinge  is  the  reference  phase.  The  lift  curve  slope  is  also  a  vectorial 
function  of  depth  but  the  non-linearity  introduced  by  the  lift  curve  slope  is  not  explored 
in  this  voliune. 


The  components  of  Eqs.  (3. 1)  -  (3.  3)  are  considered  in  detail  in  the  three 
following  sections. 
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inAIAA  Selected  Reprints,  "Aerodynamic  Flutter,"  I.  E.  Garrick, 
Editor,  March  1969. 

3-2  R.  T.  Jones,  "The  Unsteady  Lift  of  a  Wing  of  Finite  Aspect  Ratio, " 

NACA  Report  681,  1939. 
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4.  TRIM  CHARACTERISTICS 


4. 1  LIFT  CURVE  SLOPE 


Reference  4-1  is  employed  to  estimate  the  lift  and  pitching  moment  characteris¬ 
tics  of  the  foil/pod  combination.  Theoretically  this  reference  is  valid  only  for  blunt- 
based  bocdes  but  experience  has  shown  that  the  pod  does  contribute  lift  in  excess  of 
that  corresponding  to  the  shrouded  foil  area  and  that  excess  cannot  be  accounted  for 
by  the  potential  characteristics  of  the  streai  dined  body.  Similarly,  Reference  4-1 
has  provided  correlation  for  the  pitching  moment  characteristics  for  configurations 
similar  to  the  AG(EH).  Experience  indicates,  then,  that  the  blimt  based  pod  provides 
a  practical  approximation  for  the  streamlined  pod  with  its  wake  for  the  pitch  case. 


The  section  lift  curve  slope  presents  a  more  serious  unknown,  particularly  for 
the  16-series  section.  The  16-section  characteristics  employed  here  are  from  Refer¬ 
ence  4-2  but  all  of  the  data  of  Reference  4-2  except  that  at  a  Reynolds  number  of 

g 

.  85  x  10  is  subject  to  obvious  Mach  number  effects  while  any  laminar  flow  section 
data  is  subject  to  Reynolds  number  effects  at  Reynolds  numbers  of  less  than  about 
one  million.  Thus  the  use  of  section  data  measured  at  .85  million  RN  and  model  data 
measured  at  3  -  4  million  RN  to  predict  prototype  performance  at  30  -  40  million  RN 
reflects  significant  uncertainties  throughout  the  prototype  foil  characteristics. 


Reference  4-1  presents  the  pitch  lift  curve  slope  as  the  sum  of  foil,  body,  and 
nose  slopes; 


C 


L 


ot 

nose 


S  ’^W(B)  ^  ^  s  ’^B(W)  ^ 


(4. 1. 1) 


where: 

S  =  reference  area  =  area  of  complete  foil 
S’  =  area  of  exposed  foil 


SvT  =  pod  frontal  area  =  tt  D  /4 
N 

C  ’  =  exposed  foil  lift  curve  slope 

L 

Of 

D/b  =  pod  diameter/foil  span  ratio 

=  blunt  based  pod  lift  cur%^e  slope  based  on  frontal  area  =  2  from 

a  classic  potential  theory 
B 

=  pod  lift  carry-over  to  foil,  Chart  1,  Reference  4-1 

^B(W)  ”  carry-over  to  pod,  Chart  1,  Reference  4-1 

Eq.  (4. 1. 1)  is  evaluated  in  Table  4-1  for  the  AG  (EH)  forward  and  aft  configura¬ 
tions.  The  results  differ  from  that  for  the  complete  foil  neglecting  the  pod  by  3.37% 
fwd  and  1. 5%  aft  which  is  not  significant  but  the  values  of  Eq.  (4. 1. 1)  are  assumed 
hereafter  for  consistency  with  the  moment  calculations. 


Reference  4-1  presents  the  incidence  lift  curve  slope  as  the  sum  of  foil  and  body 


slopes; 


C  =  C 


S  'Sv(B)^L^  S  *'b(W)‘^L^  (4.1.2) 

where  and  are  from  Chart  1  of  Reference  4-1 

The  incidence  lift  curve  slopes  are  also  evaluated  in  Table  4-1.  The  incidence 
li;;t  curve  slopes  differ  from  the  pitch  lift  slopes  by  16%  fwd  and  22%  aft  which  is  a 
significant  difference  and  one  which  must  be  observed  as  required.  It  is  to  be  noted, 
hc‘wever,  that  the  AG  (EH)  aft  foil  and  pod  are  rigidly  attached,  the  incidence  control 
motion  is  between  pod  and  strut,  so  that  the  incidence  lift  curve  slope  is  never  required 
for  this  c^^nflguration. 

For  the  20%  chord,  sealed  gap  flaps  proposed  for  the  fwd  foils,  a  flap  effective¬ 
ness  of  .  47  was  taken  from  Figure  18  of  Reference  4-3  and  the  flap  curve  slope 
becomes: 


c  =  —  c 

L  d6  ‘"L. 


.47  X  2.49 


=  1.17  (4,1.3) 

The  lift  curve  slopes  of  Eqs.  (4. 1. 1)  -  (4. 1.3)  are  aerodynamic;  i.e. ,  infinite 
depth.  The  hydrodynamic  pitch  lift  curve  slopes  were  measured  for  the  fwd  and  aft 
foils  at  three  depths  and  two  speeds  each  for  Reference  4-4  with  the  results  presented 
in  Figures  1-30  and  I-'Sl  of  that  Reference  and  on  Figure  4-1  of  this  report.  Also 
presented  on  Figure  4-1  are  the  free  surface  effect  predictions  of  References  4-5  -  4-8. 
These  four  predictions  present  distinctive  aerodynamic  (infinite  depth)  lift  curve  slopes 
which  have  been  adjusted  to  the  lift  curve  slope  of  Eq.  (4. 1. 1)  on  Figure  4-1;  i.e. , 
each  prediction  has  been  employed  in  the  form  C  /C 

Li  Li 

a  a 

00 

In  considering  these  slopes,  it  must  be  recognized  that  hydrodynamic  lift  curve 
slopes  can  only  be  measured  with  a  precision  on  the  order  of  +  5%  so  the  entire  spread 
of  predictions  on  Figure  4-1  is  scarcely  significant.  Nevertheless  the  50-knot  forward 
foil  results  of  Figure  4-1  appear  doubtful  and  this  reflects  on  the  other  characteristics 
associated  with  these  measurements. 

The  incidence  lift  curve  slope  was  measured  only  for  the  fwd  foil  in  Reference 
4-4  but  there  are  independent  measurements  at  two  depths  in  Reference  4-9  and  these 
are  compared  with  the  predictions  on  Figure  4-2. 

The  tapered,  swept  foil  free  surface  effect  has  received  inadequate  attention  in 
the  literature  and  the  incidence  lift  case,  particulary,  presents  a  problem  in  selecting 
the  reference  chord  for  depth  effect.  The  MAC  is  preferred  over  the  MGC  as  a  refer¬ 
ence  in  this  report  and  on  Figure  4-2  the  MAC  for  the  entire  foil  is  employed  rather 
than  the  exposed  foil  MAC  --niply  because  it  provides  a  better  comparison  with  the 
predictions. 

From  Figures  4-1  and  4-2  the  Panchenl^ov  free  surface  effect  has  been  selected 
as  the  definition  for  this  effect.  By  emplojdng  a  common  definition  for  pitch  and 
incidence  lift  free  surface  effect,  the  ratio  between  the  two  lift  curve  slopes  becomes 
independent  of  depth  and  for  the  AG(EK)  fwd  foil  is  defined  by; 
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2.49 
2.  97 


.838 


(4.1.4) 


The  lift  curve  slope  Is  Independent  of  camber  so  the  experimental  results  con¬ 
sidered  thus  far  can  be  expected  to  be  representative  of  the  prototype  except  for  tlie 
unknown  Reynolds  number  effect.  In  the  consideration  of  foil  and  flap  angles  below, 
however,  we  must  recognize  that  the  model  and  prototype  have  different  cambers. 


4.  2  RESIDUAL  LIFT 


In  the  linear  lift  range  the  lift  curve  is  conveniently  expressed  in  slope-intercept 
form  where  the  Intercept  may  be  on  the  angle  of  attack  scale  at  the  zero  lift  angle  of 
attack  or  on  the  lift  coefficient  scale  at  the  residual  lift;  i.  e. ,  the  lift  when  pitch, 
Incidence,  and  flap  are  zero.  The  slope-residual  lift  form  of  the  lift  curve  is  most 
convenient  where  there  is  more  than  one  angle  involved  and  in  this  report  the  lift 
curve  is  considered  in  the  form 


C  =  C  (y+  C_  i  +  C  6  +  C 

^  h  ^i  ^6  ^0 


=  C, 


or 


-  J 


+  C, 


+  C, 


(4.2.1) 


eff. 


““pod 


C  is  not  known,  in  general,  in  the  presence  of  a  pod  and  particularly  in 


the  presence  of  a  free  surface.  It  is  approximated  in  this  report  by 


^  =  -C  a 

h  „  h 

eff 


(4.2.2) 


where  the  section  zero  lift  angle  of  attack  for  most  sections  is  available  in  Reference 

4-10.  The  only  comparable  source  available  to  this  report  for  the  16 -series  sections 

was  Reference  4-2  and  the  zero  lift  angles  in  this  report  are  from  an  unpublished 

0 

summary  of  the  data  of  Reference  4-2  for  a  .  3  Mach  number  and  .  85  x  10  Reynolds 
number. 
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Is  the  lift  due  to  pod  asymmetry  and  to  the  pod/foil  asymmetry  in  theory 
pod 

but  includes  dlfficiencies  in  the  approximation  of  Eq.  (4. 2. 2)  in  practice. 

For  the  fwd.  foil  model  of  Reference  4-4,  Eq.  (4.2.1)  becomes 


Ct  =  C_  \a  +  - — i  “C  01  A  +  C_ 


a  +  - —  1 

a 


0  c. 


(4.2.3) 


and  at  zero  lift 


'L=  0 


(4, 2. 4) 


The  measured  zero  lift  incidence  angles  of  Reference  4-4  are  presented  in  the 

form  of  Eq.  (4. 2.4)  on  Figure  4-3.  The  difference  between  the  35  and  50  \not  results 

on  Figure  4-3  is  incredible,  particularly  for  deep  submergence,  and  the  50-knot 

results  are  discarded  because  they  are  suspect  on  Figure  4-1.  The  35-knot  results 

indicate  that  the  pod  lift  is  a  function  of  depth  in  the  same  way  that  C  is;  i.e. , 

^i 

C  /C  is  a  constant  and  the  pod  lift  can  be  expressed  as 
^pod  ^i 


where 


'L  ,  ■  C  =  0) 

P"''  ‘  "^pod 


a  c  =  0  =  -  2.25  -  (-2.92) 
pod 

=  .67° 


(4.2.5) 


and  where  this  angle  is  assumed  to  be  independent  of  fhe  foil  camber. 


For  the  aft  foil  model  of  Reference  4-4  Eq.  (4. 2. 1)  becomes 
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c 
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O' 


c .  = 


and  at  zero  lift 


O' 


(4.2.6) 


(4.2.7) 


The  measured  zero  lift  angles  of  Reference  4-4  for  the  aft  foil  model  are  pre¬ 
sented  on  Figure  4-4  in  the  form  of  Eq.  (4. 2.  7).  The  results  are  rather  badly  scat¬ 
tered  and  for  this  report  it  is  assumed  from  Figure  4-4  that  the  pod  lift  is  zero. 

The  foil  lift  curve  results  are  summarized  in  Table  4-1. 

4.3  FOIL  PITCHING  MOMENT  SLOPE 

Grumman  Aerospace  Corporation  practice  with  pods  is  to  locate  the  shoulder, 
the  first  point  of  maximum  diameter,  at  the  foil  root  chord  leading  edge  (center  of  pod). 
Whether  that  practiv-^e  is  reflected  in  the  ”as-built”  AG  (EH)  or  not  is  uncertain.  Un¬ 
defined  pod  changes  are  noted  in  the  design  history  (Reference  4-11)  and  Figures  1 
and  1. 1  of  the  model  test  report  (Reference  4-4)  conflict  on  this  point.  Prototype 
production  prints  available  to  date  do  not  clarify  the  issue.  It  must  be  noted  that  the 
methods  of  Reference  4-1  are  invalid  if  the  shoulder  is  not  forward  of  the  exposed  foil 
root  chord  leading  edge. 

In  this  report  the  terms  ’’forebody”  and  "nose"  are  considered  synonymous  and 
the  pod  is  divided  into  a  forebody,  prismatic  section,  and  afterbody.  Thus  Reference 
4-1  terminology  is  related  to  the  terminology  of  this  note  as  follows: 


It  must  be  noted  that  the  "centers  of  pressure"  of  Reference  4-1  are  centers  of 
pressure  only  tor  foils  of  symmetric  section  or  in  the  context  of  the  thin  airfoil  theory 
of  that  reference.  In  the  more  general  case  of  the  cambered  foil,  the  centers  of  pres¬ 
sure  of  Reference  4-1  must  be  designated  as  aerodynamic  centers  and  they  are  so 
regarded  in  this  report. 

Reference  4-1  defines  the  moments  acting  on  the  total  configuration  by  locating 
the  aerodynamic  center  associated  vdih  each  of  the  three  loads  identified  by  Eqs. 
(4.1.1)  and  (4.1.2). 
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Reference  4-1  states  that  the  aerodynamic  center  for  the  foil  lift  is  unknown  in 
the  general  case  and  is  best  approximated  by  employing  the  aerod3mamic  center  for 
the  exposed  foil.  Referring  to  Figures  1-2  and  1-3,  then,  the  AG(EH)  foil  lift  aero¬ 
dynamic  centers  are  located  on  the  total  foil  mean  aerodynamic  chord  at 


a.c. 


foil 


8.10  -  4.93  3.17 

- WS -  =  9733  = 


a.  c. 


4. 167  -  2. 245  1.  622 


foil 


4.81 


4.81 


=  .337 


In  general  the  nose  lift  aerodynamic  center  is  located  by 


I  =  £  -  V  /S 

N  S  N 


(4.3.1) 


(4.3.2) 


Reference  4-1  gives  the  aerodynamic  center  for  ogival  noses  in  Chart  9,  Ellip  ¬ 
soidal  noses  were  specified  for  the  AG  (EH)  and  are  assumed  here  though  the  ”as  built" 
nose  lines  were  not  available  for  this  report.  For  ellipsoidal  noses,  the  aerod3niamlc 
center  is  located  1/3  of  the  nose  length  aft  of  the  leading  edge.  Referring  to  Figures 
1-2  and  1-3,  the  AG  (EH)  nose  aerod3mamic  centers  are  located  on  the  total  foil  mean 
aerodynamic  chord  at 


a.c. 


nose. 


=  -C 


4. 93  +  2/3  X  3. 763 
9.33 


4.93  +  2.51 
9.33 


7.44 

9.33 


=  -  .797 


a.c.  nose^  = 


Y2.545  +  2/3  x  2.867\  2.545  +  1.91  _  4.555  _ 

"V  4.81  /  ”  ”  d  Ri  ~  " 


4.81 


4.81 


(4.3.3) 

The  location  for  the  body  lift  aerodynamic  center  is  a  tedious  calculation,  par¬ 
ticularly  for  low  aspect  ratio  foils.  These  calculations  are  presented  in  Appendix  A 
where  Eq.  (A  16)  gives  the  result: 


a.  c. 


pod^ 


a.  c. 


.1325 

.  163 


(4.3.4) 
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In  general  these  aerodynamic  centers  in  combination  with  the  lift  curve  slopes 
of  Table  4-1  produce  a  total  aerodynamic  center  defined  by 


a.c. 


H 

c 


(C  orC  ) 
or  C  ) 


(4.3.5) 


where  the  terms  of  the  numerator  and  denominator  are  determined  by  the  type  ot  lift 
Involved  and  by  the  application  to  be  made  of  the  aerodynamic  center  and  where  the 
terms  of  numerator  and  denominator  do  not  necessarily  correspond. 

Noting  that 

C„  =  (-  -  a.c.)  C,  orC  =  (-  -  a.c.)  C 
H  \c  'L  H.  \c  'L, 

a  a  i  i 


dC 


H  H 


=  —  -  a.c. 


Eq.  (4. 3.  5)  simply  states  that 


total  C 


^  component 

_ a 

H  ^  component  C 


O' 


(4.3.6) 


(4.3.  Y) 


The  forward  foil  hinge  moment  slopes  of  interest  to  the  AG  (EH)  are  numerically 
evaluated  in  Table  4-2  which  illustrates  the  application  of  Eq.  (4.3. 7)  for  the  case 
where  incidence  lift  control  is  achieved  by  moving  the  foil  relative  to  the  pod.  For  the 
AG  (EH)  aft  foil,  where  foil  and  pod  are  rigidly  attached,  there  is  no  ’’incidence  lift” 
case  and  Eq.  (4.3.7)  becomes  simply: 


-  a.c.  )  C  +(--a.c.  )c  +  (-  -  a.c.  )c 

\c  foil/  L  \c  pod/  L  \c  nose/ 

o .  a  _ 


foil 


pod 


nose 
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(.336  -  .337)  2.385  -h  (.336  ■  .163).  501  (.336  .  925)  .  1326 

3.015 


-  .001  X  2.385  +  .173  X  .501  +  1.261  x  .1326 

3.015 

-  .002385  ^  .0868  +  .1671  _  .2515 

3.015  "  3.015 


=  .0834  (4.3.8) 

The  hinge  moment  slopes  of  Figures  1-8  through  1-17  of  Reference  4-4  are  pre¬ 
sented  on  Figure  4-5  as  aerodynamic  center  for  this  case  from  Table  4-2.  Figure  13 
of  Reference  4-15  indicates  that  there  is  an  uncertainty  range  of  about  +  5%,  -  1%  of 
MAC  associated  with  the  aerodynamic  center  prediction  and  this  range  has  bracketed 
previous  experience  with  this  characteristic.  Figure  4-5  indicates,  however,  that  the 
prediction  error  in  this  case  barely  lies  within  even  the  more  conservative  uncertainty 
range  of  Reference  4-1. 

The  hinge  moment  aerodynamic  center  prediction  for  pitch  lift  of  Table  4-2  is 
compared  with  the  measured  slopes  of  Figures  B1  through  B6  of  Appendix  B  on  Figure 
4-6.  The  result  is  the  same  as  for  incidence  lift;  the  two  cases  present  about  the 
same  a.c. ,  as  theory  indicates,  but  both  aerodynamic  centers  are  well  forward  of  the 
expected  position.  Similarly,  the  measured  aft  foil  aerodynamic  center  measured  by 
the  slopes  of  Figures  1-4  through  1-7  of  Reference  4-4  lie  well  ahead  of  the  predicted 
position  as  shown  on  Figure  4-7. 

The  suspect  term  in  these  calculations  is  the  exposed  foil  aerodynamic  center 
and  the  ’’measured”  characteristics  of  Table  4-2  proceed  from  the  measured  charac¬ 
teristics  by  appropriate  adjustment  to  the  exposed  foil  aerodynamic  center.  It  is 
proposed  that  this  aerodynamic  center  location  be  estimated  by  integrating  the  pitch 
and  incidence  spanwise  load  distributions  outboard  the  pod.  This  is  not  a  difficult 
procedure  but  time  does  not  permit  such  a  procedure  for  this  report. 

4. 4  RESIDUAL  FOIL  PITCfflNG  MOMENT 

The  moment  about  the  aerodynamic  center  is  assumed  to  be 
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(4.4.1) 


a.  c. 


MAC*  S» 
MAC  S 


eff 


a.c. 


~  MAC*  S* 
MAC  S 


eff 


where,  from  an  ^mpublished  compilation  of  the  .  3  Mach  Number,  .  85  x  10  Re5molds 
number  data  of  Reference  4-2,  the  effective  design  lift  coefficient  for  an  8%  16-series 
is  taken  to  be  83%  of  the  design  lift  coefficient. 

Then  the  hinge  moment  for  the  AG  (EH)  aft  foil  may  be  written: 


o  0 


(4.4.2) 


where 


C„  =  C..  +  C_. 

H.  M  M  , 

0  a.c.  pod 


C^  is  the  unpredictable  pod  effect, 
pod 


To  account  for  depth  effect  on  the  lift  curve  slope,  this  equation  roust  be  employed 
in  the  form: 


C„  =  -  a.c.)  +  C„ 

H  L  \  c  /  H . 

Qi  C 


(4.4.3) 


At  zero  lift  the  aft  foil  hinge  moment  Is: 


”c  =  0  ^ 

Lj 


/H  N  /  S\ 

(e  -  !-  j-r  s, 

“  \  V  ' 
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^TT 

’  ^0  ^ 


=  C 
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+  C 
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eff 


pod 


(4.4.4) 
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This  equation  is  compared  with  the  measured  model  zero  lift  hinge  moments  on 
the  lowest  graph  of  Figure  4-8.  The  reduction  in  hinge  moment  indicated  on  that  plot 
for  a  fixed  transition  is  in  the  expected  direction  but  the  possibility  of  premature  transi¬ 
tion  on  the  prototype  has  profound  implications  for  drag  and  cavitation,  as  well  as  for 
incidence  and  flap  moments,  which  are  beyond  the  scope  of  the  present  study. 

The  systematic  variation  of  zero  lift  hinge  moment  with  depth  is  of  interest 
here.  In  an  attempt  to  account  for  this  variation,  £q.  (4.4.4)  was  considered  in  the 
form; 


=  C 
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+  C 


^L=  0 


a^  c« 


M 


eff. 


pod 


(4.4.5) 


The  result,  also  presented  on  Figure  4-8,  indicates  that  the  moment  coefficient 
about  the  aerodynamic  center  is  subject  to  the  same  free  surface  effect  as  the  lift 


curve  slope;  i.e. ,  /C 


"M 


a.  c«  a.c. 


=  C  /C  .  From  Figure  4-8  a  value  of 
L  ij_ 


O'  O' 


-.  0G55  is  taken  as  C 


M 


pod 

For  the  more  general  case  with  incidence  motion  relative  to  the  pod,  the  hinge 
moment  is 

a  i  0 


where; 


C„  O'  +  C„  1  +  C-,  + 

H  H,  M  M 

o  i  a.c.^^^  pod  0 


"M 
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(4.4.6) 
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=  1.  0  chord  line 


4-12 


C  =  ,  83  C  for  8%  16-series  section  on  a  =  1, 0 
eff 

=  unpredictable  pod  effect 
pod 

For  the  purpose  of  examining  the  model  data  Eq.  (4. 4. 6)  may  be  written: 

C  =  C_  -  a.c.  )  o  +  C  -  a.c.  )  i  +  C„ 

H  L  \c  o  /  L  \  c  i„/ 


Cl  =  0 


=  C  (—  -  a.c.  )a  +  C  (—  -  a.c.  )  i 
L  \  c  a,}  L,  \  c 


(4.4.7) 


where: 


C  =  C,  a  +  'j^  1  +  =0 

L(  L  L,  L- 
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%‘c,  =0=-\-‘^L  “ 
i  L  0  a 


,  0  O' 
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(4.4.S) 


Then  Eq.  (4. 4. 7)  may  be  written: 
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0  a.c. 

eff 


+  C  +  C  C 
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(4.4.9) 
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For  the  purpose  of  evaluating  the  residual  hinge  moment  from  the  towing  tank 
data,  Eq.  (4. 4.  9)  was  considered  in  the  form  (which  assumes  that  C  is  subject  to 

d  o 

00 


■'M 


*  ■  CT^  „ 
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C  -  (-.c.  -  a.c.  )  C  a 

V  =  0  'h  ^d 

d  L  “ 


(4.4.10) 


The  towing  tank  data  is  presented  in  the  form  of  Eq.  (4.4. 10)  on  Figure  4-9 
which  contradicts  the  results  of  Figure  4-8  with  respect  to  the  effects  of  the  pod  and 
of  the  fixed  transition.  From  Figure  4-9  a  value  of  +  .  0151  is  taken  as 

pod 

4.5  COMPARISON  WITH  PROTOTYPE  DATA 

Reference  4-14  provided  independent  measures  for  the  AG  (EH)  forward  and  aft 
lift  and  hinge  moment  curves  which  can  be  compared  with  theory  and  with  the  towing 
tank  results.  On  the  craft,  however,  the  forward  and  aft  lift  and  hinge  moment  are  all 
interdependent  and  the  buoyancy  for  the  submerged  configuration  becomes  of  some 
interest. 


Foil  buoyancy  is  given  by 


®foll 
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where  K  =  section  area/ct 


(4.5. 1) 
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The  buoyancy  coefficient  is  a  function  of  the  configuration  size,  b,  and  is 
negligible  in  model  scale  but  of  some  significance  on  the  prototype. 

For  the  AG  (EH)  the  numerical  evaluation  for  Eq.  (4.  5. 1)  is; 

8  .735  _  b  1.39 


=  1.384 


V‘ 


C„  =  1.384  X  26/V^  =  35.95/V^ 

C„  =  1.384  X  13. 416/V^  =  18.58/V^ 


(4.5.2) 


,  Much  of  the  interior  of  the  prisimatic  portion  of  the  pod,  excluding  the  foil 
structure,  is  wet.  Consequently,  the  buoyancy  for  the  prismatic  portion  of  the  pod 
is  approximated  here  by  employing  the  buoyancy  of  the  total  foil  and  neglectir.g  the 
prismatic  pod  portion  entirely. 

The  pod  buoyancy  employed  here,  then  is 
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(4.5.3) 
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By  reference  to  Sheet  22  of  Reference  4-16  and  to  Foil  Contract  Guidance  Plan 

(AG(EH)  800-2068282,  Hydrofoil  System  Geometry  -  Structural  Arrangement,  It  was 
determined  that: 


d^  =  16.32  -  height  (any  pitch) 

d^  =  14.82  -  height  (pitch  =  0)  (4.5.4) 

The  pitch  reference  axis  Is  assumed.  For  the  condlUons  of  Figures  1  -  3  of 
Reference  4-14,  then,  the  foil  depths  are: 

d^  =  11,32  =  1.215  MAC 


dg  -  9.82  +  125,87/57.3  =  12.015  =  2. 495  MAC 
The  average  submerged  strut  thickness  ratio  Is 
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(I).  ■  ■ 


1  +  ,00238 


(--f) 


Then  the  strut  buoyancy  is  given  by 
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(4. 5.  6) 
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^^"^(12.015  -  1.125)[.l  +  .00238(12.015  -  1.125)]  =38.9/V^ 


(4.5.7) 
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and  the  total  buoyancy  coefficients  become 

C„  =  (35.95  +  22.30  +  34. 15)/V^  =  92.6/V^  =  32.5/V^^ 

C„  =  (18.58  +  17.45  +  38.9)/V^  =  74.9/V^  =  26.3/V„^  (4.5.8) 

®2  ^ 

With  the  buoyancy  coefficient  defined,  the  balance  of  lift  forces  on  the  craft  may 
be  written: 
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The  moment  balance,  neglecting  the  thrust/drag  moment,  may  be  written: 
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For  many  purposes,  Eqs.  (4, 5. 9)  and  (4. 5. 10)  may  be  employed  in  a  form  wMch 
neglects  buoyancy  and  hinge  moment; 
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The  parameters  required  for  Eqs.  f4.  5.  9)  -  (4. 5. 11)  are  evaluated  numerically 
in  Table  4-4  for  the  weight  and  balance  of  Figures  1  -  3  of  Reference  4-14  and  four 
forms  of  the  equations  are  presented  in  numerical  form  for  the  same  weight  and 
balance  in  Table  4-5. 


The  variation  of  f(  rward  and  aft  incidence  angle  with  speed  was  calculated  for 
each  of  the  four  sets  of  equations  of  Table  4-5  with  the  result  presented  on  Figure 
4-]  1.  The  results  demonstrate  that  foil  moment  and  buoyancy  are  both  significant  to 
crau  of  this  size  compared  with  craft  of  the  PGH  size,  for  example,  where  neither 
effect  had  any  practical  significance.  The  measured  foil  angles  of  Reference  4-14 
are  also  shown  on  Figure  4-10  where  agreement  with  theory,  and  the  towing  tank,  is 
reasonable  forward  and  negative  aft. 

For  the  theoretical  incidence  angles  of  Figure  4-10,  with  foil  moment  and 
buoyancy,  the  corresponding  hinge  moment  variation  with  speed  was  calculated  by 
the  methods  of  Tables  4-2  and  4-3  with  the  result  presented  on  Figures  4-11  and  4-12. 
Two  theoretical  hinge  moments  appear  on  these  figures  because  the  towing  tank  aero- 
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dynamic  centers  do  not  agree  with  theory.  The  prototype  measured  hinge  moments 
are  also  shown  on  these  figures  and  the  agreement  with  theory  is  not  satisfactory  for 
the  forward  or  aft  foil. 


A  negative  increment  in  hinge  moment  due  to  buoyancy  could  be  added  to  the 
theoretical  ounces  of  Figures  4-11  and  4-12.  There  is  a  mean  buoyant  chord  just  as 
there  is  a  mean  aerodynamic  chord.  The  mean  buoyant  chord  may  be  defined  as  the 
chord  of  a  rectangular  foil  having  the  same  area  and  buoyancy  as  the  total  foil,  i.e. , 
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(4.5.13) 


It  will  be  noted  that  the  MBC  is  located  at  mid-semi-span  on  a  rectangular  foil, 
evaluated  by  taking  first  derivative  of  numerator  and  denominator  with  respect  to  X , 
and  at  1/3  the  semi-span  for  a  zero  taper  ratio  as  would  be  expected. 


The  effect  of  the  pod  nose  and  afterbody  buoyancy  on  the  aft  hinge  moment  are 
neglected  here  but  would  have  to  be  examined  on  a  larger  craft  of  this  configuration. 
Similarly  the  section  center  of  gravity  would  have  to  be  established  on  a  larger  craft 
but  it  is  assumed  to  lie  at  mid-chord  here. 
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As  the  foil  buoyancy  is  approximated  here,  therefore,  the  forward  and  aft 
centers  of  buoyancy  are  identically  located  at  48. 6%  of  the  MAC  and  contribute  hinge 
moments  of 

M„  =  qSMAC 

“b  ”b 

=  c„  -  b.c.)  q  SMAC 

foil 

OE  qc  2  6** 

=  ^^^(.385  -  .486)  X  .995V  x  225  x  9.33  =  -  7,580  =  -  .091  x  10  #  fwd 

=  i5i|i(.336  -  .486)  x  .  995V^  x  60  x  4. 81  =  -  800  =  -  .0096  x  10^”#  aft 

(4.5.14) 

Neither  of  these  results  is  significant  to  the  results  of  Figures  4-11  and  4-12. 

The  towing  tank  and  prototype  lift  and  hinge  moment  results  are  summarized  on 
Figures  4-13  and  4-14  and  on  these  figures  the  prototype  weight,  less  buoyancy,  has 
simply  been  distributed  on  the  hinge  location  to  consen^e  time.  The  forward  foil 
results  of  Figure  4-13  indicate  reasonable  agreement  between  theory,  towing  tank, 
and  prototype  for  lift  but  poor  for  hinge  moment.  The  towing  tank  hinge  moments  are 
low  and  outside  the  range  of  previous  experience  with  theory.  The  aft  foil  prototype 
lift  results  of  Figure  4-14  are  high  and  of  suspect  slope  and  zero  lift  angle.  The  pro¬ 
totype  hinge  moments  are  higher  and  the  towing  tank  moments  lower  than  theory  with 
both  sets  of  results  being  outside  the  prediction  error  results  for  previous  results 
with  the  theory. 

4.6  FLAP  HINGE  MOMENT 

In  the  nomenclature  of  this  Volume,  Theodorsen’s  Eq.  XIX  of  Reference  4-17 
becomes,  for  steady  state  incidence  and  flap  deflections: 
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£q.  (4. 6. 1)  provides  2-dimensional,  thin  airfoil  flap  hinge  moment  slopes  which 

are  compared  on  Figures  4-15  and  4-16  with  two  other  predictions  for  these  slopes 

which  provide  accountability  for  the  section  thickness.  The  DATCOM  prediction  of 

these  figures  is  from  Reference  4-12  and  Toll’s  prediction  is  from  Reference  4-3. 

The  16-series  /c^  employed  was  from  an  unpublished  compilation  of 

a  a  theoretical 

the  lowest  Mach  Number  data  of  Reference  4-2;  that  for  the  65-series  section  was 
derived  from  Section  4. 1. 1  of  the  DATCOM.  Figures  4-15  and  4-16  provide  identi¬ 
fication  for  Theodorsen’s  ”T”  parameters  in  terms  of  practical  coefficients. 


Three  procedures  for  converting  the  2-dimensional  flap  moment  slopes  to 
3-dimensional  are  considered  in  Appendix  C.  The  state-of-the-art  provides  little 
guidance  in  this  particular  area  and  Perkin's  and  Hage's  procedure  was  adopted 
simply  because  it  is  the  simplest  procedure  and  the  most  satisfying  intuitively. 


For  C,  ,  Perkins  and  Hage  simply  employ 
n 
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(4.6.2) 


o 

and  in  considerati  ?f  the  uncertainty  associated  with  this  derivative  there  is  no  real 
loss  of  accuracy  if  the  theoretical  section  slope  is  employed 
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For  C,  Perkins  and  Hage  add  the  difference  between  the  2-  and  3 -dimensional 

C,  '8,  reduced  by  the  flap  effc'^Hveness,  to  the  2-dimensional  flap  effertiveness  and 
h 
a 

that  procedure  is  adopted  here  without  further  comment: 
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For  the  residual  flap  hinge  moment  this  volume  simply  employs; 

C.  =  -  C  /2 
•*0  ^i 

eff 


(4.6.5) 


noting  that  this  value  is  very  probably  an  upper  (negative)  bound  for  the  true  value 
and  that  a  better  approximation  might  be  provided  by: 
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(for  triangular  load  distribution)  (4. 6. 6) 


It  must  be  recognized  that  the  flap  hinge  momen*,,  being  an  intimate  function  of 
the  boundary  layer  characteristics  at  the  foil  trailing  edge,  is  the  most  uncertain  of 
all  the  foil  characteristics  and  this  remains  true  with  or  without  model  test  guidance. 

The  AG  (EH)  flap  hinge  moment  is  numerically  evaluated  in  Table  4-6.  These 
characteristics  were  measured  in  the  towing  tank  but  the  results  are  not  available 
at  this  writing. 

4.  7  FOIL  PITCHING  MOMENT  DUE  TO  FLAP 

In  the  nomenclature  of  the  volume,  Theodorsen’s  Eq.  XX  of  Reference  4-17 
becomes  for  a  steady-state  flap  deflection: 


p  2  2  ^4  ^10 

fe  V  X  - ^ - 


6  +  2n  ^  c  V 


£  ,2  ^2  JH  _  1',  _10 


2  “  ’  \c  4/  n 
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(4.7.1) 


The  2n  coefficient  of  the  first  term  is  the  aerodynamic,  2- dimensional  lift 
curve  slope  and  Eq.  (4. 7. 1)  may  be  written  more  generally  as: 


C 


H 


-  -nr  + 

2'  4 


(4.7.2) 


a 

Eq.  (4. 7. 2)  presents  flap  lift  as  the  sum  of  two  chordwise  load  distributions, 
one  identical  with  the  load  distribution  associated  with  the  section  angle  of  attack  and 
having  its  center  of  pressure  at  the  section  aerodynamic  center  and  the  other  being 
concentrated  on  the  flap  hinge  line.  The  principle  is  employed  in  Section  6  of  this 
volume  to  calculate  the  pressure  distribution  over  the  flapped  foil  and  it  is  discussed 
in  good  detail  in  Reference  4-18. 


The  first  term  of  Eq.  (4. 7. 2),  then,  is  the  product  of  the  lift  coefficient  due  to 
flap  deflection  and  the  distance  between  the  foil  incidence  hinge  and  the  aerodynamic 
center  for  one  component  of  the  flap  lift.  Where  the  flap  extends  the  full  span  of  the 
exposed  foil  that  aerodynamic  center  is  identical  with  the  aerodynamic  center  for 
incidence  lift. 


Theodorsen’s  is  obviously  the  theoretical  flap  effectiveness,  da/d6, 

graphically  presented  for  example  on  Figure  18  of  Reference  4-3  which  also  presents 
practical  values  for  the  flap  effectiveness. 


As  Theodorsen  defines  T^  and  T^^  in  References  4-17,  the  second  term  of 


Eq.  (4. 7. 2)  may  be  written; 


-  i<T4  "  «  =  -  2 


h  C 


M 


a.  c. 


d  6 


(4.7.3) 


which  is  the  approximation  employed  in  Eq.  6. 1. 2. 1-a  of  Reference  4-12,  for  example, 

for  the  moment  slope.  The  approximation  is  2-dimensional  and  the  3-dimenslonal  slope 

is  unknown  but  reduction  of  the  2-dimensional  slope  by  the  ratio  C  /C  should 

L.  L 

^  ^2-D 

provide  a  reasonable  approximation  for  the  3-dimensional  effect. 
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With  these  considerations  incorporated,  a  practical  representation  for  Eq. 
(4.7.2)  becomes 


For  the  proposed  AG  (EH)  flap  system,  the  hinge  moment  slope  for  flap  deflec¬ 
tion  becomes 

-  |(-  .1424  +  .5498) 

-  .1018  (4.7.5) 


The  incidence  hinge  moment  due  to  flap  deflection  for  the  proposed  AG(EB 
flap  system  has  been  measured  in  the  towing  tank  but  the  results  are  not  available 
at  this  writing. 
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TABLE  4-1  AG  (EH)  FOIL  LIFT  CHARACTERISTICS 


TABLE  4-2  AG(EH)  FOIL  MOMENT  SLOPES 


PREDICTED 


H 

—  -  a.c. 
c 


COMPONENT  H/c  a.c.  =  C, 


PITCH  LIFT  -  FWD  FOIL 


Foil  I  .340  I  .045 

Pod  .385 

Nose  I-  .797  11.182 


Foil  +  pod  +  nose  =  C„ 

n 

®trim 

Used  to  trim  craft  moments 


FoiJ  only 

Control  force  hinge  moment  slope 


COMPONENT  H/c 


2.52 

.1135 

.374 

.0945 

.0776 

.0918 

.2998 

^trim 

Cl 

2.97 

a 

trim 

.101 

a.c. 

°^:rim 

.284 

a 

.1135 

Cl 

2.97 

a 

\ 

.0382 

a.c. 

"h 

.3468 

O 

_ 1 

i 

Ot  C, 


.2375  .0943  .291 


Foil 


Pod 


Foil  +  pod  =  C, 


INCIDENCE  LIFT  -  FWD  FOIL 


5  2. 175  I  .  0979 


.1325  .2525  .306  .0773 


trim 

Used  to  trim  craft  moments 


.1752 


2.49 


.1745  .0801  .305 


.0773 


.2518 


2.49 


.1012 


.2838 


TABLE  4-2  (Continued) 


-  a.c. 


COMPONENT  H/c  a.c.  =  C, 


Foil  only  =  C„ 

Control  force  hinge  moment  slope 


(Note  that  a.  c. ,  =  a.  c. .  ,,  for 

‘h 

incidence  lift) 


PREDICT! 

:d 

.0979 

% 

2.49 

.0393 

.3457 

MEASURED 


-  =  ,336 
c 


.1745 


2.49 


,070 


O' 

.2515 

.588 

or 

3.015 

3.015 

.0834 

.195 

a.c. 

.2526 

.141 

TABLE  4-3  AG(EH)  FOIL  MOMENTS 


TABLE  4-4  NUMERICAL  EVALUATION  FOR  TRIM  PARAMETERS 


=  13.8  MAC^  =  9.33 


THEORY/TOW  TANK 
.  lOl/.  1427 


trim. 


trim. 


.0705/.  1012 


S  =  510  t  =  125.9 


MAC, 


=  .0741 


.  0834/.  195 


MAC„ 


=  .03825 


d  =  li.32  C,  ^  2.8 

1  L 


.  1044 


.  32.5/v^ 


d  rr  12.02  C,  =  2  35 

^  L. 


.943 


C,  =  .1165 

\ 


MAC, 


=  .979  Cjj  =  -  .0573/-  .0431 


0. 


2 

MAC„ 


2.495 


C  =  -  .0490-  .0543 
°2 


57.3\V  ^  57.3  X  2240  x  295  ^  gO,  100/V  ^ 
q  S  >  *  1/ 


Craft  lift  curve  slope,  C  =  ^  C 


L  S  i.  S  L 
1  "2 


--  .882  X  2.8  ♦  .1176  x  ?.  95  =  2.47  ♦  .347  =  2.82 


Static  stability,  C.,  =  *  2.47  -  .891  x  .347  =  .2705  -  .3095  =  -  .039 

M  •'Li  *  S  L 

"l  '^2 

^  "2 

Craft  residual  lift,  C  -  .882  x  .  1044  •  .  1176  x  .  1165  =  .09208  ♦  .01.37  -  .1058 

°  ‘"2 

f  S  !  B 

11  2  2 

Craft  residual  lift  moraent,  "jT  *  ~~^L  '  1095  x  .09208  -  .891  x  .0137  =  .01008  -  .01221 

Craft  residual  momeni,  Cj^j  ‘  ‘"T^S^^Il  ''  ‘  '  ' O64.3) 


0„ 


M  1  2 

-  .  00375  -  .0002205  _  -  .00397  Theory 

-  .  00282  -  .  000244  ’  -  .  003064  Towing  Tank 


Craft  buoyancy  coefnclcnt,  C„  =  7^  C,,  •  ;pC,,  -  (.882  x  .32.5  *  .  2276  x  29. 1)/V..^  =  (2,87  *  3.42)/V'  ^ 

o  ^  15  n  15..  K  K  ,,  2 

'l  '"1  '2^2 

Craft  buoyant  moment,  C  —  Tr^'n  "  '  <*1095  x  28.  7  -  .  891  x  3.42)/V’  2 

Mjj  S  »  b  K 

(3.  143  -  3.047)/V'  ^  -  .096/V^^ 

K  K 


.00213 
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TABLE  4-5  NUMERICAL  EVALUATION  FOR  TRIM  EQUATIONS 


« 


NO  FOIL  MOMENT  OR  BUOYANCY 

LIFT:  1  “  :  I  °  --  .  C  s'’  -  57.3  C 

-  S  L,^  I  S  2  q  S  L„ 

.882  X  2.35  1  °  ♦  .1176  x  2.9.5  1°  =  26,100/V  ^  -  2.82  x  1  -  57.3  x  .  1058 
1  2  K 

2.07  -  .347  i°  =  26,  lOO/V.,'  -  2.82  -  6.07 

12  K 


J 


=  26,  100/V.,  -  8.89 

K 

/  S  IS 

r  r^L.  -  T  r^L  ‘2°  =  -  Si 

‘1  “2  “  \ 

.1095  X  2.07  I  -  .891  x  -347  12°  =  .  039  x  I  +  57.3  x  .00213 

.227  1^°  -  .3095  12°  =  .039  +  .122 
=  .161 

PREDICTED  FOIL  MOMENT  -  NO  BUOYANCY 

HII=  T\  '1°  ^  rS  '2°  =  ^  -  S 

2.07  1  °  >  .347  1  °  =  26.100/V.,^  -  8.89 
X  2  K 


1 


MOMENT; 


►  S  L«  I  J  1  S  b  I  J 


trim. 


^  -S  '2 

S,  °2 


[MAC,  S, 

- ^  -^C 

£  S 


"M 


MAC2  f2^ 


]' 


trim. 


£  S  ''H 

^l.  °2 

*■'2 


■  ‘  c,  , 

L  M 

.882  [.1095  +  .0741  x  .O705J  x  2.35  1^°  -  .  1176  1.891  -  .03825  x  .0834]  x  2.95  12° 

=  .039  -  [.0741  X  .882  X  .  101  x  2.8  +  .''*^625  x  .1176  x  .0834  x  2  95)  x  1 

-  57. 3(-  .00213  -  .00397) 

.882(.  1095  +  .0052)  x  2.35  1^^°  -  .1176(.891  -  .003)  x  2.95  12°  =  .039  -(.01852  +  .00110)  +  .349 
.832  X  .1147  x  2.35  1^°  -  .1176  x  .888  x  2.  95  12°  =  .  039  -  .01962  +  .349 

.2375  1  °  -  .308  1  °  =  .3684 
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TABLE  4-5  (Continued) 


PREDICTED  FOIL  MOMENT  WITH  BUOYANCY 


Hii;  T\  'i“  *  r^L  '2“  =  ^ 


2.07i  °  ^  .347  =  2G.100/V'  “  -  8.89  -  57.3  x  32. 12/V 


K 


=  (26,100  -  1837)/\'  -  8.89 

K 

=  24. 260/ V  ^  -  8.89 
K 

MOMENT;  SUBTRACT  57.3  FROM  RIGHT  SIDE  OF  PREVIOUS  EQUATION 

.  2375  1  °  -  .308  1  °  =  .3684  -  57.3  x  .  096/V^^ 

12  K 

=  .3684  -  5.49/V'  ^ 

K 

TOWING  TANK  FOIL  MOMENT  WITH  BUOYANCY 


LIFT:  UNCHANGED 


2.07  i°  ^  .347  1.  =  24,260/\'  ^  -  8.89 
12  K 


£  MAC 

MOMENT;  j-  ^  C^^ 


=  .  1095  +  .0741  X  .1012  =  .  1095  -  .0075  = 


trim , 


£„  MAC„ 
2  2 


£  H 


C„  =  .891  -  .03825  x  .195  =  .891  -  .00746  =  .884 


n 


MAC^  S^ 

~  s"  '  H 


MAC  S„ 

c  + - -  c  c 

_  ^L  £  S  ^L 

"1  ^L,  "2 


=  .  0741  X  .882  X  .1427  x  2. 8  +  .  0::825  x  .  1176  x  .  195  x  2.  95 
=  .0261  +  .00259  =  .02869 

.882  X  .117  X  2.35  -  .  1176  x  .884  x  2.  95  12° 

=  .039  X  1  -  .02869  x  1  -  57, 3(-  .00213  -  .003064)  -  57.3  x  .oUG/V, 

.2425  1  °  -  .306  1  °  =  .039  -  .02869  +  57.3  x  .00519  -  5.49/V  ^ 
12  K 


K 


-  .01031  +  .297  -  5.49/V, 


K 


=  .3073  -  5.49/V 


K 


.1170 
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TABLE  4-6  FLAP  HINGE  MOMENT  -  FWD.  ONLY 


=  C.  o+C.  S  +  C 

"h  “I"* 

Of  6  0 


=  "  k  -f,K  -  £h  ‘  -  "^L. 


Of 


ITS 


'h.  d6  '“h  ^ 

6  a  Of 


, •  k  -  n  k  - 

<y  L  6  '  a 


eff. 


— i  C 
2tt  h 


a'J 


6  -  C-  /2 

1  « 
eff. 


2tt  h 


C,  a  +  C 


dof 


i 

h - ij 


Of 


h,  d6  h  X  2tt 

6  Of 


6  -  /2 
eff. 


C,  is  from  Figure  4-15 
h 


Of 


C.  is  from  Figure  4-16 
‘*8 

C,  is  from  Table  4-1 
‘'i 

eff 


=  -  .  043  a  + 
2tt 


-  .54  +  .47  X  .043  U  - 
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\  ( 
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Figure  4-3.  Fwd.  Foil  Zero  Lift  Angle 
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Figure  4-14.  Prototype  Aft  Foil  Lift  and  Hinge  Moment 
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Figure  4-15.  Flap  Hinge  Moment  Derivative,  Cj^  ,  2-Dlmensional 
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5.  UNSTEADY  LOAD  CHARACTERISTICS 


Unsteady  load  effects  present  significant  attenuations  and  rather  large  phase 
shifts.  Reference  5-9  demonstrates  that  these  effects  can  extend  to  wavelengths  of 
significance  to  the  hydrofoil  craft. 

The  hydrofoil  can  present  from  nine  to  twenty  unsteady  load  responses,  depend¬ 
ing  upon  configuration.  The  AG(Ei’'.  fitted  with  flap  lift  control  and  variable  incidence, 
requires  evaluation  of  all  twenty  of  these  responses  for  autopilot  and/or  structural  con¬ 
siderations.  The  unsteady  load  review  of  this  section  is  not  intended  to  be  an  exhaustive 
theoretical  review  of  the  subject  or  to  contribute  to  the  stace-of-the-art.  It  is  orJy 
intended  to  review  the  predictions  available  for  unsteady  load  response  and  to  measure 
the  confidence  level  associated  with  those  predictions  in  order  to  provide  the  best  possi¬ 
ble  representation  for  unsteady  lift  for  the  concurrent  optimal  control  program  develop¬ 
ment  reported  in  Volume  II  of  this  report. 

In  one  area  it  was  necessary  to  proceed  beyond  classic  aircraft  practice  for 
unsteady  load.  In  aircraft  practice,  there  is  no  wave  celerity  and  wavelength  and  fre¬ 
quency  are  inversely  proportionaj,  Hydrofoils  never  encounter  a  zero  wave  celerity 
and  it  was  necessary  to  consider  the  spatial  (circulatory)  and  inertial  terms  separately. 
For  foil  moticn  responses  this  presents  no  problem  because  Theodorsen  has  displayed 
these  two  effects.  In  the  case  of  orbital  motion  only  Leehey  and  Nishayama,  of  the 
literature  reviewed,  have  recognized  the  significance  of  the  difference  between  wave¬ 
length  and  frequency  and  neither  of  them  displayed  the  spatial  and  iner  ial  terms  in  a 
manner  convenient  to  modification  to  a  particular  3-dimensional  case.  Greenberg  was 
therefore  employed  to  derive  Leehey's  result  with  an  intermediate  step  displaying  the 
two  terms  of  interest.  The  final  lift  result,  however,  presents  a  contradiction  with 
Leehey  which  remains  to  be  resolved. 

The  significance  of  the  wavelength/frequency  distinction  is  not  considered  in 
this  report.  Efforts  now  underway  are  directed  to  examining  the  significance  of  every 
unsteady  load  term  to  the  autopilot  design  in  order  to  minimize  the  complex’’ ty  which 
these  effects  introduce  into  the  optimal  control  program. 


The  particular  unsteady  load  responses,  predictions,  and  experimental  results 
considered  In  this  section  are  listed  in  Tables  5-1,  5-2,  and  5-3,  respectively.  The 
sub-section  and  figures  associated  with  each  response  are  listed  in  Table  5-4  and  a 
judgment  of  the  confidence  level  for  each  response  is  presented  in  Tabla  5-5. 

1972  AG  (EH)  flapped  model  test  results  are  not  considered  in  the  main  body  of  this 
report;  they  are  presented  in  an  appendix. 

5. 1  JONES’  FOIL  MOTION  AND  GUST  FUNCTIONS 


Jones’  Eqs.  (29)  and  (30)  of  Reference  5-1  may  be  written 
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where  C(S)  =  C,  /C 
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(5. 1. 1) 


These  are  approximate  expressions,  in  a  form  convenience  to  analysis,  which 
reconcile  conflicting  distinct  theories  at  zero  S.  A  significant  amount  of  smoothing 
is  involved  lor  distances  less  than  3/4  chord,  where  the  starting  vortex  lies  in  the 
foil,  and  perhaps  extending  to  as  much  as  two  chord  lengths. 

The  inverse  wavelength  form  of  (5. 1. 1)  is 

Cfri  -  1  -  165k  _  1.335k 

'■A  =  ®  .045  +  Ik  .3  +  ik 

Cft)  =  1  -  t.361k 

'a  =6  .381  +  Ik 


^ ,  1. 283k 

^^^A  =  3  ^  ■  .54  +  ik 


(5.1.2) 


of  which  the  last  two  are  Jones’  Eq.  (39)  divided  by  the  ratio  oi  the  2-dimensional  to 
the  finite  aspect  ratio  lift  curve  slope.  The  two-dimensional  form  of  Eq.  (5. 1.  2)  is 
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virtually  identical  with  Theodorsen’s  function.  All  three  forms  of  the  equation  present 
the  circulatory  portion  of  the  unsteady  lift;  the  inertial  portion  remains  t3  be  accounted 
for.  Eq.  (5. 1.  2)  is  presented  graphically  on  Figure  5”1  for  reference. 


Jones’  Eq.  (45)  gives  the  lift  transient  response  to  a  gust  as 

G(S)  =  1  -  .236  e  “  *  -  .513  e  “  -  .  171  e  ^ 

'  'A  =  ® 
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where  G  (8)  =  C  /C 
L  L 
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(5.1.3) 


which  are  analytically  convenient  approximations  for  the  rationally  derived  response 
curves  of  his  Figure  9.  Eq.  (5. 1.3)  assume^  y  rectangular  planform  and  8=0 
represents  the  penetration  of  a  sharp-edged  gust  by  the  leading  edge.  Vhe  inverse 
wavelength  equivalent  of  Eq.  (5. 1. 3)  is 
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The  total  lift  response,  including  inertial  effects,  is  included  in  Eq.  (5. 1.4)  but 
only  for  the  case  where  the  reduced  Inverse  wavelength  and  reduced  frequency  are 
identical;  i.e. ,  only  when  k  =  k'.  For  this  particular  case,  Jones’  2-dimensional 
result  is  virtually  identical  with  Sears’  function  and  with  BispTnghoff’s  representation 
(Eq.  5-376  of  Reference  5-2),  except  that  both  of  the  alternative  forms  reference  the 
response  phase  to  the  angle  of  attack  at  mid-chord  rather  than  to  the  angle  of  attack 
at  the  leading  edge. 
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Eq.  (5. 1. 4)  is  presented  graphically  on  Figure  5-2  for  reference  but  is  only  an 
approximation  for  the  hydrofoil  case  where  k  and  k’  cannot  be  equal.  More  appropriate 
forms  for  the  hydrofoil  case  are  considered  in  Section  5. 11. 

5.2  HEAVE  LIFT 

For  a  foil  heaving  in  smooth  'vater,  k  and  k*  are  identically  given  by  Eq.  (3.  7) 
and  Eq.  (3. 1)  degenerates  to 

C,  =  —  k\  -  i-khc  C(k) 

L,  c  c  L 
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(5.2.1) 


The  instantaneous  angle  of  attack  is  -  h  /V  and  Eq.  (5.  2. 1)  is  conveniently 


expressed  as 
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(5.2.2) 


Theodorsen's  prediction  is  aerodynamic  and  a  first  order  ap’.roximation  for  free 
suiface  effect  is  made  in  Eq.  (5.  2. 2)  by  applying  a  classic  free  surface  effect  correc¬ 
tion  to  the  steady-state  lift  curve  slope,  C  . 

L 

O' 

The  only  heave  unsteady  lift  measurements  considered  for  this  report  were  those 
of  Pattison  in  Reference  5-3.  Pattison  compares  theory  and  experiment  in  the  form 
(noting  that  the  heave  sign  convention  here  reverses  that  of  Pattison); 
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Then  Patdson's  results  are  related  *o  Theodorsen's  prediction  In  the  form  of 
Eq.  (5.  2. 2)  by 

C  =  1  TT  k  C 

l.j  Li 


a 


Pattison’s  measured  steady-state  lift  curve  slopes  are  presented  on  Figure  5-3 
where  they  are  compared  with  the  Geising-Smith  p  rediction  of  Pattison's  report,  with 
Feldman's  prediction  of  Reference  5-4,  and  the  Giobs  &  Cox  prediction  of  Reference 
5-5.  The  Feldman  and  Gibbs  &  Cox  represent  the  extreme  range  of  a  number  of 
classic  slope  predictions  available.  The  low  measured  slopes  of  Figure  5-3  reflect 
doubt  on  the  unsteady  lift  results  as  well  and  only  ore  set  of  Pattison's  results  is  pre¬ 
sented  in  this  report. 

Figures  5-4  and  5-5  compare  Pattison's  measured  heave  unsteady  lift  at  one 
chord  submergence  with  several  predicted  heave  responses  and  are  typical  of  the 
results  throughout  the  depth  range  tested  except  that  the  5^  a  effect  at  6.32  Froude 
number  is  not  apparent  at  the  other  depths. 

The  Widnall  predictions  of  Figures  5-4  and  5-5  are  from  Patt. ion's  report. 
Theodorsen's  unsteady  lift  is  compared  with  the  Theodorsen  function  oi.  the  figures  to 
display  the  wavelength  at  which  the  inertial  terms  become  significant;  that  wavelength 


5-5 


would  be  increased  in  head  seas  and  reduced  in  following  seas.  The  Smullin  &  Bender 
predictions  are  from  Figure  10  of  Refenmce  5-6.  The  Stark  program  prediction  is  from 
a  Saab  Aircraft  flutter  analysis  progranri  which  is  not  particularly  well  suited  to  this 
purpose  since  the  lift  results  are  an  intermediate  step;  a  program  modification  and  a 
substantial  hand  calculation  effort  were  required  to  display  those  lift  results.  Stark’s 
Reference  5-7  has  been  received  since  the  derivation  of  the  predictions  of  this  report 
and  would  provide  for  a  more  convenient  evaluation  of  the  aerodynamic  case.  Of  course, 
the  Widnall  and  Bolt,  Beranek,  &  Newman  work  present  the  best  potential  for  the  hydro¬ 
foil  case. 

In  view  of  the  questionable  precision  of  the  data,  the  Pattison  results  present  no 
practical  necessity  for  any  more  sophisticated  prediction  procedure  for  the  heave  lift 
case  than  that  of  Theodorsen’s  Reference  5-8. 

5.3  HEAVE  MOMENT 

For  a  foil  heaving  in  smooth  water,  k  and  k'  are  identically  given  by  £q.  (3.  7) 
and  Eq.  (3. 2)  degenerates  to 
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which  is  conveniently  evaluated  as 
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For  the  steady  state  case,  the  hinge  moment  coefficient  slope  is 
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which  is  a  classic  relationship  and  which  reduces  Eq.  (5. 3.2)  to 
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For  the  2-dimensional  case  there  is  a  further  reduction  to 
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(5.3.3) 


(5.3.4) 


(5.3.5) 


The  amplitude  for  H,  is  unbounded  for  the  general  case  and  carries  little  intuitive 
n 

significance,  therefore  Theodorsen’s  incidence  hinge  moment  is  compared  with 
Pattison’s  data  in  the  2-dimensional  form  of  Eq.  (5.3. 1): 
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and  that  comparison  is  made  with  reference  to  the  quarter-chord  station: 
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Note  that,  for  the  aerodynamic  case,  all  of  the  heave  quarter  chord  moment  is 
produced  by  the  inertial  load  acting  at  the  half-chord.  The  hydrodynamic  case  alters 
this  relationship  only  slightly  by  modifying  the  circulatory  chordwise  pressure  distri¬ 
bution. 


Theodorsen’s  unsteady  moment  is  related  to  Pattison’s  form  of  the  coefficient 
by  (again  noting  a  reversal  for  the  heave  sign  convention): 
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which,  for  Pattison’s  data,  is  related  to  the  quarter-chord  station  by 
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A  typical  comparison  between  the  Pattison  data  and  the  predictions  considered 
here  is  shown  on  Figures  5-5  and  5-6.  The  validity  of  these  comparisons  for  the 
Stark  result  is  doubtful  because  the  Stark  result  was  derived  for  the  AG  (EH)  foil  plan- 
form  and  adjusted  to  the  2-dimePoional  case  only  by  changing  the  reference  lift  curve 
slope.  As  for  the  heave  lift,  the  Pattison  data  will  not  support  a  necessity  for  any 
more  sophisticated  analysis  than  tbit  of  Theodorsen. 
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5.4  PITCH  LIFT 


For  a  foil  pitching  in  smooth  water,  k  and  k’  are  identically  given  by  Eq.  (3.  7) 
and  Eq.  (3, 1)  degenerates  to: 
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Further  reduction  depends  upon  whether  it  Is  the  pitch  or  incidence  lift  case  which 


is  presented: 
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(5.4.3) 

The  difference  between  the  coefficients  for  the  inertial  terms  is  scarcely  signifi¬ 
cant  since  S'/S  is  a  reasonably  good  approximation  for  C  /C  ;  e.g. ,  the  area  ratio 

L,  L 
i  a 

for  the  AG  (EH)  forward  foil  is  .  813  and  the  lift  curve  slope  ratio  is  .  838.  The  hinge 
location,  however,  has  not  been  defined  for  the  incidence  lift  case.  Intuitively  it  seems 
that  the  hinge  should  be  located  on  the  exposed  foil  MAC  for  this  case  rather  than  on 
the  total  foil  MAC  and  it  is  so  defined  in  this  report.  For  this  definition  Eqs.  (5.4.2) 
and  (5.4.3)  are  virtually  identical  and  it  is  the  incidence  lift  form  which  is  presented 
on  Figures  5-8  through  5-11  with  spot  comparisons  with  the  pitch  lift  resuit  at  two 
wavelengths. 

Submergence  enters  Eqs.  (5.4.2)  and  (5,4.3)  as  a  steady-state  lift  curve  slope 


but  the  significance  is  slight.  The  1. 14  MGC  submergence  slope  is  presented  in  the 
figures  with  spot  comparisons  with  the  .57  MGC  submergence  result  at  two  wavelengths. 


The  test  applied  to  Theodorsen's  pitch  lift,  which  is  aerodynamic,  for  the  hydro- 
dynamic  case  is  that  provided  by  O' Neill's  data  of  Reference  5-9  which  presents  the 
incidence  lift  case.  Two  other  predictions  are  compared  with  Theodorsen  and  with 
O'Neill's  data  on  Figures  5-8  through  5-11. 

•  Wldnall's  hydrodynamic  lifting  surface  program.  This  is  not  the  2-dimen- 
sional  prediction  of  Reference  5-9  but  a  later  S-dimensional  prediction 
supplied  by  the  Contracting  Officer 

•  The  Stark  aerodynamic  lifting  surface  program  previously  mentioned. 

The  Smullin  and  Bender  hydrodynamic  lifting  surface  program  is  of  particular 
interest  because  it  provides  flap  effects  but  the  unsteady  pitch  lift  of  Figure  13  of  Ref¬ 
erence  5-6  is  2-dimensional  and  referenced  to  a  qua?  ter-chord  hinge  location;  there¬ 
fore  not  comparable  with  O'Neill's  data.  The  Smullin  and  Bender  result  is  related  to 
O'Neiirs  data  here,  though  remotely,  by  comparison  with  the  2-dimensional  form  of 
Eq.  (5.4. 1)  with  a  quarter-chord  hinge  on  Figures  5-12  and  5-13.  For  this  case  the 
Theodor  sen  equation  reduces  to; 

Cl  /e 

A  -  — ^  =  (1  +  ik)C(k)  + 

L 

or 

=  (1  +  ik)C,M  +  1.208  -  (-.  I  +  i)  forj  =  |  (5.4.4) 

Figures  5-8  through  5-11  do  not  indicate  any  practical  requirement  for  resort 
to  the  added  complexity  of  the  lifting  surface  theory. 

5.  5  PITCH  MOMENT 

For  a  foil  pitching  in  smooth  water,  k  and  k'  are  identically  given  by  Eq.  (3.  7} 
and  Eq.  (3.2)  degenerates  to 
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where  It  must  be  understood  that  the  final  term  is  inferred  in  Eq.  (3.  2) 


All  of  the  parenthetical  terms  containing  the  hinge  location  are  from  2-dimensional 
theory  and  none  have  been  defined  for  the  practical  case.  In  one  instance  the  (~  -  “  j 
term  becomes  identifiable  if  the  circulatory  term  is  written  in  the  form: 
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where  the  measured  's  of  Table  4-2  are  employed  in  this  study. 

It  is  not  so  obvious,  and  in  fact  it  is  unlikely,  that  the  second  occurence  of 
~  )  is  .  Throughout  this  study,  except  for  the  one  use  of  ,  the 

hinge  has  been  related  to  the  total  foil  MAC  for  pitch  lift  and  to  the  exposed  foil  MAC 
for  incidence  lift.  The  chord  reference  stations  —  1/4,  1/2,  and  3/4  —  are  employed 
unchanged.  Actually  for  the  3-dimensional  case  the  fractional  terms  require  redefini¬ 
tion. 


O’Neill’s  data  is  compared  with  Theodorsen,  Widnall,  and  Stark  on  Figures  5-14 
through  5-17  where  the  following  comments  are  in  order. 

Theodorsen’s  prediction  correlates  well  with  O’Neill’s  data  throughout  except  for 
phase  for  the  327r  hinge  on  the  shortest  waves.  The  indication  is  that  the  camber  term 
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of  the  circulatory  terms  or,  more  likely,  the  pitch  acceleration  term  needs  more 
attention.  Note  that  it  is  the  experimental  C,_  of  Section  4  which  correlates  the 

data  of  Figures  5-14  and  5-15.  The  theoretical  C  of  that  section  does  not  produce 

as  good  a  correlation,  inferring  good  consistency  between  the  measured  steady-state 

and  unsteady  hinge  moments.  The  relative  insignificance  of  the  circulatory  term  is 

encouraging  because  elimination  of  this  term  would  eliminate  a  troublesome  cross 

product,  01  C  C(k),  in  the  equations  of  motion. 
a 

The  Widnall  prediction  has  overemphasized  the  steady-state  moment  slope  for 
the  32%  hinge  with  a  consequential  effect  on  the  phase.  The  effect  is  not  so  evident 
for  the  6.2%  hinge  where  C  is  not  critical  but  the  practical  case  is  the  low  , 

preferably  zero,  case.  Note  that  one  disadvantage  of  the  lifting  surface  programs  is 
the  accessability  of  individual  characteristics  for  modification. 

The  Stark  program  also  missed  the  steady-state  hinge  moment  and  is  presented 
on  the  figures  as  the  product  of  gain  and  the  measured  steady-state  hinge  moment. 

The  procedure  is  not  really  valid  for  amplitude,  however,  and  does  not  affect  phase 
so  the  result  is  not  satisfactory. 

Smullin  and  Bender’s  prediction  is  not  available  in  3~dimensions  and  is  compared 
with  Theodorsen's  2-dimensional  result  on  Figures  5-18  and  5-19.  The  aerodynamic 
result  only  indicates  that  Smullin  and  Bender  p/oduce  Theodorsen  at  infinite  submer¬ 
gence. 

For  the  2-dimensional  case,  Theodorsen's  unsteady  moment  about  the  quarter- 
chord  reduces  to 

C  -  ^  k(.375  k  -  i)  (5.5.3) 

having  no  circulatory  contribution.  What  the  lifting  surface  programs  can  provide, 
which  is  not  yet  available  in  analytic  form,  is  the  steady-state  moment,  i.e. ,  the 
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chordwise  load  distribution,  which  must  vanish  relative  to  the  quarter-chord  in  1 

Theodorsen’s  unsteady  lift.  Comparison  of  Smullin  and  Bender’s  Figures  13  and  14  i 

indicate  an  aerodynam.ic  center  shift  of  almost  2\  of  the  chord  forward  at  a  steady- 
state  submergence  of  1/2  chord.  The  effect  on  the  moment  amplitude  on  Figure  5-18 
is  small,  but  there  is  a  significant  effect  on  the  phases  of  Figure  5-19. 

In  summary,  Theodorsen  provides  reasonably  good  correlation  for  O’Neill's 
hinge  moment  data  and  better  correlation  than  is  available  at  present  with  the  lift 
surface  theories,  though  Theodorsen  does  not  provide  the  understanding  of  free  surface 
effect  chat  will  become  available  with  a  fully-developed  lifting  surface  theor>'. 

5.6  FLAP  LIFT 


For  flap  motion  on  a  2-di:nensional  airfoil  in  smooth  water,  k  and  k’  are  identi¬ 
cally  given  by  Eq.  (3.7)  and  Eq.  (3.  1)  degenerates  to 
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(5.6.1) 


For  the  theoretical,  thin  airfoil,  case  do'/d6  is  Theodorsen’s  given  here 

in  Table  3-1.  For  the  practical  case,  the  flap  effectiveness  may  be  taken  from  Toll’s 
Reference  5-10  or  equivalent  or  from  model  tests.  Eq.  (5.  6. 1)  is  aerodynamic  but  a 
first  order  approximation  to  the  hydrodynamic  case  may  be  obtained  by  making  the 
classic  adjustment  to  the  lift  curve  slope. 

Eq.  (5.  6.  1)  is  compared  with  Smullin  and  Bender's  Figure  16  on  Figures  5-20 
and  5-2T.  Also  shown  on  these  figures  is  Theodorsen’s  equation  for  the  AG  (EH)  for¬ 
ward  foil  at  l/2  ehord  submergence: 


.iIUk)c  CB  .  nSkDk  . 

1 


=  (l  +  ^  Al  l  2.06  -"23.221^  +  i  .  1424) 


=  (1  +  i.3165k)C(k)  -  .0611k(k  -  i  6. 14) 


(5.6.2) 


The  AG  (EH)  foil  response  has  been  measured  in  the  towing  tanl^  but  the  results 
are  still  being  analyzed  at  this  v/riting  and  will  be  reviewed  in  an  appendix  to  this 
volume. 

Figiires  5-20  and  5-21  indicate  that  this  response  is  probably  insignificant  to  the 
AG  (EH)  though  significant  for  foils  of  higher  aspect  ratio. 

5.  7  FLAP  HINGE  MOMENT 


For  flap  motion  on  a  2-dirrierisional  airfoil  in  smooth  water,  k  and  k’  are  identi¬ 
cally  given  by  Eq.  (3.7)  and  Eq.  (3.3)  degenerates  to: 
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(5.7.1) 


For  the  theoretical,  thn  airfoil,  case  the  flap  hinge  moment  derivative  C,  is 

n 

Of 

related  to  Theodorsen's  T  and  may  be  taken  from  Table  3-1.  The  practical  case  is 
more  difficult  and  uncertain.  For  this  report  the  2-dimensional  value  of  tliis  derivative 
has  been  evaluated  by  the  methods  of  References  5-10  and  5-11,  which  are  quite  similar 
for  two  sections  of  interest  to  hydrofoils  with  the  results  presented  on  Figure  5-22. 

The  3-dimensional  derivative  may  practically  be  evaluated  as: 
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(5.7.2) 
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which  is  essentially  the  3 '■dimensional  evaluation  of  Perkins  fr  Huge  in  Reference  5-12. 

The  thin  airfoil,  2-dimensional  flap  hinge  moment  derivative  C,  can  be  related 

•■‘5 

to  Theodorsen’s  ’’T”  coefficients  by  writing  the  steady-state  form  of  Eq.  XIX)  of  Ref¬ 
erence  5-8  for  flap  deflection: 

M  2  ^12 
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(5.7.3) 
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which  is  compared  with  the  practical  values  of  Toll  and  the  Datcom  on  Figure  5-23. 

Toll  does  not  give  a  conversion  to  the  3-dimensional  case  for  this  derivative  and 
the  Datcom  conversion  is  not  ver>  convenient.  The  following  conversion  is  adapted 
from  Perkins  and  Hage  and  is  presented  without  theoretical  justification. 

C  =  C  +  %  (  C  -  c  ) 
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Smullin  and  Bender's  flap  moment  response  is  compared  with  Theodorsen’s  on 
Figures  5-24  and  5-25  which  indicate  little  practical  difference  between  the  two  predic¬ 
tions.  The  comparison  is  worsened  when  the  Smiillin  and  Bender  steady-state  deriva¬ 
tives  are  substituted  for  Thcodorsen's  in  the  Theodorsen  equation.  Presumably  this 
is  because  Theodorsen’s  coefficients  are  so  interrelated  that  changes  in  any  coefficient 
reflect  changes  in  others.  It  would  be  desirable  to  identify  all  of  Theodorsen’s  "T” 
coefficients  in  order  to  adapt  his  equations  to  the  practical  case  with  more  confidence. 
The  30^  chord  Theodorsen  result  is  shown  on  Figures  5-24  and  5-25  to  relate  these 
figures  to  the  following  two  figures. 


5-14A 


f  i 

it  w 


i  t 

%  9 


r 

% 


The  Stark  program,  which  was  for  the  AG  (EH)  foil,  is  compared  with  Theodorsen 
on  Figures  5-24  and  5-25.  The  Stark  3-dirnensional  and  Smullin  and  Bender  2-dimen¬ 
sional  results  are  quite  similar  and  both  differ  rather  significantly  from  the  3-diii*eii- 
sional  Theodorsen  result. 

The  AG  (EH)  flap  moment  response  has  been  measured  in  the  towing  tank  but  the 
results  are  still  being  analyzed  at  this  writing  and  '✓ill  be  reviewed  in  an  appendix  to 
this  volume. 

5.  3  FLAP  MOTION  FOIL  HINGE  MOMENT 

For  flap  motion  in  smooth  water,  k  and  k’  are  identically  given  by  Eq.  (3.  7)  and 
Eq.  (3.2)  degenerates  to*. 
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which  reduces  for  the  steady-state  case  to: 
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which  holds  academic  as  well  as  practical  interest  and  which  is  discussed  further  in 
Section  4.  The  numerical  evaluations  considered  here  employ  for  j  )  the 

measured  incidence  lift  of  Table  4-2. 

^L 

Two  forms  of  Eq.  (5.  8.  1)  are  compared  with  Figure  IS  of  Reference  5-6  on 
Figures  5-28  and  5-29;  the  Stark  prograin  output  is  not  available  for  this  response. 
There  are  no  significant  differences  between  any  of  the  responses  considered  and 
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tills  particular  response  may  be  negligible  for  the  practical  cas^,  depending  upon  the 
sensitivity  of  the  closed  loop  system  to  phase  shift. 

This  response  has  been  measured  in  the  towing  tank  on  the  AG  (EH)  model  but 
the  results  are  still  being  analyzed  at  this  writing  and  will  be  reviewed  In  an  appendix 
to  this  volume. 


5.  9  FOIL  HEAVE  FLAP  HINGE  MOMENT 

For  a  foil  heaving  in  smooth  water,  k  and  k’  are  identically  given  by  Eq.  (3.  7) 
and  Eq.  (3.3)  degenerates  to: 
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Eq.  (5.  9. 1)  may  be  written 
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For  the  theoretical,  thin  airfoil  case  che  denominator  of  the  second  term  may  be 
replaced  by  Its  Theodorsen  equivalent  of  Note  2  of  Table  3-1: 
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Theodorsen’s  response  is  compared  with  Figure  12  of  Reference  5-6  on  Figures 

5-30  and  5-31  which  indicate  that  this  will  be  a  troublesome  response.  The  two  sets 

of  results  differ  substantially  in  the  hinge  moment  slope  with  respect  to  foil  angle  of 

attack,  C,  ,  and  in  the  inertial,  quadrature  component.  The  quadrature  component 
h 

O' 

can  be  examined  by  writing  Eq.  (5.  9. 3)  in  the  form 


(5.9.4) 


where  the  first  form  may  be  expected  to  represent  actual  practice  while  the  second 
form  is  Theodorsen's  2-dimensional,  thin  airfoil  result.  The  left  side  of  the  equation 
is  appropriate  for  Smullin  and  Bender’s  aerodynamic  case  but  only  approximate  for  the 
finite  depth  case. 

Smullin  and  Bender’s  predictions  are  plotted  in  the  form  of  Eq.  (5.  9.4)  on  Figure 
5-32.  The  corresponding  phases  are  nearly  in  quadrature  and  inconsequential  to  this 
analysis.  The  Theodorsen  predictions  in  this  form  are  straight  lines  of  slope  deter¬ 
mined  by  the  relative  flap  chord  and  by  the  denominator,  of  Eq.  (5.  9.4),  employed. 
Obviously,  Theodorsen  and  Smullin  &  Bender  differ  substantially  in  quadrature  com¬ 
ponent,  as  well  as  in-phase  component,  and  the  problem  is  going  to  be  particular y 
difficult  for  the  practical  case  where  the  sections  employed  produce  low  and  uncertain 
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No  data  has  been  located  to  test  this  response  and  the  Stark  program  was  not 
employed  for  this  case. 


5-17 


5. 10  FOIL  PITCH  FLAP  HINGE  MOMENT 


For  a  foil  pitching  in  smooth  water,  k  and  k’  are  identically  given  by  Eq.  (3.  7) 
and  Eq.  (3.3)  degenerates  to 


where  the  incidence  lift  case  is  not  being  considered  here. 

Eq.  (5. 10. 1)  is  compared  with  Figure  15  of  Reference  5-6  on  Figures  5-33  and 
5-34  which  indicate  much  the  same  difficulty  associated  with  the  flap  hinge  moment 
response  to  foil  heaving  motion  except  that  there  are  several  more  effects  associated 
with  the  pitching  motion. 

No  data  has  been  located  to  this  response  and  the  Stark  program  was  not  employed 
for  this  case. 

5. 11  ORBITAL  MOTION  LIFT 

In  the  notation  of  this  volume,  Eq.  (12)  of  Greenberg’s  Reference  5.  13  is 
L  w 

-  2n  p|  C(k;  [.Tp(k)  +  J^(k,]  (5.11.1) 

where  the  subscript,  C,  indicates  circulatory  lift. 

The  2tt  coefficient  is  the  2-dimensional  lift  curve  slope  and  Eq.  (5. 11. 1)  t'S 
arbitrarily  generalized  by  writing 


=  C(k)  (5.11.2) 
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where  the  prime  indicates  that  the  lift  phase  reference  is  the  or’.jital  angle  of  attack  at  the 
half-chord  position. 


Greenberg's  Eq.  (BG)  in  the  notation  of  this  volume  may  be  written 


(5.11.3) 


where  the  subscript,  NC,  indicates  non-circulatory  (inertial)  lift. 


Greeiiberg  reduced  this  equation  for  the  case  of  a  zero  wave  celerity  but  a  more 


general  reduction  follows  from  the  relationship: 
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Then  Eq.  (5. 11.  3)  may  be  written 
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This  is  an  inertial  term  and  the  rr  coefficient  is  not  related  to  the  lift  curve  slope 
but  to  a  fluid  volume,  therefore  the  reducl  n  proceeds: 
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Then  the  total  lift  is  given  by 
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=  [jQ(k)  -  i.J^(k)]  C(k)  +  ^  J^(k) 

L 

Oi 


(5.  XL  7) 


For  the  2-dimensional  case  Eq.  (5. 11.  7)  reduces  to 

t 

g'  =  [jo(k)  -  lj^(k)]  C(k)  +  i|-J^(k)  (5.  XI.  8) 

which  is  10  be  compared  vith  Eq.  (8. 26)  of  Leehey’s  Reference  5. 14  and  with  the 
identical  Eq.  (6)  of  Steele’s  Reference  5. 15* 

? 

g'  =  [j(,(k)  -  ij,(k)]  C(k')  +  J^(k)  (5.11.9) 

which  makes  Theodorsen’s  function  a  temporal  rather  than  a  spatial  function. 

For  waves  of  zero  celerity,  Eqs.  (5. 11.  8)  and  (5. 11.  9)  both  reduce  to 

g’  =  [jjj(k)  -  iJ^(k)]  C(k)  +  iJ^(k)  (5.11.10) 

which  is  Eq.  (5-376)  of  Bisplinghoff’s  Reference  5-2  and  which  is  virtually  identical 
numerically  with  Sear’s  function. 

Eqs.  (5. 11.  7),  (5. 11.  8),  and  (5.  11. 10)  may  be  shifted  in  phase  to  reference  the 
orbital  angle  at  the  leading  edge  by 
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and  for  this  reference  phase  Eq.  (5. 11. 10)  is  virtually  identical  numerically  with  Jones’ 
2«dimensionai  ^st  function;  only  Jones  offers  a  3-dimensional  version  of  this  function 
though  in  the  volume  Eq.  (5. 11.  7)  is  employed  with  Jones’  3-dimensional  C(k)  and  with 
the  3-dimensional  hydrodynamic  lift  curve  slope  to  approximate  3-dimensional  and 
free-surfacc  effects. 


O’Neill  offers  an  adjustment  to  the  results  above  which  accounts  for  the  variation 
in  dynamic  pressure  through  the  wave  cycle.  From  Reference  5-9  that  adjustment  has 
the  form 
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(5.  11.12) 


where  C  is  the  steady-state  lift  coefficient.  This  q  effect  increment  was  tested  in 

L  j 

the  preparation  of  this  volume  hut  the  assumption  of  a  fixed  reference  C  should  be 

Li 

reexamined  in  the  future. 


In  deep  water  gravity  waves  the  ratio  of  the  wave  celerity  to  craft  velocity,  for 
head  and  following  waves,  is  given  bv 
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so  the  lift  response,  G,  is  a  rather  complex  function  of  Froude  Number,  depth,  and 
wavelength. 

For  Reference  5-15,  Steele  ran  an  aspect  ratio  6-1/4,  rectangular  foil  under 
head  and  following  seas  at  the  i^cro  lift  angle  of  attack.  Wave  length/height  ratios 
were  very  long,  ranging  from  24  to  226.  Measured  wave  celerities  checked  well  with 
theoretical.  Orbital  velocities  varied  from  .  8 ^  to  7.5^  of  model  speed  introducing 
the  possibility  of  q  variations  to  157c. 

For  this  report,  Steele’s  data  was  reduced  directly  from  his  Tables  2-5  by: 

-  (k  i  (5. 11.14) 


The  following  numerical  values  were  employed 

p  =  1.938 
S  =  25  sq.  in. 

Wadlin’s  free-surface  effect 


k  ^  at  4”  depth 


k  »  at  2”  depth  (5.11.15) 

\ _ w 

Steele’s  data  is  compared  with  Eq.  (5. 11.7),  shifted  in  phase  by  Eq. 

(5. 11.  11),  on  Figures  5-35  and  5-36.  Jones'  gust  function  for  aspect  ratio  6,  which 
Is  an  aerodynamic  result,  is  also  shown  on  these  figures. 
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Leehey  notes  that  Eq.  (5. 11. 10)  approaches  the  value  of  C(k)  for  sufficiently  long 
waves  and  likens  the  orbital  lift  response  to  the  heave  response  on  long  waves.  The 
heave  response  is  therefore  compared  wi^h  Steele's  data  and  with  the  predictions  of 
Figures  5-35  and  5-36.  The  analogy  is  not  good  in  the  wavelengths  of  interest  to  the 
control  system  and  is  theoretically  weak  anyway;  actually,  all  the  unsteady  load  res¬ 
ponses  approach  each  other  on  long  wave©  because  they  all  approach  unity. 

O’NeiU's  form  of  the  q  effect  is  not  applicable  to  Steele’s  data.  The  effect  of  the 
sinusoidal  variation  in  submerg  ace  was  considered  briefly  and  should  be  considered 
in  more  detail  in  the  future.  The  effect  on  the  steady-state  lift  curve  slope  could  pro¬ 
duce  a  substantial  phase  shift  but  only  a  5^  change  in  amplitude.  The  effect  on  the 
instantaneous  orbital  angle  was  not  examined.  In  summary,  no  accountability  was 
found  for  Steele's  following  sea  amplitude  or  head  sea  phase.  The  low  Froude  Number 
for  Steele's  data  should  be  noted. 

The  data  presented  by  Conolly  in  Reference  5-16  was  examined  but  found  incon¬ 
clusive  because  of  data  scatter.  Coaolly's  following  sea  lift  responses  on  24  chord- 
length  waves,  corrected  for  q  effect,  present  attenuations  to  1. 1  with  a  mean  at  about 
.  9.  The  head  sea  responses  range  to  an  attenuation  of  1.  45  with  a  mean  at  about  1.  1 
so  the  data  qualitatively  reverses  the  indication  of  Steele’s  data  but  with  a  low  level  of 
confidence. 

For  Reference  5-9,  O'Neill  ran  the  AG  (EH)  foil  model  under  head  seas  at  two 
speeds  (with  distinct  angles  of  attack),  two  foil  hinge  locations  (insignificant  to  lift 
data),  and  at  a  single  depth.  Wave  length  to  height  ratios  were  very  long,  ranging 
from  30  to  125.  Orbital  velocities  do  not  exceed  1.  5%  of  model  speed  but  the  model 
was  at  prototype  angle  of  attack,  introducing  the  possibility  of  q  effect. 

c 

O'Neill's  data  is  in  terms  of  A  77  and  is  referenced  to  the  wave  peak  at  the 

L  nt. 

mean  geometric  quarter  chord.  O’Neill's  data,  then,  is  considered  here  in  the  form 
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Orbital  motion  acts  on  the  pitch  lift  curve  slope  which  was  not  measured  in  Ref¬ 
erence  5-9.  The  steady-state  lift  curve  slope  employed  here  was  the  2.  75  value  which 
appears  at  1.07  MAC  on  Figure  4-1.  The  data  reduction  for  this  report  proceeded 
from  O’Neill’s  tabulated  values  and  the  attenuations  do  not  correspond  to  O’Neill’s 
gains.  O’Neill’s  q  effect  has  been  subtracted  from,  the  data  and  from  the  Leehey  predic¬ 
tion  appearing  on  his  curves  and  the  correction  tended  to  reduce  the  scatter  in  his  data. 


O’Neill’s  data  is  compared  with  several  predictions  on  Figures  5-39  and  5-40, 

The  comparison  generally  presents  one  group  of  2 -dimensional  predictions  and  another 
group  of  S-dimensional  predictions  spaced  neatly  on  opposite  sides  of  the  data.  Basi¬ 
cally  the  difference  between  the  2-D  and  3-D  predictions  can  be  seen  in  the  Jones  foil 
motion  function  of  Figure  5-1  and  the  Jones  gust  functions  of  Figure  5-2.  O’Neill’s 
data,  then,  simply  indicates  inadequate  accountability  for  planfcrm  effects  in  the 
presence  of  a  free  surface.  The  accountability  for  planform  and  depth  are  obviously 
inadequate  in  Eq.  (5. 11.  7)  since  the  inertial  term  is  assumed  to  be  independent  of 
both.  Figures  5-39  and  5-40  indicate,  by  comparing  Eq.  (5. 11.  7)  and  the  Jones  gust 
function,  that  this  particular  inadequancy  has  significance  in  phase  but  not  in  amplitude. 

It  is  interesting  to  note  that  if  Sinullin  and  Bender’s  aerodynamic  curve  is 
matched  with  the  three  dimensional  curves,  their  hydrodynamic  curve  falls  through 
the  data;  i.  e. ,  the  prospects  for  a  3-dimensional  version  of  the  Smullin  and  Bender 
program  appear  promising. 


In  summary,  the  lift  response  to  the  orbital  angle  still  presents  theoretical  and 
experimental  contradictions  of  considerable  academic  significance.  The  question  of  . 
the  sensitivity  of  the  autopilot  design  to  this  response  therefore  assumes  considerable 
practical  importance. 
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5. 12  ORBITAL  MOTION  MOMENT 

In  the  notation  of  this  volume,  Eq.  (13)  of  Greenberg’s  Reference  5-13  is 


2  vv 

—  2-  pv^rf  [i  - 


P  2  2 

=  -  -  V  C  X  n 
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■^[1  ■  ^2) 
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C  -  -  n 

M  V 

"c 


[2  ■  "  2  (5.12.1) 


If  the  hinge  moment  reference  is  expressed  as  a  chord  station  and  identified  as 
the  foil  h-nge,  "a”  becomes 


H  2 


-  1 


(5.  12.2) 


and  (5. 12.  .1)  becomes 
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It  is  Interesting  to  note  that  for  the  steady-state  case  this  coefficient  vanishes  for 
a  hinge  at  half-chord;  i.e. ,  the  center  of  pressure  for  the  steady-state  circulatory  lift 
is  at  the  half-chord  position. 


Greenberg's  Eq.  (B15)  in  the  notation  of  this  volume  may  be  written 

2  _  w  3  w 
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(5.12,4) 


At  this  point  Greenberg  assumed  a  zero  v/ave  celerity,  which  is  still  aerodynamic 
practice  but  which  cannot  be  the  case  for  gravity  waves.  For  the  hydrodynamic  case 


and  (5. 12.4)  may  be  written 
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From  the  recurrence  relationship  for  Bessel  functions: 
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Referenced  to  a  hinge  chord  station  bv  Eq.  (5. 12. 2): 
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Note  that  in  smooth  water  this  equation  reduces  to 
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which  is  a  pure  couple  or  a  moment  about  the  aerodynamic  center.  Added  to  the  smooth 
water  form  of  the  circulatory  moment,  Eq.  (5. 12.3),  Eq.  (5. 12.  8)  produces  the  clas¬ 
sical  quarter-chord  center  of  pressure  for  the  aerodynamic  section: 
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Eqs.  (5. 12.3)  and  (5.12.7)  give  the  total  moment  lor  the  unsteady  case: 
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The  reduction  is  trivial  but  is  given  in  dotail  to  conserve  the  time  of  the  inter¬ 


ested  reader: 
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With  reference  to  the  aerodynamic  center  the  lift  term  drops: 


(5.12.11) 
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and  for  the  2-dimensional  foil  there  is  a  further  reduction  to 
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Leehey's  Eq.  (8.29’)  in  the  notation  of  this  volume  becomes 
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which  is  identical  with  (5. 12. 13).  The  uncertainty  with  regard  to  the  proper  form  for 
Theodorsen's  function  does  not  enter  the  moment  prediction  then,  except  to  transfer 
the  moment  to  hinge  positions  away  from  the  aerodynamic  center. 

For  zero  wave  celerities  the  moment  about  the  aerodynamic  center  vanishes, 
hence  Bisplinghoff,  e.g. ,  concludes  that  the  lift  in  sinusoidal  gusts  acts  through  the 
aerodynamic  center  (quarter-chord).  Leehey  appears  to  be  the  only  author  who  has 
recognized  the  significance  of  the  wave  celerity  to  the  moment  about  the  aerodynamic 
center  and  this  prediction  must  be  credited  to  him.  Note,  however,  that  Leehey's 
result  is  only  for  the  2-dimensional,  aerodynamic  case  and  that  Eq.  (5. 12. 12)  is  the 
more  general  form. 

There  are  many  pi'oblems  associated  with  Eq.  (5. 12. 12)  for  the  3-dimensional 
case  ?nd  very  little  experimental  guidance.  None  of  the  coefficients,  the  rr  terms,  or 
the  Bessel  function  terms  have  been  examined  for  three  dimensions  or  in  the  presence 
of  a  free  surface.  For  this  volume  only  the  lift  c  :rve  slope  is  modified  to  suit  the 
hydrodynamic,  3-dimensional  case. 
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Only  O'Neill’c  experimental  results  are  examined  for  this  volume.  His  results 
for  the  two  hinge  positions  hav^  been  shifted  to  the  aerodynamic  center  by  the  relation¬ 
ship: 

Si  =  S  -(f  -  S  <5- '2- 15) 

a  a  a 

a.  c. 

where  the  .  0801  C  of  Table  4-2  was  employed  for  ~  -  a.  c.  for  the  .32  MGC  hinge 

S 

position  and  where  that  distance  for  the  .  062  MGC  hinge  was 

T  MGC 

t  -  .062)  -~ 

=  .0801  -  .2.58  X  |;f 

-  -  .1589 

The  data  has  been  compared  with  several  predictions  in  the  form  C^^  /C^ 

O'  a 

a.  c.  w 

on  Figures  5-41  -  5.43.  This  form  of  presentation  has  the  disadvantage  tha'  the  lift 
response  uncertainties  are  brought  into  the  moment  comparison  but  it  has  the  advantage 
of  relating  the  moment  about  the  aerodynamic  center  to  the  moment,  about  any  other 
hinge  position,  produced  by  the  lift;  the  ratio  of  moment/lift  is,  in  fact,  a  distance  as 
a  fraction  of  the  chord. 

The  moment/'lift  ratio  for  the  data  was  derived  directly  from  O’Neill’s  data  in 
the  form 
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The  predictions  presented  are: 

•  The  Leehey  prediction  of  Reference  5-9,  read  from  the  graphs  of  that  refer¬ 
ence.  Presumably  these  are  2-d’‘mensional  predictions  differing  from  the 
Leehey  prediction  below  principally  in  inclusion  of  q  effect 
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•  The  2-D  Leehey  prediction  of  Eqs.  (5. 11.  9)  and  (5. 12. 14)  with  the  k’ A 
ratio  corresponding  to  an  of  5  and  a  one  chord  depth 

•  The  3-D  Leehey  prediction  of  Eqs.  (5. 11.  7)  and  (5. 12. 12)  for  F,  =  5, 
d/c  =  1,  and  C  =  2. 75 

Li 

•  The  Smullin  and  Bender  prediction  of  Figures  7  and  9  of  Reference  5-6. 

This  is  a  2-dimcnsional  prediction  for  F.  =  10  and  d/c  =  1/2. 

n 

The  comparison  of  Figures  5-41  through  5-43  are  not  encouraging.  The  principal 
difficulty  appears  to  lie  in  a  classic  small  difference  between  two  large  numbers,  best 
illustrated  by  example.  For  a  wavelength  of  32  chords,  k  is  .  0982  and  for  an  F^^  of  5 
at  one  chord  depth: 


k’ 


7.07  vT 


=  .451 


Then  Eq.  (5. 12. 14)  becomes 


i  (r  -  0  ■  i 


=  i "  ^  -  ‘-"H 

=  .  708  [.  998  -  .  99796  -  i  .  049o] 

This  is  the  characteristic  arithmetic  and,  in  fact,  for  wavelengths  greater  than 
10  chords  the  moment  slope  may  be  written 
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which  is  sufficiently  small  so  that  practical  considerations;  such  as  3  dimensions, 
sweep,  q  effect,  etc. ;  overwhelm  the  effect  of  Eq.  (5, 12. 14). 
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The  momert  problem  appears  so  forbidding  that  it  would  seem  most  advisable  in 
this  case  to  establish  the  precision  required  first,  by  trending  studies  from  the  auto¬ 
pilot  end  of  the  loop,  before  further  attacks  upon  the  moment  itself. 

In  the  meantime  Figures  5-41  and  5-42  indicate  that  for  the  AG(EH)  a 

a.  c 

of  109(  of  the  lift  can  be  employed.  That  is  in  quadrature  with  the  lift,  presum- 

a.  c. 

ably  leading  on  head  seas  and  lagging  on  following  seas. 

5.  13  ORBITAL  MOTION  FLAP  MOMENT 

Since  steady-state  flap  hinge  moment  is  proportional  to  lift,  the  aerodynamic  foil 
lift  attenuation,  G,  might  be  expected  to  furnish  an  approximation  for  the  flap  moment 
attenuation  in  orbital  motion.  This  presumption  is  supported  by  comparison  of  Jones’ . 
2-dimensional  gust  function  with  the  aerodynamic  flap  moment  response  of  Figure  7  of 
Reference  5-6  as  shown  on  Figures  5-44  and  5-45.  Three  effects  are  not  accounted  for, 
however: 

1.  The  effect  of  wave  celerity  for  which  we  have  only  the  2-D,  aerodynamic 
effect  of  Leehey  and  that  analysis  does  not  include  Rap  load.  The  two 
previous  sections  indicate  poor  agreement  with  experiment,  particularly 
for  foil  moment  which  is  closely  related  with  flap  moment. 

2.  Free  surface  effect  which  further  redistributes  the  chordwise  load.  The 
effect  in  two  dimensions  is  indicated  by  the  two  Smullin  and  Bender  curves 
of  Figures  5-44  and  5-45. 

3.  Planform  effect  which,  presumably,  modifies  the  relative  importance  of  the 
circulatory  and  inertial  terms.  This  effect  is  illustrated  by  the  two  Jones’ 
functions  of  Figures  5-44  and  5-45. 

Those  effects  are  not  available  in  combination  yet  by  an  theoretical  means  and  no 
experimental  data  on  this  effect  has  been  examined.  For  the  present  the  Jones’  3-dimen- 
Lional  gust  function,  operating  on  the  steady-state  hinge  moment  slope,  is  being  employed 
to  represent  the  unsteady  flap  hinge  moment. 

5.  14  SUMMARY 

It  is  to  be  noted  thai  one  of  the  most  significant  limitations  upon  the  confidence 
level  for  moment  response  to  unsteady  motion  is  the  steady-state  response  to  displace- 
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ment.  Errors  here  can  hav^e  a  substantial  effect  upon  the  response  phase  and  a  lesser 
effect  upon  the  amplitude.  This  is  true  of  the  incidence  moment  where  the  steady-state 
moment  is  minimized  to  reduce  control  load  and  of  the  16-series  section  flap  which 
presents  a  vanishingly  small  and  unpredictable  moment  siope.  A  second  limitation,  on 
the  rational  analytical  (non-lifting  surface)  methods  is  a  convenient  accountability  for 
chordwise  load  redistribution  under  a  free  surface.  Both  of  these  limitations  are  signi¬ 
ficant  to  the  steady-state  case  as  well  as  to  the  unsteady  case  and  both  are  currently 
under  review. 

ri'e  measure  of  unsteady  load  prediction  accuracy  is  its  significance  to  the  auto¬ 
pilot  desi.Ti  problem.  Therefore,  current  effort  is  directed  to  approaching  Table  5-5 
from  the  autopilot  end  of  the  control  loop,  employing  the  optimal  control  program  to 
establish  the  unsteady  load  accuracy  required  for  confident  autopilot  design  practice. 
This  will  limit  further  effort  to  improve  the  confidence  levels  of  Table  5-5  to  those 
responses  whero  good  confidence  levels  are  a  practical  necessity.  At  the  same  time, 
the  necessity  for  accounting  for  the  circulatory  and/or  inertial  terms  of  each  response 
and  for  distinguishing  between  wavelength  and  frequency  will  be  examined.  Wherever 
possible,  the  circulatory  portion  of  the  response  will  be  dropped  and  the  inertial  portion 
accounted  for  by  modifying  the  craft  mass  and  inertia. 

5.15  REFERENCES 

5-1  R.  T.  Jones,  ’’The  Unsteady  Lift  of  a  Wing  of  Finite  Aspect  Ratio", 

NACA  Rpt.  681,  1939 

5-2  Bisplinghoff,  Aihley,  &  Halfman,  "Aeroelasticity",  Addison-Wesley, 

1957 

5-3  J.  H.  Pattison,  "Unsteady  Hydrody  imic  Loads  on  a  Two-Dimensional 
Hydrofoil",  NSRDC  Report  3245,  No.  mber  1970 

5-4  J.  Feldman,  "Experimental  Investigation  of  Near-Surface  Hydrodynamic 
Force  Cocfiicients  for  a  Systematic  Series  of  Tee  Hydrofoils,  DTMB 
Series  HF-1",  DTMB  Report  1801,  December  1963 

5-5  "Hydrofoil  Handbook",  Vol.  II,  Gibbs  &  Cox,  Inc.,  1954 


5-33 


5-6  Smullin  &  Bender,  "The  Unsteady  Hydrodynamics  and  Control  of  Hydro¬ 
foils  Near  A  Free  Surface",  Bolt,  Beranek  and  Newman  Report  No.  1970, 
January  1971 

5-7  Stark,  "A  Subsonic  Oscillating-Surface  Theory  for  Wings  with  Partial - 
Span  Controls",  AIAA  Paper  72-61,  January  1972 

5-8  T.  Theodorsen,  "General  Theory  of  Aerodynamic  Instability  and  the 
Mechanism*  of  Flutter",  NACA  Report  496,  1935.  Currently  available 
in  AIAA  Selected  Reprints,  "Aerodynamic  Flutter",  I.  E,  Garrick, 

Editor,  March  1969 

5-9  W.  C.  O’Neill,  "Unsteady  Lift  and  Hinge  Moment  Characteristics  of  the 
AG  (EH)  Main  Foil  and  Strut  Assembly",  NSRDC  Report  2805,  July  1968 

5-10  Toll,  "Summary  of  Late  rial  Control  Research",  NACA  Rpt.  868,  1947 

5-11  "USAF  Stability  and  Control  Methods"  (DATCOM),  Douglas  Aircraft  , 
1965 

5-12  Perkins  &  Hage,  "Airplane  Performance,  Stability,  and  Control",  Wiley, 
1956 

5-13  J.  M.  Greenberg,  "Some  Considerations  on  an  Airfoil  in  an  Oscillating 
Stream",  NACA  TN  1372,  August  1947 

5-14  P.  Leehey,  "The  Hilbert  Problem  for  an  Airfoil  in  Unsteady  Flow", 

DTMB  Report  1077,  January  1957 

5-15  J,  M.  Steele,  "The  Unsteady  Lift  Force  on  a  Restrained  Hydrofoil  in 
Regular  Waves",  NSRDC  Report  3386,  November  1970 

5-16  A.  C.  Conolly,  " Flapped  Hydrofoils  in  Waves,  Subcavitating  Mow", 
General  Dynamics /Con vair  Report  63-032,  May  1963 


5-34 


RESPONSES  I  i  iNClDENCE  MOr4ENT 


TABLE  5-2  UNSTEADY  LOAD  RESPONSE  PREDICTIONS 


1.  THEODORSEN,  REFERENCE  5-8 

2- D,  aerodynamic,  thin  airfoil,  =  0  lift,  foil  and  flap  moment  response  to 

heave,  pitch,  and  flap  motion.  Incorporation  of  practical  steady-state  moment!? 
essential.  Inertial  terms  displayed.  Extend  to  3-D,  hydrod3mamic  with  appro¬ 
priate  steady-state  derivatives;  inertial  accountability  uncertain.  Employ  k’ 
with  inertial  terms  for  C^.  For  3-D  employ  Jones’  C(k). 

2.  JONES,  REFERENCE  5-1 

3- D,  aerodynamic,  =  0  lift  response  to  orbital  motion, 

3.  LEEHEY,  REFERENCE  5-14 

2- D,  aerodynamic  lift  and  foil  moment  response  to  orbital  motion  with  C 
accountability. 

4.  WIDNALL  LIFTING  SURFACE  PROGRAM 

3- D,  hydrodynamic.  2-D  lift  and  moment  response  to  heave  and  pitch  over 
wide  range  of  d/C  and  Fj^  appear  in  Reference  5. 3.  2-D  lift  and  moment  res¬ 
ponse  to  pitch  appear  in  Reference  5.  9.  3-D  AG  (EH)  lift  and  moment  response 
to  pitch  available  but  unpublished. 

5.  SMULLIN  AND  BENDER  LIFTING  SURFACE  PROGRAM 

2-D,  hydrodynamic,  =  0,  with  flap.  All  responses  to  all  motions  pre¬ 
sented  in  Reference  5. 6  for  aerodynamic  and  one  hydrodynamic  case. 

6.  STARK  LIFTING  SURFACE  PROGRAM 


3-D,  aerodynamic,  with  flap  but  without  orbital  motion.  Predictions  of  this 
note  are  intermediate  output  of  flutter  program;  see  Reference  5.  7  for  more 
suitable  program. 
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TABLE  5-3  UNSTEADY  LOAD  DATA  REVIEWED 


1.  PATTISON  -  REFERENCE  5-3 

Lift,  moment,  and  drag  response  to  heave  and  pitch  for  model  between  walls 
(2-D)  without  flap.  Steady-state  lift  curve  slope  is  not  2-D. 

2.  STEELE  -  REFERENCE  5-15 

Lift  response  to  head  and  following  seas,  2  depths,  several  speeds.  Aspect 
ratio  6. 25  foil  without  flap.  Lift  amplificiations  to  1. 5  in  30  chord  length 
following  seas. 

3.  CONOLLY  -  REFERENCE  5-16 

Lift  response  to  head  and  following  seas  for  various  flap  configurations.  Number 
of  wavelengths  measured  too  limited  for  data  scatter  encountered. 

4.  O’NEILL  -  REFERENCE  5-9 


Lift  and  incidence  moment  response  to  incidence  motion  and  to  head  seas  for 
AG  (EH)  fwd  foil /pod/strut.  Good  data  precision  and  steady-state  characteris¬ 
tics. 

NOTE;  This  i'^  not  an  exhaustive  listing  of  the  data  available  in  the  literature.  Time 
did  not  permit  further  data  analysis. 


UNSTEADY  LOAD  COMPARISON 


A 


Theodorsen  Theodorsen  Theodorsen 

Smullin  &  Bender  Smullin  &  Bender  Smullln  &  Bender 

Stark  Stark 


TABLE  5-5  (Continued) 
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Figure  5-4.  Heave  Lift  Attenuation,  2-D,  d/c  =  .  995 
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Figure  5-5.  Heave  Lift  Phase,  2-D,  d/c  =  .995 
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Figure  5-6.  Heave  Moment  Amplitude,  2-D,  d/c  =  .  405 
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Figure  5-12.  Pitch  Lift  Attenuation,  H/c  =  1/4,  2-D 


Figure  5-14.  Pitch  Moment  Slope,  H/MGC  =  .32,  AG(EH)  Fwd.  Foil 


Figure  5-16.  Pitch  Moment  Slope,  H/MGC  =■  .  062,  AG(EH)  Fwd.  Foil 
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Figure  5-18.  Pitch  Moment  Amplitude,  2-D,  H/c  =  1/4 
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Figure  5-19.  Pitch  Moment  Phase,  2-D,  H/c  =  l/4 


Figure  5-20.  Flap  Lift  Attenuation 
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Figure  5^23.  Flap  Hinge  Moment  Derivative  C, 
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Figure  5-27.  Flap  Moment  Phase,  AG  (EH)  Fwd.  Foil 
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Figure  5-28.  Flap  Foil  Hinge  Moment  Attenuat 
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Figure  5-29.  Flap  Foil  Hinge  Moment  Phase 


Figure  5-30.  Heave  Flap  Hinge  Moment  Amplitude 


Figure  5-31.  Heave  Flap  Hinge  Moment  Phase 
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Figure  5-32.  Heave  Flap  Hinge  Moment  Analysis 
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Figure  5-34.  Pitch  Flap  Hinge  Moment  Phase 


Figure  5-35.  Orbital  Motion,  Steele’s  Data,  Attenuation 
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Figure  5-36.  Orbital  Motion,  Steele’s  Data,  Phase 
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gure  5-37.  Orbital  Motion,  Steele’s  Data,  Attenuation 


I 


Figure  5-3S.  Orbital  Motion,  Steele’s  Data,  Phase 
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e  5-40.  Oribtal  Motion,  O’Nelirs  Data,  Phase 


gure  5-41.  Orbital  Motion  Foil  Pitching  Moment  Amplitude 
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Figure  5-42.  Orbital  Motion  Foil  Pitching  Moment  Phase 
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Figure  5-43.  Orbital  Motion  Foil  Pitching  Moment 


lire  5-44.  Orbital  Motion  Flap  Hinge  Moment  Attenuation 


6.  CAVITATION  CHARACTERISTICS 


6.1  INCIDENCE  LIFT  CAVITATION 

The  cavitation  bucket  prediction  employed  in  this  report  is  given  in  some  detail 
because  it  is  not  entirely  classic.  In  particula  ,  it  should  be  noted  that  no  theoretical 
accountability  is  taken  of  the  flow  conditions  at  the  plane  of  symmetry  or  ”kink"  region. 


The  theoretical  spanwise  lift  *1’  ,.»ributions  for  the  AG(EH)  fwd  foil  for  the  pitch 
lift  and  incidence  lift  cases  are  prest  Hed  in  Figure  6. 1.  These  distributions  were 
derived  by  the  methods  of  DeYoung’s  Reference  6.1  where  the  pitch  lift  case  becomes 
a  trivial  special  case  and  the  incidence  lift  case  can  be  presented  as  a  full  chord, 
partial  span  flap.  The  two  lift  cases  produce  significantly,  measurably  different  cav¬ 
itation  buckets.  The  maximum  and  minimum  (at  pod)  lift  coeffieients  on  this  figure 
provide  the  station  and  foil  loading  for  initial  upper  and  lower  surface  cavitation  respec 
tively  for  the  two  lift  cases. 


The  total  velocity  ratio  for  the  section  is 

Av  .  ^ 


J  + 

is  V  -  V- 


where  S  =  pressure  coeffici jnt,  1  -+ 


V  * 

Pq  .  p  eff 


V  =  increment  due  to  thickness,  function  of  chord  statio:i 


V 

Av  =  increment  due  to  camber 

V 


(6.1.1) 


A'^^  “  increment  due  to  angle  of  attack,  function  of 

eff  chord  station 

y  ,  ^  taken  from  Reference  6-2 

V  V 

Upper  and  lower  signs  refer  to  upper  arid  lower  surface  respectively.  The 
effective  design  lift  coefficient  is  a  section  function  and  is  a  little  more  uncertain. 
The  value  employed  here,  839^  of  the  design  lift  coefficient,  is  from  an  unpublished 
analysis  of  the  data  o'  Reference  6"3. 


{ 


The  square  root  of  the  pressure  coefficient  is  a  linear  function  of  the  section 


lift  coefficient: 


V  ^  ^  S  AVa\i  AVa 


s  =  vi  V  ^ 


=  Y  +  AVa  c 


(6.1.2) 


and  a  family  of  such  lines  outlines  the  section  cavitation  bucket  as  shown  on  Figure 

(6.2). 


Cavitation  takes  place  v/hen  the  local  pressure  drops  to  vapor  pressure;  i.e. , 


when: 


p  -  p  p  -  p 

0  0  V 

S  -1  =  - - -  =  a 

q  q 

g  _  ^cav  =  ^cav _ 

•f  V  ^  X2.85  V' 2 


=  .592 


J- —  /  ^cav 
.689^  ^  (S-l) 

,592  1^ 

1  ^  (S-l 


(6.1.3) 


The  most  convenient  way  to  proceed  for  pitch  and  incidence  lift  control  is  to 

calculate  the  cavitation  speed,  V'  ,  for  arbitrary  lift  coefficients  on  the  section  bucket 

K 

of  Figure  (6.2).  Then  the  corresponding  foil  lift  coefficient  is  )  for  the  max- 

imum  or  mini.ium  local  c  ratio  as  appropriate.  Then  the  foil  loading  is  C  q’.  Finally, 

Z  L 

the  entire  process  is  identified  as  having  taken  place  in  the  plane  perpendicular  to  the 

quarter-chord  line  and  the  craft  speed  is  identified  as  V  '  /  cos  A. 

K 

This  process  proceeds  quite  conveniently  in  tabular  manner  under  the  following 
headings: 


1 


--  •  - ;  .  ;  ^1.. 


1.  List  arbitrary  c  ‘s  over  the  range  of  Figure  6-2,  particularly  select 

X 

the  values  at  the  corners  of  the  bucket. 

2.  List  the  corresponding  coefficients,  ^ S 

3.  Compute  V’  from  Eq.  6.1.3 

K 

4.  Compute  C  =  c  ic.  /C,  )  .  from  Figure  6-1 

L  1  i  L  max  or  mm 

5.  Compute  foil  loading,  W/S  =  (W/S)’  =  q'  C 

Xj 

6.  Compute  craft  speed,  V  =  V  ’  /  cos  A 

K  K 

The  result  is  the  foil  cavitation  bucket  of  Figure  6-3. 

Alternatively,  the  relationships  of  Eq.  (6.1.3)  may  be  expressed  as: 

S  -1  -  c  • 

Cf  +  c  )  ^  -1  =  0 ' 

-  -T" 

‘'y  +  c  =  #»!  +  a’ 


-  C  = 


^  +  a^  -'f 


W  W 


,  ^  Kl  +ct’  -Y 

L  ±  Ava  /c 


=  q-CL  = 


V  \  C 


q  r - - 

±lla7v\ 


^  Fl 


But  q'  in  the  plane  perpendicular  to  the  quarter  chord  is: 


q’  =  q  cos  A  =  2  x  2.85  cos  A  V 


-  1.  425  cos^  A  P  V 

K 


and  (G.1.4)  may  bo  written 


(6.1,4) 


(6.1.5) 


6-3 


where  fhe  (+)  sign  goes  with  upper  surface  velocity  distributions  and  with  the  (c  /C  ) 

t  Lj' 
ms' 

and  the  (-)  sign  goes  with  lower  surface  velocity  distributions  and  with  the 
min 

This  equation  provides  the  variation  of  incipient  cavitation  speed  with  foil  loading 
for  any  chord  station. 


The  effect  of  cav  tation  on  foil  lift  as  measured  in  the  Grumman  Whirling  Tank  is 
shown  on  Figure  6-4.  Figure  6-4  was  measured  with  the  AG(EH)  fwd  foil  model ; 
the  PGH  foil  model  results  were  similar  though  fewer  speeds  were  tested.  The  PGH 
model  exhibited  the  zero  lift  angle  shift  at  higher  speeds  shown  on  Figure  6-4.  The 
source  of  this  characteristic  has  not  yet  been  identified.  The  central  and  lower  region 
of  these  characteristics  has  not  yet  been  adequately  mapped. 


The  significant  characteristics  of  Figure  6-4  are  the  initiation  of  the  cavitation 
lift  effect,  which  is  important  to  control  considerations,  and  the  lift  "plateau".  For  the 
astject  ratio  6  PGH  foil  model  the  lift  "plateau"  was  a  well  defined  peak  lift  which  was 
not  again  attained  for  many  degrees  angle  of  attack.  That  peak  lift  coefficient  could 
confidently  be  employed  as  a  limit  load  for  the  PGH.  There  appear  to  be  no  significant 
lift  reversals  on  the  AG(EH)  curves,  but  in  consideration  of  tbe  cavitation  drag,  it 
appears  likely  that  the  plateau  would  still  present  a  limit  load  condition. 

The  effect  of  cavitation  on  foil  drag  is  presented  on  Fij;ure  6-5  which  is  character¬ 
istic.  The  extreme  drags  associated  with  cavitation  insure  that  operation  will  not  extend 
much  beyond  the  initiation  of  this  effect  under  any  conditions. 

Figures  6-4  and  6-5  are  summarized  and  compared  with  the  theoretical  cavitation 


inception  bucket  on  Figure  6-6.  Tnere  appears  to  be  no  difficulty  with  operation  to  the 
effective  cavitation  boundary.  The  lift  effect  boundary  marks  the  end  of  the  linear 
portion  of  the  lift  curve  and  is  sign-'H^runt  to  the  control  problem,  the  drag  boundary 
limits  normal  operation,  and  the  lift  plateau  marks  the  probable  limit  load.  It  is  to 
be  noted  that  there  is  apparently  no  significant  effect,  except  visible  caviation ,  assoc¬ 
iated  with  the  cavitation  inception  but  that  there  is  at  present  no  way  to  predict  the 
position  of  the  effective  cavitation  boundaries. 

It  is  to  be  recognized  that  there  are  objections  to  the  benefit  taken  of  sweep  in 

Eq.  (6.1.6)  and  that  Reference  6-4  presents  conclusions  in  contradiction  of  Figure  6-6 

for  an  identical  planform.  The  derivation  of  Eq.  (6.1.6)  is  essentially  that  of  Reference 

6-5  and  the  method  has  served  well  in  predicting  the  C for  highly  swept  airfoils. 

max 

Appendix  D  of  this  report  draws  support  for  Figure  6-6  from  the  same  data  from  w’hich 
Spangler  drew  a  contradictory  conclusion  in  Reference  6-4.  Nevertheless,  more  exper¬ 
imental  data  directed  to  the  question  of  sweep  effect  is  required. 


6.2  FLAP  UFT  CAVITATION 


In  Reference  6-6,  Allen  presents  the  effect  of  the  flap  on  the  velocity  distribution 
over  the  section  as  an  increase  in  the  increment  due  to  angle  of  attack  plus  an  increment 
in  the  basic  distribution  similar  to  that  associated  with  camber.  He  presents  the  flap 
basic  loading  distribution  in  his  Table  III  for  various  flap  configurations  and  displacement 
ranges.  Allen’s  parameter  is  employed  here  with  the  designation: 


Av 


(6.2.1) 


Allen’s  chordwise  stations  are  somewhat  awkward  to  use,  so  his  distribution  for 
the  209f  chord  flap  at  small  displacements  is  presented  graphically  on  Figure  6-7  for 
convenience. 


Allen  tabulates  the  centroids  for  the  flap  basic  load  distributions,  relative  to  the 
section  aerodynamic  center,  as  G  in  his  Table  IV.  Therefore: 

I 


6-5 


G  c,  =  c 

^  % 

=  d  c 


a.c.  5 


'»5  ^ 


The  flap  angle  is  related  to  the  flap  lift  by 


c ,  =  c,  ^  , 

d6  * 


c  ^ 
d6 


and  Eq.  (6.2.2)  may  be  written: 


G  c ,  do  c 


c  =  c 
b  5 


Gc  do  F 

^adl 


It  i£>  convenient  to  define  the  parameter: 


and  write  (6. 2.4)  as 


C  =  Ch  /  G  ^ 
5  a  d6 


(6.2.2) 


(6.2.3) 


(6.2.4) 


(6.2.5) 


(6.2.6) 


where  r  serves  to  proportion  the  total  flap  lift  between  that  part  producing  the  basic 
load  distribution  ana  that  part  producing  additional  (type)  load  distribution. 

It  must  be  noted  that  the  parameter  C  presents  several  uncertainties  and  thus 


6-6 


m  ^ 

H 


(1 


i  > 


n 

n 


«  It 


n 


becomes  critical  to  the  prediction  of  flapped  foil  cavitation  chai  acteristics.  The  foU 
moment  derivitive  has  already  introduced  uncertainties  into  the  steady-state  incidence 
moment  at  Eq.  (4.  7.3)  and  into  the  unsteady  lift  response  at  Eq.  (5.8.2).  There  is 
some  uncertainty  with  respect  to  all  four  components  of  C  in  2-dimensions  and  a  great 
de^i  of  uncertainty  with  respect  to  their  application  in  3 -dimensions.  For  the  purpose 
of  this  note,  and  pending  some  experimental  evidence  for  flapped  foil  cavitation,  the 
following  evaluations  are  employed: 


c 


H 


6 


-  the  2-dimensional,  thin  airfoil  theory  value. 


G  -  Allen's  value  is  employed 

c  -  the  theoretical,  2-dimensional  value,  2tt,  is  employed 

a 

da  /d6  -  the  practical  value,  .47  for  20%  chord,  is  employed 


The  remainder  of  the  lift  coefficient  due  to  flap  produces  additional  type  loading; 


c  -  c  -  c 
\ 


(C.2.7) 


=  (1-C)c^ 

Now  Eq.  (6.1.1)  can  be  written  in  the  more  generri  form: 

I  „  V  Av  Av  c  c  c  c  ^ 

p  =  YtTt  <T)f  ^'b±(  V 

eff 


Av  Av 

=  Y  +  (^“)  c  +  (c  +  c  )  — = 

X  a 


-  '  V  'f  i. 

b 


(6.2.8) 


where  c  designates  the  lift  coefficient  due  to  angle  of  attack  and  w-here  the  total  lift 

Xj 

X 

coefficient  will  be  designated  without  subscript;  i.e. 


c=c  +c  =c  -:c  +(c)+c 

i  i  '  i  . 

X  F  w  ^0  1  F 


(6.2.9) 


It  is  convenient  to  express  Eq.  (6.  2.  8)  in  terms  of  the  lift  c  and  total  lift: 

JC 

X 
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S=Y+(^)  Cc,  i  [c  +(1-C)c  ]-^ 

F  ^F 


=  yi(-^)  C(c  -c  )Mc^  )]Y 

F  X  X  X 


Av  Av  ^v„  Av„ 

"i'  +  r  (^  )  c  ■i'  r  {—  )  c  +  C  — ^  c  +  (1  “C  )  — -  c 

•  -  ^  '  V  '  i  '  V  '  it  -■  ^  V  V  f 

F  X  F  X 


«vt  cc^-(f )  ]C^  ^  [c-^  -C(^)_-^]  c 


AV. 


Vn  -  ,  AV  ,  _  ft's 


=  Yi  cCV-  -<  v>  ^  ^  Ce[^-^-  '(f  )]-;^]c 


(6.2.10) 


The  parameter  uj  is  defined  for  convenience: 

y>  C  [  Y  “  ^  Y  ' 

F 


Then  Eq.  (6. 2. 10)  becomes 


I  -  AVo 

^S=^i  oiC  4-(a3--^  )  c^ 


(6.2.11) 


(6.2.12) 


For  this  more  general  case,  Eq.  (6.1.4)  becomes: 


S  -1  =a’ 

AVj,  2 

[w+ujC^  +  (rjD-  -^  )  C^]  -l=a' 

X 


V  +  (i)  c  + 

X 


^  (  W  -  -  ^1  +  CT^ 


(6.2.13) 


At  this  point,  it  is  necessary  to  provide  for  orbital  motion  and  to  distinguish  the 
lift  curve  slopes  upon  which  the  incidence  control,  flap  control,  and  orbital  motion  act: 
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Then  the  reduction  of  Eq.  (6.2.13)  proceeds: 


which  may  be  written 


y^r  ) 

;  Lp  w 


(6.2.16) 


(6.2.17) 


The  parameter  ^  is  defined  for  convenience; 


i  -1 


(c  ./C, ) 


i  L' 


F 


?  - 
* 


»  >1 


(6.2.18) 


Then  Eq.  (6.2.17)  will  reduce  to 


i  \  .  ,  , 


Av, 


§.)c;  ^  )  c  ’ 

V  Of  L, 


V  i  L 


w 


L  j,  "  w 


Av, 


W 

s 


=  q-C  =[i(jrT?^.)q'.^)  .(,,.^5.)  (I)  .(,„  +  ^i  )C^'  q'l/(^- 

^  ref  w 


(6.2.19) 


W 

where  (  -  )  is  related  to  incidence  angle  by 
^  ref 

W. 


(-)  =(C  i  +  C.  )q  =  C;  q’ 

ref  h  "^0  ^ 

/  2  ,  2 

C  ’  =  (C  i  +  C  )/  cos  A=  C  /cos  /\ 
L  L.  L-  L  . 

1  0  ref 


^  L  w 

implying  that  C’  =  a  g 

w  2  ,  V 

cos  A 


(6.2.20) 


For  the  case  where  the  incidence  and  flap  lift  curves  slopes  are  the  same,  as  for 
the  AG(EH),  Eq.  (6.2.19)  reduces  to: 


•t'.'  --i 

w 


(6.2.21) 


Eq.  (6.2.21)  provides  the  variation  of  incipient  cavitation  foil  loading  with  speed 
for  any  chord  station  with  flap  lift  control  on  an  arbitrary  incidence  foil  loading.  The 
parameters  required  for  the  Au{EH)  models  for  the  stations  of  interest  are  given  in 
Table  6-1. 

Eq.  (6.2.21)  is  not  as  awkward  nuniorically  as  algerbraically.  For  example,  the 
foil  loading  for  upper  surface  hinge  line  cavitation  in  smooth  water  with  the  incidence 
fixed  at  1435  psf  foil  loading  at  50  knots  is 


=  1975 


(6.2.22) 


For  this  particular  case  the  leading  edge  is  more  restricted  at  1470  psf. 
6.  3  FLAP  OPERATION  MODES 


Where  there  is  '  ariable  incidence  supplied  with  the  flap,  as  for  the  AG(EH),  there 
is  the  question  of  the  optimum  schedule  for  the  incidence.  The  first  schedule  examined 
here  was  that  for  the  design  foil  loading;  i.e. ,  the  incidence  is  varied  with  the  speed  so 
that  the  flap  is  always  neutral  for  Ig  operation.  The  corresponding  cavitation  bucket 
is  shown  on  Figure  6.8. 

The  speed  for  incipient  cavitation  for  the  basic  incidence  control  bucket  is  the 
corss-over  point  when  comparing  the  incidence  control  and  flap  buckets  for  the  fixed 
incidence  loading  mode  of  operation.  Above  this  speed,  the  upper  surface  leading 
edge  boundary  is  expanded  because  the  hinge  line  is  carrying  a  greater  share  of  the  load. 
Below  this  point,  the  leading  edge  boundary  is  reduced  because  the  flap  is  not  as  effective 
as  incidence  in  unloading  the  leading  edge.  Since  operation  at  these  lower  speeds  is  out¬ 
side  the  bucket  in  either  case,  it  is  not  clear  whether  this  is  a  disadvantage  or  not. 

A  second  possible  mode  of  operation  is  to  employ  the  flap  as  a  cavitation  control 
device,  i.e, ,  the  incidence  can  be  scheduled  so  that,  at  any  speed,  when  cavitation 
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occurs,  it  will  occur  at  the  leading  edge  and  hinge  line  simultaneously.  On  Figure  6-8, 
for  example,  it  is  evident  that  the  best  possible  boundary  at  55  knots  occurs  when  the 
incidence  is  set  to  carry  a  1435  psf  foil  loading. 

The  simultaneous  solution  for  Eq.  (6.2.21)  written  for  the  leading  edge  and  hinge 
line  will  yield  the  foil  loading  and  reference  foil  loading  for  this  optimum  bucket  and  the 
result  for  the  AG(EH)  fwd.  foil  is  shown  on  Figure  6-9.  Note,  however,  that  Eq. 

(6. 2. 21)  may  be  written  in  the  form: 


-i. )  s 


W  I -  a’ 

-)  =(VrT7  -  )  --  -7c  7 

ref  ;  L 


;  )  C  +  u)  =  iJTTP’- T)  /  (c  /C  > 


)C^ +  ,,[{€ J  ]  =  -Y)/(c  /C^) 

^  i  0  ^  Y 

A  Vo  Vl+ri*’  -  V 

.,.,c 

1  Jt  L  p  0 


(6.3.1) 


It  follows  that  the  "optimum”  incidence  angle  is  Lndependent  of  speed  except  for 
the  variation  of  the  cavitation  number  with  speed  and,  in  fact,  the  foil  loading  variation 
for  a  constant  lift  coefficient  is  practically  indistinguishable  from  the  optimum,  as 
shown  on  Figure  6-9. 

The  .167  C  of  Figure  6-9  requires  flap  at  50  knots,  but  it  would  be  desirable 
Vef 

to  carry  neutral  or  negative  flap  at  high  speed  in  order  to  reduce  the  flap  required  at 

low  speed.  For  this  reason  Figure  6-10  was  prepared  to  show  the  effect  of  an  increased 

C  on  the  cavitation  bucket.  The  C  selected  for  the  comparison  was  .202, 
ref  ref 

providing  the  design  1435  ai  CO  knots  (trim  flap). 

6.4  ORBITA  L  MOTION  i'JFFECT 

For  the  flapped  foil  operating  at  optimum  incidence  angle,  both  cycles  of  the 
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orbital  angle  in  seas  must  reduce  the  right  boimdary  of  the  cavitation  bucket  because 
they  both  produce  non-optimum  angles  of  attack.  The  ’’optimum"  incidence  angle  is 
defined  for  the  right  boundary  of  the  bucket  without  regard  to  the  incidence  effect  upon 
the  upper-mid-chord  and  lower  surface  bound"iries;  therefore  the  -rbital  angle  does 
not  necessarily  restrict  these  two  boundaries.  This  inefficient  method  for  neutralyzing 
the  orbital  angle  is  the  most  distinctive  hydrodynamic  characteristic  of  the  flap  lift 
control  system  and  is  the  characteristic  which  produces  the  only  significant  reservation 
about  the  system.  It  is  a  very  difficult  characteristic  to  judge  in  real  seas  and,  in  fact, 
quantitative  judgement  does  not  yet  exist. 

The  effect  of  orbital  motion  on  the  flap  lift  control  cavitation  bucket  is  illustrated 
on  Figure  6.11  which  was  derived  for  the  AG(EH)  fwd  foil  from  Eq.  (6.2.21)  and  Tables 
4.1,  4.3,  and  6.1.  The  10  X  200  ft.  wave  train  was  arbitrarily  selected  for  illustration. 
The  classic  orbital  velocity  for  a  10-ft.  depth  was  assumed. 

Platforming  (Ig)  the  10  X  200  ft.  train  of  regular  waves,  the  flapped  foil  cavitation 
bucket  cycles  from  the  smooth  water  bucket  through  each  of  the  extreme  buckets  alter¬ 
nately.  Flying  through  such  a  wave  train  at  50  knots,  the  cavitation  pattern  on  the  foil 
would  cycle  from  clear  through  upper  surface  leading  edge  and  hinge  line  cavitation 
alternately.  In  heavier  seas,  the  margin  to  the  cavitation  boundary  is  the  lift  margin 
available  for  the  normal  accelerations  required  to  negotitate  the  seas  before  cavitation 
takes  place. 

As  for  the  incidence  lift  control,  the  exisience  of  cavitation  is  not  significant,  it 
is  the  lift  and  drag  effect  boundaries  which  are  significant  and  the-e  boundaries  are  not 
predictable.  Towing  tank  cavitation  tests  of  tne  AG(EH)  flapped  foil  conducted  in  connec¬ 
tion  with  this  study  were  intended  to  measure  the  cavitation  eflect  boundary  for  several 
angles  of  attack  and  the  results  are  discussed  in  Appendix  E  of  this  report. 

The  incidence  lift  control  bucket  is  shown  in  Figure  6.11  for  comparison.  This 
bucket  is  not  affected  by  seas;  the  cavitation  bucket  is  a  function  only  of  section  geometry. 
The  comparison  of  the  incidence  lift  bucket  and  smooth  water  flap  bucket  on  Figure  6.11 
indicates  that  in  incidence  lift  control  system  employing  a  flap  only  for  cavitation  control 
would  make  the;  most  effective  possible  lift  control  system. 


The  effect  of  increasing  the  fixed  incidence  angle  to  reduce  the  flap  deflections 
required  is  illustrated  on  Figure  6-12.  If  the  incidence  angle  is  set  to  provide  a  trimmed 
flap  at  50  knots,  a  significant  lift  and  drag  effect  on  the  positive  orbital  angles  is  in¬ 
dicated  by  this  figure.  Further,  there  are  only  2  degrees  of  flap  deflection  to  be  gained 
with  the  increased  angle. 

Figure  6.13  presents  the  variation  of  flap  angle  for  1  +  1/8  g  with  speed  for  the 
two  fixed  incidence  angles  considered.  Thcoe  curves  assume  a  .  47  flap  effectiveness 
throughout.  The  curves  and  the  flap  cavitation  bucket,  through  C ,  are  both  adversely 
affected  by  separation.  The  aerodynamic  literature  anticipates  such  separation  at  flap 
angles  of  5  -  20  degrees  and  the  hydrodynamic  case  is  no  more  predictable,  thus  this 
question  also  Introduces  a  significant  undertainty  into  the  characteristics  of  the  flapped 
foil.  The  AG  (EH)  flapped  foil  towing  tank  model  tests  were  also  directed  to  this 
questions  but  it  must  be  noted  that  this  separation,  even  more  than  that  associated  with 
angle  of  attack,  is  strongly  influenced  by  Reynolds  Number. 
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TABLE  6-1 


CAVITATION  PARAMETERS 
16  -  (.353)  08  SECTION 


STATION 

. 0125c 

.  025c 

.6c 

.7c 

.8c 

v/V 

1.025 

1.051 

1.092 

1.087 

1.067 

AvA^ 

.073 

c 

eff 

.293 

Av/V 

1.346 

.970 

.131 

.103 

.076 

Y 

.704 

1.346 

.840 

1.262 

1.127 

1.  057 

1.130 

1.044 

1.118 

1.016 

(Av/V)p 

.015 

.0275 

.2825 

.4025 

1.1 

c 

•  ’  =  =  545 

Gc  -.397x2t:x.47 

a 

a 

1  -£ 

•  .455 

U) 

.726 

.514 

-.0824 

-.1635 

-.559  j 

COS  A 

.817 

2 

cos  A 

.668 

q' 

P  cos^  A  =  2.84  X  .068  =  1.895 

2 

p 

cav 

-P^  +  pgh  =  2116-72  +  64x8.5  =  2588 

des.  W/S 

1435 

_ a  “  u) 

V 

.620 

.456 

.2134 

.2665 

.635 
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Table  6-1 


STATION 

.0125c  1  .025c  1  .6c 

.7c 

(1) 

a 

1.16S 

.880 

/  c  \ 

Wl 

1 

1.310 

.710 

1.310 

.710 

(c  /C  )  0/ 
i  ^ 

.891 

1.240 

"l 

O' 

2.75 

C  ^  /cos  A 

3.37 

O' 

a’ 

P  /q'  =  2588/1.895  =  1365/  V..^ 

cav  K  K 

. Ill  (Prototype  from  Table  4. 1 

For  dual  numbers  L)\  uvper  number  is  for  upper  surface,  lower  number  is  for 
b 

lower  surface. 
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Figure  6-1.  Spanwise  Lift  Distribution 
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Figure  6-2.  Pressure  Coefficient  Distribution 
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Figure  6-3.  Foil  Cavitation  Buckets 
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Figure  6-9.  Optimum  Flap  Cavitation  Bucket 
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Figure  6-11.  Flapped  Cjivitation  Bucket  in  Seas,  C  =  •  167 
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Figure  6-12.  Flapped  Cavitation  Bucket  in  Seas,  C  =  .202 
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Figure  6-i3.  nap  Deflection  vs.  Speed 


APPENDIX  A 


DOPY  LIFT  AERODYNAMIC  CENTER 

For  foil  aspect  ratios  greater  than  or  equal  to  4,  the  body  aerodynamic  center  is 
located  on  the  exposed  foil  root  chord  by  Eq.  (D6)  of  Reference  4-1  which  may  be  written 


a  +4^112 _ tanAf(5)  AorA'>4 

c'  4  4  l-D(l-i,)  A  or  A  > 

r  r 


(Al) 


D 

where  f  (-)  is 
b 


1  b  D/b  D^  2 


D 

b 


(A2) 


Note  that  Eq.  (Al)  can  be  referred  to  the  sweep  of  any  chord  station  on  the  foil: 


B(W)  ^  1  ,  A  1  +> 


4  4 


-  (tanA  /  c  +  i.x  \  f(5) 


(A3) 


and,  in  particular,  if  the  chord  line  at  x/c  has  no  sweep: 


^B(W)  ^  1 

c'  4 
r 


A  or  A’  >  4 


(A4) 


Thus  in  the  general  case  given  by  Eq.  (Al)  the  body  lift  aerodynamic  center  is  a 
function  of  foil  aspect  ratio  but  for  the  particular  fan  ily  of  foils  laid  out  on  a  straight 
x/c  chord  line  the  body  lift  aerodynamic  center  is  independent  of  aspect  ratio.  It  is  this 
particular  family  of  planforms,  significant  only  to  the  interpolation  procedure  below, 
which  is  referred  to  in  Reference  4-1  as  having  body  lift  aerodynamic  centers  inde¬ 
pendent  of  aspect  ratio. 


The  total  foil  taper  ratio  is  related  to  the  exposed  foil  taper  ratio  by 


A-1 


(1 

X  - - ^ - 

1-5  X. 

b 

SO  that  in  terms  of  the  exposed  foil  geometry  Eq.  (A4)  would  be  written 


1  .  1  -  ..D . . 

=  ^7^-  f<b>  AorA'>4 

r  A=o  1--  (Ah) 

and  this  is  the  form  presented  on  Chart  16  of  Reference  4-12.  This  is  implied  in  the 
text  for  Reference  4-1.  it  is  specified  for  aspect  ratio  in  Reference  4-12  but  there  is 
a  contrary  implication  there  in  the  use  of  ^be  symbol  "  ^  "  where  "  >  would  be  more 
suitable. 

For  the  aspect  ratios  less  than  about  4  the  body  lift  aerodynamic  center  must  be 
interpolated  between  the  value  given  by  (A4)  or  (A6)  and  the  theoretical  location  for  zero 
aspect  ratio,  given  in  Figure  4. 3. 2.1  -  10b  of  Reference  4-12  as 


t  A-  (?.  +X')tanA  j  A’  (1  +  V)  tan  A  «1 

r 

=  .5  ;  j  A'  (1  +  >.')  tan  A  >  1  (A7) 

For  a  foil  having  a  straight  chord  line  x/c,  the  leading  edge  sweep  is  given  by 


L.E.  =  i 


A  =  0  1  +  X 


so  that  it  may  be  more  convenient  to  employ  Eq.  (A7)  in  the  form 


A  =0  A  =0 


:i  -  ^')>i 


where  any  foil  has  a  straight  chord  line  at 


A  =  o 


1  ^  A 
4  4 


1  +  X 
1  -  X 


tan  ^ 


(AlO) 


The  required  interpolation  nniay  be  guided  by  the  aerodynamic  center  for  the 
root  chord  loading  for  the  DA)  =  0  case;  i.e.,  for  the  foil  alone.  Only  the  root  chord 
loadings  of  Reference  4-1,  presented  here  on  Figure  A2*  were  available  for  this  report 
though  Reference  4-1  reports  that  the  subject  is  treated  in  some  detail  in  Reference 
4-13. 


For  the  AG(EH)  foils  the  chord  line  of  zero  sweep  is 


X  .708  =  1.234 

A  =0 


(All) 


,D.  . 


The  functio:i  f  (g  )  is  presented  graphically  on  Figure  Al.  Referring  to  Figures 
1-2,  1-3,  and  Al,  Eq.  (A6)  becomes  for  the  AG(EH)  fwd  and  aft  foils: 


r 

The  parameter  of  Eq.  (A9)  becomes: 


A  or  A’  <'  4 


(A12) 


FWD  (")  (  1  -  X')  =  1.234  X  .671  =  .828 

^  A=  o 


AFT  (-  )  (1  -  V)  =  1.234  X  .  661  =  .816 

^  A  =  o 

and  the  zero  aspect  ratio  body  lift  aerodynamic  centers  become 


(A13) 


r  -B(W)  1  = 

c  ' 
r 

r  -B(W)  1 

«  »  ^ 


.414 


=  .408  A  =  0 


(A14) 
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The  interpolations  required  are  shown  on  Figure  A2  where  the  final  aerodynamic 
centers  selecter  Vvcre 

FWD  ’‘b(W)  =  .371 
c  ’ 


AFT  _Bim  =  .385 


(A15) 


which  are  located  on  the  total  foil  mean  aerodynamic  chord  at 

1.065  ^.371  X  12.12  -  4.93  1.235 

FWD:  a.  c. 


=  .1325 


p(xl  9.33  9.33 

AFT:a.c.  ,  =  .994  +  .385  x  6.07  -  2.545  .784  _ 

pod - ^ — - -  .  --  -  •  tO'3 


4.81 
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Figure  Al. 
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Figure  A2.  Body  Lift  Aerodynami;;  Center,  AG  (EH)  Foils 
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APPENDIX  B 


AG(EH)  FWD  FOIL  DATA  CROSSPLQTS 

The  forward  foil  data  of  Fij^ures  1-8  through  1-17  of  Reference  4-4  arc  presented 
as  lift,  drag,  and  moment  curves  over  a  range  of  incidence  angles  for  five  pod  pitch 
angles.  That  data  is  cross-plotted  on  Figures  B1-B6  to  examine  the  lift  and  moment 
slopes  as  a  function  of  pitch  angle  for  various  incidence  angles.  The  cross-plots  are 
smoothed  to  the  extent  that  not  all  of  the  data  was  measured  at  integral  incidence  angles 
and  the  coefficients  of  the  cross-plots  are  taken  from  the  curves  of  Reference  4-4 
where  necessary  and  in  most  instances;  some  of  the  moment  curves  of  Reference 
4-4  were  reinterpreted  for  these  cross-plots.  The  cross-plots  do  not  include  the 
cavitated  points  of  Reference  4-4. 

The  predicted  lift  curve  slopes  shown  on  the  cross-plots  are  the  Panr'V.#^nkov 
slopes  of  Figure  4-4  which  are  in  reasonable  agreement  with  the  data  throughout  the 
speed  and  depth  ranges.  It  is  possible  that  the  50-knot  pitch  lift  curve  slopes  of 
Reference  4-4  were  measured  in  the  vicinity  of  the  lower  end  of  the  linear  portion  of 
the  lift  curve  where  there  is  a  characteristic  slope  increase. 

The  residual  lift  of  Figures  B1-B6  is  the  .  098  of  Table  4-1.  To  that  has  been 
added  the  Panchenkov  incidence  lift  slope  of  Figure  4-2  to  tick  off  the  zero  pitch  for 
each  incidence  angle.  These  residual  lifts  (zero  pitch  lifts)  agree  reasonable  well 
with  the  data  only  for  the  two  deepest  35-knot  plots.  The  50-lmots  disagreement  is 
considered  to  be  the  pod  cavitation  effect  discussed  in  Section  6. 

These  cross-plots  were  prepared  in  support  of  the  moment  study  of  Section  4 
of  this  report  and  are  considered  further  in  that  section.  These  cross-plots  did  no«^ 
influence  the  lift  considerations  of  Section  4  which  were  based  entirely  upon  the  plots 
presented  in  Reference  4-4. 
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APPENDIX  C 


EVALUATION  OF  FLAP  HINGE  MOMENT  DERIVATIVES 

C  AND  a 
n  n 

01  _ ^ 

GENERAL 

Three  methods  for  predicting  the  three  dimensional  derivatives  from  the  two 
dimensional  derivatives  are  compared: 

1  Perkins  &  Hage,  Reference  4-19,  which  is  the  simplest, 

2  Toll,  Reference  4-3, 

3  DATCOM,  Reference  4-12 

The  AG(EH)  configuration  is  assumed.  Figures  4-15  and  4-16  present  two  new 
sets  of  section  derivitives  for  this  configuration  and  Perkins  and  Hage  would  contribute 
a  third.  For  the  purpose  of  comparing  the  conversions  to  throe  dimensions,  however, 
the  DATCOM  section  derivitives  were  employed: 
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Toll  suggests  model  tests  for  C.  in  Reference  4-3  through  he  was  an  author  for 
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one  of  the  DATCOM  sources  and  the  DATCOM  C,  is  essentially  Toll’s. 
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SUMMARY 
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The  highly  empirical  nature  o  hese  procedures  and  the  scatter  illustrated 
here  generates  little  confidence.  ’  Perkins  &  Hage  procedure  is  adopted  in  this 
volume  because  of  its  simplicity. 
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ANALYSIS  OF  HYSTAD  MODEL  CAVITATION  CHARACTERISTICS 

DTMB  Report  2138.  ’’Model  Test  Results  for  the  Hystad  Strut-Nacelle- Foil  Com¬ 
bination”  by  Peter  Spangler,  January  1966  presents  on  Figure  3  a  "measured  crvitation 
bucket”  at  variance  with  that  predicted  at  Grumman  and  with  that  measured  in  the 
NSRDC  towing  tank  and  in  the  Grumman  WTiirling  Tank  for  the  AG(EH)  planform.  The 
Hystad  and  AG(EH)  planforms  are  identical  so  this  appendix  examines  the  Hystad  results 
to  account  for  the  apparent  discrepancy. 

Note  that  the  Hystad  and  AG(EH)  planforms  are  identical  and  that  the  Hystad  foil 
pitches  relative  to  the  pod  so  that  the  Hystad  lift  is  "incidence  lift, "  comparable  with 
the  NSRDC  AG(EH)  forward  foil  model.  No  incidence  lift  version  of  the  AG(EH)  foil 
was  tested  in  the  Whirling  Tank.  Note,  too,  that  the  Hystad  section  is  16-407  while 
the  AG(EH)  model  sections  are  16-  (.353)08. 

The  data  of  Figure  3  of  the  Spangler  report  are  plotted  as  lift  curves  for  varioi:p 
speeds  on  Figure  Dl.  This  plot  is  analogous  with  Figure  6-4  of  the  main  body  of  this 
volume,  though  Figure  6-4  presents  pitch  lift  chi-iracteristics.  By  reference  to  Figure 
4-2  of  this  volume,  it  may  be  seen  that  the  Hystad  lift  curve  slopes  are  substantially 
in  agreement  with  those  measured  at  NSRDC  and  in  the  Whirling  Tank. 

The  AG(EH)  and  Hystad  lift  curve  slopes  should  be  identical,  but  the  effects  of 
cavitation  on  the  lift  curve  slope  are  not  because  of  the  different  design  lift  coefficients. 

2 

Figure  D2  presents  the  Hystad  drag  curves  plotted  against  C  and  with  the  es- 

L 

timated  separation  drag  removed  to  linearize  the  fully  wetted  portion  of  the  drag  curve. 
The  wetted  slopes  of  Figures  6-5  and  D2  should  be  similar  but  differ  by  25%.  The 
subject  of  induced  drags  lies  outside  the  scope  of  this  report  and  no  attempt  is  made  here 
to  reconcile  the  two  sets  of  drag  curves. 

The  significant  cavitation  characteristics  of  Figures  Dl  and  D2  are  the  initial 
effects  on  the  lift  and  drag  and  the  maximum  lift  and  these  effects  are  plotted  as  speed 
versus  foil  loading  on  Figure  D3  where  they  are  compared  with  the  theorectical  incipient 
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cavitation  bucket  for  the  model.  The  results  are  consistent  with  previous  experience 
with  these  characteristics;  the  upper  surface  leading  edge  effects  leave  the  theoretical 
bucket  in  the  viciaity  of  its  right-hand  corner,  bow  out  to  the  right,  and  tend  to  suck 
back  in  at  lower  speeds.  Lower  surface  cavitation  effect  lies  outside  and  generally 
parallel  with  the  theoretical  bucket. 

The  prototype  incipient  cavitation  bucket  is  also  shown  on  Figure  D3.  Since  the 
effective  boundaries  are  functions  of  lift  coefficient  and  cavitation  number,  they  can  be 
expected  to  shift  tc  the  same  position  relative  to  the  theoretical  prototype  bucket  that 
they  hold  for  the  model  test  condition,  i.e.,  the  significant  scales  on  Figure  D3  are 
cavitation  number  and  C  which  may  be  expressed  as  model  speed  and  model  foil  load- 
ing  or  as  prototype  speed  and  prototype  foil  loading. 

Spangler's  observed  cavitation  bucket  of  his  Figure  3  is  also  shown  on  Figure  D3. 
Note  that  Spangler  has  identified  66  fps  as  first  observed  cavitation  at  4P ,  which  lies 
between  photographs  (a)  and  (b)  of  his  Figure  9  and  was  the  first  speed  at  which  the 
leading  edge  cavitation  on  the  one  semi-span  was  observed.  The  difference  between 
the  two  semi-spans  is  evidently  due  to  the  difference  in  depths  betw^een  the  two  semi¬ 
spans  and  evidently  the  model  is  rolled  to  port  rather  than  starboard  as  the  text  states. 
This  roll  angle  obscures  the  identification  of  first  observed  cavitation  to  some  extent 
but  that  effect  is  not  significant  to  our  discussion  here. 

The  significant  point  is  that  the  starboard  cavitation  of  Figure  9,  like  the 
cavitation  of  both  semi-spans  on  Figure  8.  is  classic  leading  edge  cavitation  and 
corresponds  to  the  right  boundary'  of  the  incipient  cavitation  bucket.  Now,  on  Figure 
10.  note  that  Spangler  has  identified  first  observed  cavitation  with  a  disturbance  at 
the  foil/pod  intersection  which  has  no  relationship  with  the  incipient  cavitation  bucket 
and  which  looks  more  like  the  result  of  a  surface  imperfection. 

The  foil  trailing  edge  disturbances  of  Figure  10  have  no  obvious  explanation, 
though  they  do  have  the  classic  shape  for  cavitation  originating  at  a  point  surface 
imperfection.  Only  if  those  disturbances  are  located  at  80-90%  of  the  semi-span 
could  thcv  be  associated  with  classic  leading  edge  or  mid-chord  cavitation.  Then  the 
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indication  of  Spangler’s  Figures  10-13  is  that  no  classic  leading  edge  or  mid-chord 
cavitation  was  encountered  at  0-3^  just  as  the  theoretical  bucket  would  indicate.  Such 
a  conclusion  cannot  be  certain,  however,  without  extending  the  low  angle  of  attack 
range  to  a  speed  which  would  conclusively  display  the  appearance  of  the  mid-chord 
cavitation. 

If  the  observed  cavitation  boundaries  are  redrawn  from  the  photographs  for  the 
classic  cavitation  pattern,  the  ’’preferred”  observed  cavitation  bucket  of  Figure  D3  is 
obtained  and  it  is  interesting  to  note  that  there  are  cases  where  the  lift  and  drag  are 
affected  before  there  is  visible  evidence  of  classic  cavitation,  A  case  in  point  is  the 
86.4  fps  photograph  of  Spangler’s  Figure  10  for  which  the  li'c  curve  slope  is  significantly 
reduced  though  there  is  no  visible  leading  edge  cavitation  ev on  in  the  next  photograph. 

Note  on  Spangler’s  Figures  14-16  that  the  lower  surface  cavitation  is  all  of 
classic  inboard,  leading  edge  type. 

Finally,  the  theoretical  incipient  bucket  for  the  AG  (EH)  prototype  10  ft.  depth 
is  compared  with  that  for  Spangler’s  model  test  depth  to  illustrate  the  substantial  differ¬ 
ence  in  top- of- the -bucket  speeds  and  foil  loading  for  the  two  cases.  This  difference  is 
particularly  significant  to  the  tow  tank;  the  Whirling  Tank  speeds  may  be  greater  or 
less  than  prototype  speed  and  never  differ  from  prototype  cavitation  speed  by  any  sub¬ 
stantial  amount.  Discussions  of  cavitation  speeds,  then,  must  differentiate  model  and 
prototype  speeds  carefully. 
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SUMMARY 


The  measurement  of  unsteady  load  responses  requires  light,  friction-free  iriech- 
anisms  throughout,  while  the  measurement  of  cavitation  characteristics  require  sturdy 
structure  and  positive  sealing  against  air  leakage.  The  attempt  of  this  program  to 
satisfy  both  requirements  with  the  same  model  has  demonstrated  the  impracticality  of 
such  an  approach.  The  air  sealing  provisions  were  not  adequate  but  those  seals,  with 
the  hinge  required  for  cavitation  scaled  loads,  presented  excessive  friction  during  the 
measurement  of  unsteady  load  responses. 

The  measured  cavitation  characteristics  indicate  a  substantial  effective  "wash-in” 
outboard  on  the  foil  and  are  consistent  with  a  root  wash-out  by  air  leakage.  The  foil 
pitching  moments  contradicted  previous  results  with  the  same  model  throughout  and 
were  consistent  with  the  presence  of  abnormally  high  pressure  (i.e. ,  air  leakage)  aft 
of  the  foil  hinge.  These  indications  persisted  down  to  Froude  scaled  speeds  and  pre¬ 
clude  confident  acceptance  of  any  of  the  data  obtained. 

Analysis  of  the  data  obtained  from  this  program  indicate  areas,  principally  in 
data  reduction,  in  which  improved  procedures  are  required  to  improve  the  effectiveness 
of  future  programs  of  this  type. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  most  obvious  observations  from  this  test  program  are  the  remarkable  sensi¬ 
tivity  of  the  AG(EH)  configuration  to  air  leakage  from  the  strut  and  the  difficulties  intro¬ 
duced  by  friction  into  the  interpretation  of  unsteady  load  measurements.  Both  observa¬ 
tions  indicate  that  cavitatioi  scaled  and  Froude-scaled  testing  should  proceed  with  sep¬ 
arate  models  specifically  designed  to  the  requirements  peculiar  to  these  two  objectives. 
Analysis  of  the  cavitation  scaled  data  obtained  during  this  program  further  indicates  that 
the  interpretation  of  such  data  and  the  probability  of  displaying  all  the  characteristics  of 
interest  would  be  much  improved  by  substantially  expanding  the  test  matrix.  It  is  there¬ 
fore  recommended  that  the  mapping  of  civitation  characteristics  be  conducted  in  contin¬ 
uous  flow  facilities  capable  of  providing  10-100  times  the  high-speed  data  available  to 
this  program. 

Consistent  with  the  observation  that  it  is  impractical  to  measure  cavitation  char¬ 
acteristics  and  unsteady  load  responses  with  the  same  model,  the  current  plan  is  to 
modify  themodel  of  this  test  program  for  the  measurement  of  unsteady  load  responses 
at  Froude-scaled  speeds.  To  aid  the  conduct  of  that  test  program,  the  recommendations 
below  summarize  the  experience  gained  with  the  past  program  in  some  detail. 

1.  The  indications  are  that  the  air  leakage  problem  persists  down  to  the 
Froude-scaled  speeds  and  with  significant  effect.  Further  improvement 
in  the  provisions  for  sealing  the  strut  are  therefore  required. 

2.  There  can  be  no  assurance  that  the  air  calibration  presents  the  submerged 
friction  forces  and  all  of  the  friction  forces  confuse  interpretation  of  the 
air  calibration  and  distort  the  flap  motion.  The  planned  isolation  of  the 
flap  moment  strain  gages  from  the  flap  hinge  will  aid  this  problem  but 
every  advantage  possible  must  be  taken  of  the  reduced  load  requirement  to 
reduce  the  friction  forces  throughout  the  system. 

3.  There  was  confusion  about  the  phase  relationships  throughout  this  program. 
Some  of  this  confusion  resulted  from  an  internal  misunderstanding  at  Grumman 
which  identified  positive  phase  with  lag  and  that  identification  appears  through- 
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out  Reference  El.  Further  confusion  resulted  from  data  transmission 
through  illdefined  procedures  and,  possibly,  from  two  questionable 
practices  which  have  appeared  in  the  general  literature: 

a.  Spv;cifying  the  phase  associated  with  a  negative  amplitude 

b.  Relating  the  phase  to  a  derived  quantity;  e.g. ,  a  rate  or 
acceleration. 

It  is  recommended  that  all  responses  to  flap  motion  be  related  in  phase 
to  the  flap  deflection  and  that  all  responses  to  waves  be  related  in 
phase  to  the  instantaneous  surface  height  (positive  above  MWL)  at 
any  convenient,  but  specified,  point  on  the  foil  chord.  It  is  under¬ 
stood  that  amplitudes  are  half-peak-to-peak,  absolute  (plus)  quantities 
and  that  plus  phases  are  leading. 

4.  Data  reduction,  transmittal,  and  analysis  time  must  be  reduced.  Each 
run  of  the  proposed  program  represents  several  objectives  and  there  is 
considerable  risk  of  losing  a  run  which  has  not  been  analyzed  within 
twenty -four  hours. 

5.  The  Grumman  pre-test  report  for  the  new  test  program  must  define  the 
transmittal  of  data  in  an  explicit  and  detailed  manner.  The  data  require¬ 
ments  of  the  Aero  Test  and  Marine  groups  must  be  distinguished  as  to 
content  and  timing;  a  data  identification  system  must  be  established;  data 
transmittal  duplications,  as  for  data  reduction  corrections,  must  be 
explained. 

6.  The  air  calibration  smoothing  procedure  of  Reference  El  was  not  satis¬ 
factory  because  it  reproduces  the  data  scatter  and  because  it  does  not 
distinguish  the  friction  and  inertial  sources  in  the  calibration.  Each  of 
the  responses  for  these  calibrations  should  be  considered  individually, 
their  inertial  component  estimated,  the  results  plotted  in  appropriate 
manner  and  smoothed  by  eye.  For  this  purpose  the  inertias  of  all  moving 
parts  should  be  estimated  and  documented  in  the  pre-test  report.  The  air 
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calibrations  are  the  first  and  best  test  the  data  reduction  for 
phase  relationships. 

7.  The  programmed  static  (submerged  zero  speed)  runs  are  intended  to 
test  Theodorsen^s  inertial  coefficients  but  the  validity  of  this  test  is 
not  yet  certain;  this  data  is  academic  rather  than  practical. 

8.  If  the  phase  relationships  are  well  defined  and  the  air  calibrations 
are  reasonable,  the  analysis  of  the  response  data  is  routine. 

9.  Some  consideration,  not  yet  defined,  should  be  given  to  the  unsteady 
responses  at  harmonics  of  the  flap  motion  to  test  the  flap  motion 
waveform  achieved,  test  the  superposition  of  the  responses,  and 
extend  the  test  frequency  range. 
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DISCUSSION 


Steady  State  Data  -  General 

The  reduced  steady  state  data  available  at  this  writing  is  presented  graphically 
on  Figures  El  -  E20  and  the  corresponding  test  conditions  are  tabulated  in  Table  El. 

All  reduced  data  is  from  Reference  El  except  the  zero  pitch  cavitation  data  of  Figures 
E5  -  E9.  That  data  is  incorrect  in  Reference  El  and  hand  reduced  data  recorded  on 
the  carriage  is  employed  here. 

It  must  be  recognized  that  this  data  presents  fragments  of  a  very  complex  and 
unfamiliar  set  of  characteristics;  Figure  E21  is  from  data  obtained  from  the  Grumman 
Whirling  Tank  for  a  high  aspect  ratio  flapped  foil  model  and  the  test  conditions  were 
extended  well  into  the  upper  and  lower  surface  cavitation  ranges  to  aid  interpretation. 
This  data  is  cross-plotted  from  fixed  flap  angle  pitch  curves  which  have  not  yet  been 
analyzed  and  may  change  in  final  form.  Note  that  the  combination  of  high  aspect  ratio 
and  double  struts  has  produced  almost  2-dimensional  litt  characteristics  for  this  model. 

Note  that  no  portion  of  Figure  E21  is  predictable  at  present.  There  is  a  predict¬ 
able  slope  between  2°  and  6°  flap  angle  at  the  35-knot  speed.  At  the  lower  speeds  there 
is  a  characteristic  but  unpredictable  slope  increase  and  at  higher  speeds  a  character¬ 
istic  slope  leduction  due  to  cavitation.  The  zero  lift  angle  for  the  normal  slope  has 
been  radically  altered  by  some  still  unknown  effect  so  that  even  this  portion  of  the  lift 
curve  is  unpredictable.  The  separated  flow  slope  for  positive  flap  angles  at  low  speed 
is  normal;  the  supercavitated  slopes  are  higher  than  normal.  The  positive  flap  sep¬ 
arated  and  supercavitated  slopes  are  repeated  for  negative  flap  at  25  and  55  knots;  at 
intermediate  speeds  the  negative  flap  slopes  are  still  in  transition.  A  fragmentary 
view  of  these  characteristics  would  pose  a  substantial  data  analysis  problem. 

The  data  of  Figures  El  -  E20  present  such  a  problem  and  one  which  has  been 
complicated  by  the  fact  that  these  are  abnormal  results.  The  lift,  drag,  and  moments 
are  interrelated  and  are  all  related  to  the  cavitation  characteristics.  Thus  interpreta¬ 
tions  given  the  lift  data  of  Figures  El  -  E20  are  hypothetical  until  they  have  been  shown 
consistent  with  the  other  measured  characteristics  at  their  own  and  at  all  the  other  test 
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conditions.  The  final  results  of  this  Appendix  evolved  from  several  iterations  of  this 
process  but  remain  hypothetical  in  nature  because  of  the  uncertain  character  of  the 
abnormality  involved. 

Figure  E22  was  employed  for  guidance  in  the  interpretation  of  the  lift  results. 

It  was  assumed  that  separation  effect  was  of  rhe  order  of  the  difference  between 
small  and  large  deflections  of  a  plain  flap  wn  le  flap  cavitation  was  of  the  order 
of  the  split  flap  effect  on  Figure  E22. 

Lift  Characteristics 

The  total  lift  coefficient  m  ay  be  expressed  in  the  form*. 
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which,  for  a  trimmed  flap,  becomes; 
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and  from  the  results  of  Section  4,  for  the  AG(EH)  forward  foil  this  becomes  (at  10  ft. 
submergence); 


C  =  2.75  (  o  +  ,838i  )  +  .  090 
^6=0 


(E.3) 


If  various  combinations  of  pitch  and  incidence  are  tested,  Eq.  (E.3)will  test  pitch 
lift  curve  slope,  incidence  lift  curve  slope,  and  residual  lift  simultaneously.  Further, 
the  flap  lift  curve  slope  is  tested  by  the  C  at  zero  flap.  The  lift  curve  results  of 

Lj 

Figure  El  -  E20  are  tested  in  this  manner  on  Figure  E23. 
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The  slopes  upon  which  Figure  E23  are  based  are  tabulated  in  Table  El.  It  was 
not  possible  to  employ  a  common  slope  through  the  lift  data  and  the  slopes  of  Table  El 
and  Figure  E23  are  the  solid  lift  curves  of  Figures  El  -  E20.  The  corresponding  flap 
effectivenesses  are  listed  in  Table  El  and  are  disappointing;  only  the  lowest  speed 
Froude  scaled  tests  presented  a  normal  flap  effectiveness.  Flap  effectiveness  appears 
to  be  subject  to  incidence,  pitch,  and  speed  but  not  in  any  usefully  systematic  manner. 

The  residual  lift  discrepancy  of  Figure  E23  was  noted  early  in  the  test  program 
and  initiated  a  substantial  effort  to  control  air  leakage  over  the  foil.  Elimination  of 
the  air  leakage,  visually  and  under  pressure  test,  increased  the  residual  lift  but  not  to 
the  expected  value.  Investigation  disclosed  that  the  incidence  angles  of  Reference  E2 
were  relative.  The  subject  is  not  specified  in  Reference  E3  though  it  is  standard  opera¬ 
ting  procedure  for  the  Grumman  model  shop  to  provide  a  template  locating  the  foil  plane. 
Nevertheless  these  tests  proceeded  on  what  now  appears  to  have  been  a  mistaken  assump¬ 
tion  that  this  was  the  first  absolute  measure  of  the  residual  lift  for  the  model.  Note  on 
Figure  E23  that  the  residual  lift  is,  in  some  cases,  speed  dependent  though  the  residual 
lift  for  the  Froude  scaled  tests  is  as  low  as  that  for  the  cavitation  scaled  tests.  Only 
the  lowest,  zero  pitch  speed  presented  a  residual  lift  approaching  that  of  the  1962  model 
tests. 

While  the  assumption  of  a  common  slope  and  intercept  does  violence  to  the  meas¬ 
ured  lift  data,  it  is  essential  that  this  be  done  for  performance  and  control  studies. 

The  flap  effectivenesses  for  the  cavitation  data  of  Table  El  seem  grouped  about  values 
of  .37  and  .31,  which  may  indicate  distinct  boundary  layer  conditions.  For  the  purpose 
of  characterizing  these  results  an  effectiveness  of  .34  and  a  residual  lift  of  .04  were 
employed  and  the  resulting  lift  curves  are  shown  as  dotted  lines  on  Figures  E5  -  E20. 

The  corresponding  lift  description  is 

C  =C  (y  +  C  i+C  do  6  +  C 
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The  result  provides  a  rather  good  correUticii  for  the  cavitation  lift  data. 

Foil  Hinge  Moment 

Foil  hinge  moments  were  measured  only  on  the  Froude  scaled  runs  and  do  not 

provide  a  measure  of  and  .  If  the  towing  tank  derivatives  of  Table  4-3  are 

a  i 

assumed,  the  foil  hinge  moments  of  Figures  El  -  E4  can  be  compared  with  each  other 
as  on  Figure  E24.  Only  one  set  of  four  points  on  Figure  E24  display  the  expected  slope 
and  they  are  associated  with  a  profound  change  in  residual  moment.  From  Figures 
E23  and  E24  the  1962  and  1972  tests  can  be  compared  as: 

AC 
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That  is,  the  loss  in  residual  lift  apparently  had  a  centroid  507c  aft  of  the  hinge  on 
the  total  foil  MAC  or  at  88.57  of  the  MAC  for  the  four  points  of  normal  3lope.  For 
most  of  the  points,  however,  that  centroid  is  at  65.57;  of  the  MAC  with  an  indication, 
discussed  below,  that  it  is  located  later? ily  at  the  pod. 


For  the  flap  effectiveness  above  the  C 
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’s  of  Table  4-3  produce  foil  hinge 
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moment  slopes  of: 


"  .0402  X.34  C.,  =  .  01367  C„ 

=  .  01367(-.  0985)  =  -.  C0i347/deg  (theoretical) 

=  .  013G7( -.  0689)  =  -.  000942/dcg  (tow  tank)  (E.6) 

and  these  slopes  are  shown  on  Figure  f^24  with  their  associated  residual  hinge  moments 
of  Table  4-3. 

On  some  ’’effective  mean  aerodynamic  chord”  for  pitch  moment  Ihe  flap  load 
distribution  is  described  by: 
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where  G*  denotes  the  absolute  value  of  Allen’s  G  to  avoid  associating  a  sign  with  distance. 
A  problem  in  symbols  arises  here  because 
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MAC  „  is  the  chord  station  for  the  hinge  on  the  effective  MAC,  given  by 
eii 

The  ratio,  MAC /MAC  „,  is  given  by 

eii 

MAC  ^  1  -  ^MAC  (1  -  X) 
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.  332-1.  304(-.  0985) .322  +  .1285  =  .^50  (theoretical) 

=  .  322-1. 304(-.  0689)  =  .322  +  .09  =  .412  (1962  towing  tank 

result  and  four  points 
of  1972  result) 

.  322-1. 304(-.  0032/.  0137)  =  .  322-1. 304(-.  2335)  =  .322  +  .301 
=  .626  (1972  towing  tank  results) 

There  are  too  many  uncertainties  invclved  here  for  confident,  quantitative  inter¬ 
pretation  but  the  qualitative  indication  that  the  effective  MAC  for  pitch  moment  has  been 
moved  outbuard  some  20%  of  the  semi-span  is  ve-'y  significant-  The  effect  of  the  7% 

shift  of  MAC  shift  of  xMAC  on  Fig^ire  1-2  on  uie  fc  /C  )  max  is  shown  on  Figure  6-1. 

i  b 

The  effect  her^^  is  one  which  "washes-out"  some  portion  of  the  foil  at  the  root. 

If  the  upper  four  foil  hinge  moment  points  of  Figure  E4  really  present  a  normal  condition, 
rather  than  scatter,  the  lift  curve  of  that  fig:ure  indicaies  that  the  lift  curve  slope  is  not 
affected.  Thus  Figure  indicates  that  this  "wash  out”  is  independent  of  flap  aiigle  and 
all  of  the  lift  curves  demonstrate  that  the  effect  i«  practically  independent  of  speed. 
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Therefore  the  effect  is  one  of  foil  twist,  in  a  direction  to  wash  out  the  center  section. 

In  summary,  these  results  present  two  possible  foil  hinge  moment  curves.* 

=  .00384(y°  +  .002815io  -.000942^0  -.01 

.00384q,°  +  .002815i°‘  -.00326°  -.0215 

wnere  neither  residual  moment  and  only  the  lower  C„  are  in  agreement  with  the 

1962  model  test  results.  Both  of  these  curves  are  shown  on  Figures  El  -  E4;  the  first 
as  a  solid  curve  and  the  second  as  a  dashed  curve. 

Flap  Hinge  Moment 

The  flap  hinge  moment  data  for  common  pitch  and  incidence  angles  is  presented 
on  Figures  E27-E33.  On  these  figures  the  moment  slopes  and  trimmed  flap  hinge 
moments  as  judged  by  eye  were  plotted  against  j)itch  angle  on  Figures  E25  and  E26. 

No  quantitative  effort  is  made  to  judge  these  results  because  the  lift  disturbance  is 
centered  in  the  region  of  the  flap  and  must  have  had  a  direct  effect  on  the  residual  flap 
hinge  moment  by  unloading  the  flap.  A  further  effect  on  the  residual  flap  hinge  moment 
and  upon  the  hinge  moment  slope  would  result  from  the  apparent  redistribution  of  the 
spanwise  loading  associated  with  the  lift  distu’^bance.  There  are  significant  qualitative 
indications  to  be  considered  in  these  results  however. 

It  is  to  be  noted  that  the  steep  slopes  of  Figure  E25  are  associated  with  the  extreme 
angles  of  attack;  i.e. ,  the  effect  is  opposite  that  to  be  expected  as  a  result  of  flow  separ¬ 
ation.  The  view  is  clarified  only  a  little  by  presenting  the  hinge  moments  plotted  directly 
against  flap  angle  as  on  Figure  E34.  Here  it  was  necessary  to  present  faired  results; 
the  actual  curves  are  much  more  complex  because  an  interruption  in  the  lift  or  hinge 

moment  curve  porduces  an  effect  in  the  C  -6  curve, 

L 

Figure  E34  presents  either  a  non-symmetric  slope;  i.e.,  different  slopes  for 
positive  and  negative  displacements,  or  a  symmetric  discontinuity  about  a  shallow  slope 
in  the  region  of  the  test  data.  Interpretation  is  hampered  by  the  lift  disturbance  known 
to  be  in  the  data  and  by  the  foil  camber.  In  either  event,  the  slope  discontinuities  of 
Figure  E34  are  of  considerable  significance  to  the  control  design. 


These  slope  discontinuities  are  modified  by  the  iii\  disturbance  but  they  are  not 
introduced  by  it.  Similar  discontinuities  are  evident  in  Figures  48-51  of  Reference  £6 
which  presents  data  for  a  16-303  section.  It  is  to  be  noted  that  the  camber  for  the 
AG(EH)  model  is  the  equivalent  of  20. 7  degrees  of  fliq)  which  introduces  a  substantial 
assymetry  into  whatever  boundary  layer  effects  are  at  work.  Note  that  most  of  the 
literature  on  plain  flap  hinge  moments  is  for  symmetric  sections.  It  may  be  necessary 
to  obtain  comparable  hydrodynamic  data  before  attempting  to  quantitatively  analyze 
hydrodynamic  data  for  cambered  sections. 

This  hinge  moment  provides  no  useful  quantitative  result  for  prototype  application 

and  the  flap  hinge  moment  prediction  of  Section  4  will  continue  to  be  employed  except 

that  it  does  now  seem  to  be  advisable  to  reduce  the  residual  hinge  moment  to  C.  /6 

i  eff 

rather  than  the  C.  /2  of  Section  4.  This  prediction  is  compared  with  the  faired 
i  eff 

measured  hinge  moments  on  Figures  E1-E20. 

Cavitation 


For  zero  pitch  Eq.  (6.2.19)  may  be  written 
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where  C  is  proportional  to  the  local  lift  coefficient  increment  at  75%  of  the  semi-span 
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due  to  thp  lift  disturbance. 

Eq.  (E.14)  may  be  written 
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where  C  is  the  foil  lift  coefficient  for  incipient  cavitation. 

Li 
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For  the  upper  suri  ^?ce  leading  adge  Eq.  (E15)  becomes 
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The  incipient  cavitation  coefficients  of  Figures  E5  thru  E9  are  plotted  in  this 

manner  on  Figure  E35.  Assuming  the  validity  of  the  cavitation  prediction,  Figure  E35 

indicates  that  the  lift  disturbance  is  not  a  stable,  or  consistent^  condition  though  it  is 

assumed  so  in  the  calculations  to  follow.  There  is  an  indication  that  the  lift  disturbance 

does  not  affect  the  spanwise  distribution  associated  with  angle  of  attack  since  it  could 

only  be  expected  to  increase  c  /C  and  this  would  be  inconsistent  with  the  results  on 

I  L 

Figure  E35. 


From  Figure  E35  it  is  assumed  that  the  C  ref  is  effectively  increased  by  .118 

L 

at  75?o  of  the  semi-span  which  is  equivalent  to  about  2-1/2  degrees  of  ’’wash-in"  at 

this  point.  Only  one  test  condition,  shown  on  Figure  E15,  presented  a  measure  of  the 

inboard  effective  twist.  A  direct  solution  of  Eq.  (E14)  for  this  condition  indicated  a 

C  of  00344  which  is  quite  low  and  indicates  either  that  the  lift  disturbance  is  wide- 
T 

spread  or  that  it  is  concentrated  on  the  pod  at  the  strut  intersection. 

Figure  E36  compares  the  zero  pitch  cavitation  characteristics  with  the  expected 
cavitation  bucket  to  indicate  the  severity  of  the  problem  encountered  in  this  data.  The 
incidence  control  bucket  is  compared  with  the  flap  bucket  on  this  figure  and  the  upper 
surface  leading  edge  boundaries  are  identical  for  the  two  buckets  when  the  flap  is  trim¬ 
med.  The  incidence  control  bucket  has  been  confirmed  by  the  1962  tests  so  at  the  test 
speed  of  Figure  E36  for  a  45.7  knot  prototype  speed  there  is  a  direct  comparison  of  the 
1962  and  1972  model  cavitation  characteristics.  The  low  speed  case  of  Figure  E36 
presents  a  fertile  field  for  speculation.  If  this  is  a  normal  case  flaps  will  present 
rather  poor  drag  characteristics  and,  from  Figure  E29,  a  monumental  non-linearity 
problem.  It  is  more  likely,  however,  that  at  1150  prototype  W/S  on  Figure  E36  air 
displaced  a  separated  boundary  layer  over  some  portion  of  the  flap  span  and  produced 
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the  drag  and  flap  hinge  moment  disturbances. 

Figures  E37-E41  compare  the  incipient  cavitation  bucket  for  the  assumed  foil 
twist  with  the  observed  cavitation  characteristics  at  each  of  the  pitch  angles  tested  and 
Indicate  a  reasonable  consistency  for  the  assumed  nature  of  the  difficulty.  Figure  E42 
presents  a  typical  drag  result.  Tr.**  one  out-of-line  hand  reduced  point  is  clearly  ident¬ 
ified  on  Figure  E9  as  ventilation  and  demonstrates  that,  at  least  intermittently,  there 
was  air  in  the  strut  available  to  the  foil.  Figure  E43  presents  one  of  two  cases  where 
the  minimum  drag  was  affected  without  effect  on  the  induced  drags  and  might  present 
a  case  where  an  entire  run  was  affected  as  was  one  point  on  Figure  E42. 

Figures  E37-E41  indicate  remarkable  penetration  of  the  leading  edge  cavitation 
boundary  with  a  minimal  lift  and  drag  disturbance.  This  characteristic  of  the  effectively 
twisted  AG(EH)  foil  is  similar  to  the  transit  foil  and,  in  fact,  physical  twist  would  be 
another  way  to  achieve  transit  foil  characteristics.  The  cavitation  objective  for  this 
test  program,  however,  was  to  demonstrate  the  cavitation  characteristics  when  the 
chordwise  and  spanwise  load  distributions  were  more  uniform  than  was  available  for 
these  results.  Figure  E39  presents  the  only  anticipated  case  which  was  achieved  and 
then  not  at  the  anticipated  pitch  angle. 

In  summary,  the  cavitation  test  results  do  not  achieve  the  test  objective  but  in¬ 
stead  demonstrate  a  remarkable  sensitivity  of  those  characteristics  to  foil  root  ventil¬ 
ation. 

Dynamic  Hinge  Friction  Characteristics 

With  reference  to  the  flap  angle  as  a  sine  wave,  the  sliding  hinge  friction  is 
assumed  to  have  tlie  shape  of  a  square  wave: 


Where  only  one  term  of  the  series  defining  the  friction  is  employed,  as  in  these 
results,  the  sliding  friction  is  a  cosine  wave  of  amplitude: 

2  .  fn/2  r  Tt 

a^  =  “  [  1  (-f  )  cos  X  d  X  +  I  f  cos  x  d  x  J 
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so  that  this  idealized  sliding  friction,  as  the  dynamic  data  is  reduced,  lags  the  flap 
deflection  by  90° . 


For  the  same  reference,  the  breakaway  friction  is  assumed  to  have  the  shape  of 


a  pulse: 


The  breakaway  friction,  then,  can  be  defined  as  a  sine  function  having  a  first  term 


of  amplitude: 
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which  is  in  phase  with  the  flap  deflection.  The  inertial  loads  are  in  phase  with  the  flap 
deflection  so  that  the  moments  measured  in  the  air  calibration  have  the  relationship. 
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The  flap  hinge  moments  measured  during  the  air  calibration  are  listed  in  Table 
E2  and  it  was  necessary  to  invert  the  phases  supplied  by  NSRDC  to  obtain  an  interpre¬ 
tation  with  physical  significance. 

The  real  part  of  the  measured  hinge  moment  while  cycling  the  flaps  in  air  is 
shown  on  Figure  E44  which  displays  tY  i  inertial  effect.  Removing  the  inertial  effect 
from  the  real  part  of  the  measured  hinge  moment  produces  the  measured  friction 
moments  of  Figure  E45  which  indicates  that  the  'Treeing  up"  process  on  the  flap  hinge 
actually  increased  the  hinge  friction.  Figure  E45  indicates  friction  moment  amplitudes 
of  2.47  in.  lbs.  before  freeing  up  the  hinge  and  3.97  in.  lbs.  after  freeing  up  the  hinge. 
These  are  interpretations,  of  course,  which  place  a  considerabie  strain  upon  the  data 
precision  for  these  measurements. 

Figures  E44  and  E45  indicate  hinge  moments  in  air  of 

2 

h  -  ,  000485UJ  ^  -.l-i.l8  before  free-up 

a 

2 

-  .  000485fi)  6  -.l-i.315  after  free-up  (E.19) 

Figure  E4G  compares  Eq.  (E.19)  with  the  smoothed  calibration  of  Table  II  of 
Reference  El.  The  smoothing  process  of  Reference  El  produces  a  better  set  of  results 
for  a  given  set  of  points  but  with  no  relation  to  the  physical  mechanisms  involved  and 
with  no  flexibility  for  new  deflections  and  frequencies. 


Air  Calibration 


Eq.  (E.19)  constitutes  the  air  calibration  assumed  in  this  note  for  the  flap  hinge 
moment. 

The  entire  flap  control  mechanism  rests  upon  the  lift  gage  blocks  so  the  friction 
forces  do  not  appear  in  the  measured  lift  which  should  present  an  amplitude  proportional 
to  and  180^  out  of  phase  with  the  flap  displacement.  The  measured  amplitudes  are 
shown  on  Figure  E47  and  the  change  in  slope  after  freeing  up  the  hinge  has  not  been 
explained.  The  inertia  producing  the  lift  force  of  Figure  E47  is  due  to  all  the  moving 
parts  of  the  flap  control  system  of  Figure  2  of  Reference  El.  The  inertia  has  not  been 
estimated  though  it  is  obviously  substantially  greater  than  that  of  Figure  E44  which 
presents  only  the  inertia  below  the  flap  moment  dynamometer.  The  weights  of  Figure 
E47  assume  all  of  the  moving  mass  to  be  concentrated  on  a  one  inch  arm  for  qualitative 
comparison  with  the  dotted  line  which  presents  the  estimated  three  pound  weight  below 
the  flap  dynamometer.  No  masses  were  changed  when  the  flap  was  freed  up  and  the 
indication  is  that  one  or  the  other  of  the  two  sets  of  calibrations  is  in  error. 


The  measured  lift  phases  are  presented  on  Figure  E48  which  does  not  provide  any 
reliable  variation  of  phase  with  frequency.  Note  that  it  was  necessary  to  invert  the 
reported  phase  to  give  physical  meaning  to  the  result. 


From  Figures  E47  and  E48  the  calibration  assumed  was 
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.  0224m>  6  161°  before  freeing-up  hinge 
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.  05  0)  6  161°  after  freeing-up  hinge 


(E.20) 


which  is  compared  with  the  calibration  after  freeing  the  hinge  and  with  the  smoothed 
calibration  of  Table  II  of  Reference  El  on  Figure  E49.  In  this  instance  the  machine 
smoothed  calibration  was  led  into  significant  error  by  the  inclusion  of  wild  points. 


As  with  lift,  the  friction  forces  do  not  register  on  the  drag  gages  but  the  data 
reduction  procedure  eliminates  most  of  the  drag  response  to  flap  motion.  The  analysis 
here  proceeded  only  far  enough  to  indicate  a  form  for  the  air  calibration. 
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Only  the  flap  mass  is  considered  for  illustration  though  Figure  2  of  Reference  El 
indicates  that  other  masses  are  in  motion,  some  countering  the  effect  of  the  flap  mass. 
Assuming  the  flap  mass  concentrated  at  mid-flap-chord,  that  mass  is  located  on  x  at 


so  (E.21)  may  be  written 


^  =  i  ^f  ■  2 


where  ^  =  6  ^  sin  o)  t 


-  c^  cos  (6^  sin  t) 


The  acceleration  on  x  may  be  written 
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One  extreme  is  given  at  t  =  0: 
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The  other  extreme  is  given  at  sin  u)t  =  +1,  cos  (nt  =  0: 
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If  cos  is  approximated  by  unity  and  sin  6^  is  approximated  by  6^.  the  acceleration 


is  bounded  by: 
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Eq.  (E:23)  is  quite  a  complex  function  having  a  fundamental  frequency  double  that 

2 

of  the  flap  motion,  indicated  by  the  cos  term,  and  a  distorted  sinusoidal  wave  form. 

2 

From  the  form  of  Eq.  (E.  26)  the  air  calibration  was  first  examined  in  terms  of  D/6 

on  Figure  E50  but  the  systematic  variation  with  flap  amplitude  indicates  that  this  is  an 

unwise  procedure.  Evidently  the  amplitude  at  the  fundamental  frequency  is  a  function 
2 

of  5  while  the  amplitude  at  the  flap  frequency,  provided  by  the  data  reduction  procedure, 
is  more  a  function  of  6 . 

At  the  flap  frequency,  the  first  harmonic  of  the  basic  inertial  frequency  for  drag, 
only  the  distortions  in  the  fundamental  wave  form  are  viewed  and  measurement  precision 
can  be  expected  to  suffer  in  amplitude  and,  especially,  in  phase.  Therefore  the  calibra¬ 
tion  is  defined  here  purely  in  empericai  terms.  The  amplitude  is  presented  on  Figure 

'  2 
E51  in  terms  of  D/6  and  has  arbitrarily  been  fit  to  a  quadratic  in  m  for  convenience. 

The  phase  is  presented  on  Figure  E52  where  the  low  frequency  phases  are  questionable 

because  of  the  low  amplitudes  with  which  they  are  associated.  Two  forms  of  curve  fit 

were  examined  and  that  which  neglects  the  lowest  frequency  points  is  preferred  because 

it  appears  to  provide  better  correlation  for  those  points  of  significant  amplitude. 

From  Figures  E51  and  E52  the  drag  calibration  in  air  which  was  adopted  was 
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The  calibrations  of  this  appendix  and  of  Reference  E.l  are  compared  with  the  measured 
points  on  Figure  E53.  There  is  no  objective  measure  of  the  relative  validity  of  the  two 
smoothed  phase  calibrations. 


Note  that  time  does  not  permit  the  very  elaborate  analysis  required  to  provide 
physical  significance  to  the  drag  phase  at  flap  frequency:  the  measured  phases  have 
been  reversed  here  only  bec.iuse  that  reversal  was  physically  required  for  the  other 
three  air  calibrations. 

Time  limitations  also  prevent  analysis  of  the  physical  significance  of  the  foil  hinge 
moment  calibration.  Normally  the  foil  hinge  moment  would  be  related  to  a  flap  load 
and  moment  by 


H  =  h  -c(““-)L 
a  a  c  c 


(E.28) 


but  in  this  case  there  are  masses  and  frictions  in  the  fcil  hinge  moment  case  which  are 
not  in  the  flap  moment  calibration  and  the  analysis  of  the  air  calibration  is  complicated 
by  the  model  linkage  arrangement  provided  to  transmit  hydrodynamic  flap  load  and 
moment  to  the  foil  pivot.  Eq.  (E.28)  has  been  employed,  then  only  to  provide  a  logical 
form  in  which  to  display  the  foil  moment  calibration. 


Related  to  the  exposed  foil  mean  aerodynamic  chord,  for  convenience,  Eq.  (E.28) 


becomes 


H  =  h  -  .348  L 
a  a  a 


(E.29) 


and  from  Eq.  (E.19)  and  (E.20)  the  numerical  evaluation  is 

H  =  [.01694  6  -.1  ]  +  i  [ -.  00566  6  - 


where  .18  -  before  freeing- up  hinge 


(E.30) 


.315  =  after  freeing-up  hinge 

The  real  part  of  the  air  calibration  is  shown  on  Figure  E54  where  the  scatter 
practically  precludes  interpretation.  From  the  results  after  freeing-up  the  hinge,  the 
slope  of  Eq.  (E.  30)  is  employed  for  Figure  E54  and  curves  are  drawn  for  two  of  the 
calibrations.  The  scatter  is  still  worse  for  the  imaginary  terms  on  Figure  E55  and  one 
curve  is  drawn  to  represent  two  of  the  calibrations. 
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The  solid  lines  of  Figures  E54  and  E55  represent  the  calibrations: 


H  =  C.  00224(1)^  6  -  ]  -  i  [.00568  u)^  6  +  .5] 

a  1.25 


where  1.84  =  before  fresing-up  hinge 

1.25  =  after  freeing-up  hinge  (E.31) 

which  are  compared  with  their  respective  measured  calibrations  on  Figures  E56  and 
E57.  The  result  is  most  unsatisfactory  on  Figure  E56  but  that  measured  case  appears 
to  be  in  conflict  with  the  flap  moment  and  lift  calibrations,  individually,  and  with  the 
foil  moment  case  of  Figure  E.57. 

It  should  be  noted  that  calibrations  over  frequency  and  amplitude  ranges  simul¬ 
taneously  hamper  considerations  of  the  result.  A  better  procedure  would  be  a  range 
over  frequency  at  two  or  three  amplitudes  followed  by  a  range  over  amplitude  at  two 
or  three  frequencies. 

In  summary,  the  air  calibrations  assumed  for  the  freed-up  hinge  are 
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The  measured  flap  hinge  moments  in  air  are  tabulated  in  Table  E2.  The  measured 
foil  hinge  moment,  lift,  and  drag  in  air  are  presented  in  Table  E3. 


static  Calibrations 


No  zero  speed,  submerged  runs  were  made  after  the  hinge  was  freed  up  and  the 
results  of  only  one  of  the  three  static  calibrations  made  was  available  to  this  study  with 
the  results  tabulated  in  Table  E4. 

For  a  zero  speed,  Theodorsen’s  equation  for  flap  hinge  moment  reduces  to: 

,  tt  „  3^  2. 

h  =  -  7  oSc  T-Cps  6 
4  3 

^  ■  17^; 


where  p  =  1.9384  (assumed) 

S  =  1.27  -  exposed  foil  area 

c  =  .  719  =  exposed  MGC 

T  =  -  ,  007n  from  Table  3-1 
3 


For  these  values  Eq.  (E.33)  reduces  to 


h  =  .000396(1)  6 


(E.33) 


(E.34) 


This  equation  is  compared  with  the  raw  measured  data  and  with  that  data  corrected 
for  the  friction  of  Table  II  of  Reference  E.l  and  of  Eq.  (E.19)  of  this  appendix  on  Figure 
E58.  The  applicability  of  Theodorsen’s  equation  to  the  3-dimensional  case  is  question¬ 
able,  particularly  at  finite  sumbergence,  so  the  scatter  of  Figure  E58  precludes  con¬ 
clusive  judgement.  Figure  E38  might  indicate  that  Theodorsen’s  T  coefficient  does 

O 

not  apply  to  this  case  or  that  the  hinge  friction  was  increased  some  17%  in  water.  No 
significant  difference  between  the  two  friction  calibrations  is  indicated  but  the  result 
does  confirm  that  the  Phases  supplied  are  inverted. 

For  a  zero  speed,  Theodorsen's  equation  for  foil  hinge  moment  reduces  to: 


E-22 


where  P,  S,  and  c  are  as  for  £q.  (£.33) 


=  .00428TT  from  Table  3.1 
=  -.02322tt  from  Table  3.1 


(£.35) 


For  these  values  £q.  (£.35)  reduces  to 


H  =  .00523(„  6 


(£.36) 


which  is  compared  with  the  raw  data  and  with  two  sets  of  corrected  data  on  Figure  £59. 
As  for  flap  hinge  moment,  the  applicability  of  Theodorsen’s  coefficients  in  3-dimensions 
at  finite  submergence  is  in  question  but  Figure  £59  is  inconclusive  because  of  data 
scatter.  The  Reference  £1  calibration  appears  to  be  preferable  for  amplitude  but 
neither  air  calibration  is  adequate  for  phase  unless  there  is  a  substantial  free  surface 
effect. 


For  a  zero  speed,  Theodorsen’s  equation  for  foil  lift  reduces  to; 

T  1  oK 

L  =  “  pSc  T  cu  0 
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which,  for  the  values  of  £qs.  (£.33)  and  (£.35),  reduces  to 
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=  .01161^  6  /180° 


(£.38) 


£q.  (£.  38)  is  compared  with  the  raw  data  and  with  two  sets  of  corrected  data  on 
Figure  £.  60.  Little  choice  between  the  air  calibrations  is  i.idicated  but  there  is  an 
indication  here  that  the  Theodorsen  coefficient  does  not  describe  the  case  tested.  Note 
that  for  this  case  it  was  necessary  to  employ  the  phase  as  reported  to  derive  a  cred¬ 
itable  result. 

Unsteady  Load  Responses  -  General 

The  steady-state  characteristics  considered  in  the  cavitation  study  were  not  of 
adequate  confidence  level  to  provide  a  basis  for  ll.v,  ...c ^surement  of  unsteady  load  re¬ 
sponse,  Those  responses  are  considered  here  L'^ly  to  provide  guid.v'ce  for  a  new  un¬ 
steady  load  test  series.  Two  subjects  are  of  particular  interest  in  this  connection: 
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The  reliability  of  the  data  reduction  and  transmittal  procedures, 
with  particular  reference  to  the  phase  angles, 

The  significance  of  t?ie  air  calibration,  including  friction  effects, 
to  the  measured  hydrodynamic  responses. 


Flap  Lift  Response 

The  lift  responses  transmitted  from  NSRDC  through  Grumman  Aero  Test  are 
presented  in  Table  E5.  The  phases  have  been  inverted  and  the  sign  reversed.  Revers¬ 
ing  the  sign  accounts  for  the  difference  in  sign  convention  between  Reference  E.l  and 
this  report  but  the  source  for  the  phase  inversion  has  not  been  identified. 

The  data  supplied  has  had  the  air  calibration  removed  in  accordance  with  the 
relationship: 
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where,  from  Eq.  (E.32) 
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where,  from  Table  I  of  Reference  El,  C  is  .0198/deg.  for  the  low  speed  run  and 
.0167/deg.  for  the  high  speed  run. 

The  data  with  and  without  the  air  calibration  removed  are  compared  with  the  final 
result  of  Reference  El  on  Figures  E61  and  E.62  which  should  also  be  compared  with 
Figure  22  of  Reference  El.  The  solid  line  of  Figures  E61  and  E62  attempts  to  represent 
both  test  speeds  after  freeing  up  the  hinge  and  does  so  on  a  gain/frequency  plot  which 
exaggerates  the  scatter  at  high  frequency.  The  attenuation/wavelength  presentation  is 
preferred  here  because  it  identifies  that  scatter  with  the  unbounded  inertial  effect  at 
wavelengths  of  questionable  significance.  Figures  E61  and  E62  do,  nevertheless,  idi« 
cate  that  the  air  calibrations  of  Reference  El  and  this  appendix  do  differ  substantially 
in  their  effect  on  the  measured  amplitudes  though  not  on  phase. 

Figures  E63  and  E64  compare  the  flap  lift  hydrodynamic  responses  (measured 
less  air  calibration)  for  four  runs  with  the  result  of  Reference  El  and  with  the  most 
appropriate  predictions  of  Section  5  of  this  report.  Figures  E63  and  E64  should  be 
examined  in  connection  with  Figures  5-20  and  5-21.  The  persistent  attenuation  to  long 
wavelengths  may  indicate  the  presence  of  air  leakage.  The  experimental  phase  agree¬ 
ment  with  Theodorsen  on  Figure  E64  is  about  as  good  as  could  be  hoped  for  with  this 
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type  measurement.  The  various  Theodorsen  predictions  of  Figure  5-21  indicate  that  a 
3-dimensional  version  of  Smuiliri  &  Bender’s  program  might  well  coincide  with  the 
measurements  and  with  Theodorsen  on  Figure  E64. 

Flap  Hinge  Moment  Response 

The  flap  hinge  moment  responses  transmitted  from  NSRDC  through  Grumman 
Aero  Test  are  presented  in  TalDle  E6.  For  thvs  case  it  was  necessary  to  employ  the 
phase  as  reported  to  obtain  a  reasonable  result. 


The  data  supplied  has  had  the  air  calibration  removed  in  accordance  with  the 
relationship: 

C, 
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where,  from  Eq.  (E.32): 
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where,  from  Table  I  of  Reference  El,  C.  is  -.01175/rad.  for  the  low  speed  run 

*'6 

SS 

and  -.0109/rad.  for  the  high-speed  runs. 

The  data  with  and  without  the  air  calibration  removed  are  compared  with  the  final 
result  of  Reference  El  on  Figures  E65  and  E66  which  should  also  be  compared  with 
Figure  23  of  Reference  El.  In  this  case,  again,  a  significant  difference  between  the  air 
calibrations  of  Reference  El  and  of  this  appendix  appears  in  amplitude  but  not  in  phase. 
For  the  flap  hinge  moment,  however,  the  air  calibration  has  a  very  substantial  effect 
upon  the  response  indicating  excessive  friction  in  comparison  with  the  measurement. 

Figures  E67  and  E68  compare  the  flap  hinge  moment  hydrodynamic  responses 
(measured  less  air  calibration)  for  four  runs  with  the  result  of  Reference  El  and  with 
the  Theodorsen  prediction  of  Section  5  of  this  report.  Figures  E67  and  E68  should  be 
compared  with  Figures  5-24  -  5-27.  The  results  are  totally  inconclusive  because  of 
the  persistent  attenuation  at  low  frequencies  which  might  be  due  to  differing  hinge 
frictions  in  air  and  water  (also  indicated  on  Figure  E58  or  to  air  leakage  over  the  flap 
during  the  steady-state  and/or  unsteady  runs.  No  confident  interpretation  cBsi  be  given 
the  unsteady  result  unless  the  air  calibration  and  static  submerged  calibration  can  be 
interpreted  confidently. 

Foil  Pitching  Moment  Response 

The  foil  pitching  moment  responses  transmitted  from  NSRDC  through  Grumman 
Aero  Test  are  presented  in  Table  E7.  For  this  case  it  was  necessary  to  employ  the 
phase  as  reported  to  obtain  a  reasonable  result. 


The  data  supplied  has  had  the  air  calibration  removed  in  Jiccordance  with  the 


relationship: 
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where,  from  Eq.  (E.32): 
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and  where,  from  Table  I  of  Reference  El,  C  is  -.2825/rad.  for  the  low  speed  run 

H 

and  -.275/rad.  for  the  high  speed  runs. 

The  data  with  and  without  the  air  calibration  removed  are  compared  with  the  final 
result  of  Reference  El  on  Figures  E-69  and  E-70  which  should  also  be  compared  with 
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Figure  21  ot  Reference  El  and  with  Figure  £56  of  this  appendix.  For  this  case  the  air 
calibration  has  a  very  substantial  effect  upon  the  hydrodynamic  response  and  the  smooth¬ 
ed  calibrations  of  Referc^c^  £1  and  of  this  note  produce  substantially  different  hydro- 
dynamic  amplitudes  and  phases. 

Figures  £71  aiKi  £72  compare  the  foil  pitching  moment  hydrodynamic  responses 
(measured  less  i\ir  calibations)  for  four  runs  with  the  result  of  Reference  El  and  with 
the  Theodorsen  prediction  of  Section  5  of  this  report.  Figures  £71  and  £72  should  be 
compared  with  Figures  5-28  and  5-29.  The  results  are  not  contributory  in  any  sense 
and  only  present  a  scatter  about  a  steady-state  value  which  has  been  offset  by  random 
friction  values  (in  calibration  and  data  runs)  and/or  by  air  leakage.  The  very  poor 
results  of  Figures  £71  and  £72  give  some  indication  that  the  foil  pitching  moment  re¬ 
sponse  to  flap  motion  presents  the  most  difficult  measurement  problem  of  those  pre¬ 
sented  though  it  is  also  possible  that  the  phase  angle  confusion  throughout  has  produced 
an  erroneous  final  result. 
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NSKDC  Phase  Has  Been  Inverted  (0  -  18(F  ) 
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Figure  E  l.  Fr;  ie  Scaled  Test  Results,  l4.  y  ixiivois,  Eost  sandpaper 


Figure  E^.  Froude  Scaled  Test  Results,  14.  7  Knots,  No  Sandpaper 


Fi^rc  E3.  Froude  Scaled  Test  Results 
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Figure  E5.  Cavitation  Scaled  Test  Results.  27.  5  Knots 
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Figure  E6.  Cavitation  Scaled  Test  Results 
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vication  Scaled  Test  Results,  36.  5  Knots, 


Figure  E8.  Cavitation  Scaled  Test  Results,  40.  9  Knots 
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Figure  ElO.  Cavitation  Scaled  Test  Results,  32  Knots, 
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Figure  E12.  Cavitation  Scaled  Test  Results 
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Figure  E15.  Cavitation  Scaled  Test  Results,  44.  9  Knots, 


Figure  Elt,  Cavifation  Scaled  Teet  Results,  32  Knots,  a  =1.50 
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igure  E19.  Cavitation  Scaled  Test  Results,  32  Knots 
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Fi^re  E21.  High  Aspect  Ratio  Foil  Flap  Cavitation 
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Figure  E22.  Flap  Effectiveness 
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Figure  E23.  Residual  Lift 
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Figure  E26.  Residual  Flap  Hinge  Moment 
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Figure  E32.  Flap  Hinge  Moment,  or  =  1. 5  ,  i  =  4. 5 


d's^  fOf^r,  p^Torm^ 
yppcn  c^f- 
Lomn  Cif-.  .009*^  . 


.(if  itPTe^ct->fiPfnPvsp*ttp  ^Jt^ae. 
Ft^ 


o  WO^AMt  «»i.  APpeAFPHCe, 

tk-  Vi-%  Cfi^nATfOH .. 

^  <rtc  /tw  j 

J-OfKB^F^Adt-^tTArfW . 

V  VH  e^iTjAr/C^j  .__.  .■  .  .,[ 

jPMW  cjMwtndf  sfMiiBt.  itdKATes  ii6nrai  ptfPifMtrrsM 
X  K^Mfftfus  6f  FPnmrtw  pmtt  eiPAecr 
Ja^kt-AT y.lM^_M>Of*ieJiAja/fM  ^k«*wrv/r<r/r/. „. 


M.  lA  .  A  J; 


o.too. 


Km 


D  2£;«:>  4C<7  <^QO  atp  \tiOO  IZJOO  1400  I6CD  1300  dG)0 
FOIL  FOAaiH4fyW/3  i; 


Figure  E37.  Cavitation  Characteristics,  o'  =  0 
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Figure  E38.  Cavitation  Characteristics,  cy  =  -1.  5 
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Figure  E40.  Cavitation  C'naracteristics,  a  =  1.  5 
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Figure  E41.  Cavitation  Characteristics,  cy  =  3® 
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Cavitation  Drag  Effect,  44.9  Knots,  a  =  0,  AG(EH)  Fwd.  Foil 
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Figure  E44.  Flap  System  Inertia  In  Air 


Figure  E45.  Flap  Hinge  Friction  In  Air 
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Figure  E50.  Drag  Amplitude  In  Air,  As  Function  of  « 
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Figure  E51.  Drag  Amplitude  In  Air 
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Figure  E52.  Drag  Phase  In  Air 


Figure  E53.  Foil  Drag  Air  Calibration 
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Figure  K57.  Foil  Hinge  Mor.iant  Air  Calibration,  Before  Freelng-Up  Hinge 
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Figure  ESS.  Static  Calibration,  Flap  Hinge  Moment 
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Figure  E60.  Static  Calibration,  Lift 
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Figure  E61.  Calibration  Effect,  Flap  Lift  Amplitude 
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Figure  E63.  Flap  Lift  Amplitude 
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Figure  E64.  Flap  Lift  Phase 
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Figure  E65.  Calibration  Effect,  Flap  Hinge  Moment  Amplitude 
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Figure  E66.  Calibration  Effeot,  Flap  Hinge  Moment  Phase 
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Figure  E67.  Flap  Hinge  Moment  Amplitude 
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Figure  E68.  Flap  Hinge  Moment  Phase 
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Figure  E69.  Calibration  Effect,  Foil  Pitching  Moment  Amplitude 
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Figure  E70.  Calibration  Effect,  Foil  Pitching  Moment  Phase 


Foil  Pi tc bins  Moment  Amplitude 


. . —  ■ 

;b  t6s  '  /.W^ I  FtveiKia-yrifmec 

^ .  zo.r .  s/j.  9jr  j  .  :  ; 

C  tos  X^frfgff  pffpgifN^.tfP  WH^e 

□  ;+.2.  u*iit7} 

-  '  MF, ie./> Vcj/^.2  ^ io.ffi>»Tij^ix^ FFetj>M6-iF' m/^ae. 

. - imbocmeH 


--aE'i— 


1  ^  ' 

^ 

;  ^  ^ 

RIH 

1  .  .  . 

a 

□  ;  ■  o 

a - - - ^ - 

